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SUMMARY 


This report describes the theoretical development and application of an 
analysis for predicting the major static and rotary stability derivatives for 
a complete airplane. The analysis utilizes potential flow theory to compute 
the surface flow fields and pressures on any configuration that can be syn- 
thesized from arbitrary lifting bodies and nonplanar thick lifting panels. 
The pressures are integrated to obtain section and total configuration loads 
and moments due side slip, angle of attack, pitching motion, rolling motion, 
yawing motion, and control surface deflection. Subcritical compressibility 
is accotmted for by means of the Gothert similarity rule. 
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Within the scope of predicting the total configuration stability 
derivatives it was necessary to study the problem of computing the spanwise 
variation of potential form drag due to panel lift and thickness. Included 
in appendix F is a solution to this problem. Also, in solving the potential 
form drag probiem the work of Woodward and Wagner was thoroughly analyzed. A 
complete derivation of Woodard's influence equations is given in appendix D 
and a comprehensive review of Wagner's lifting surface theory is included in 
appendix E, 
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INTRODUCTION 


In order to develop a procedure for predicting total configuration 
ctability derivatives, the perturbation flow due to the vehicle had to be 


represented by a grid of efficient aerodynamic finite elements general enough 


to satisfy the boundary conditions over the surfaces of diverse shapes and 
also nroduce the correct resultant forces and moments. the quadrilateral 
vortex was selected to represent the perturbation velocity due to the bodies 
because of its numerical efficiency. net force producing capability, limited 
raige of influence,and reiaiionship to. the horseshoe vortex which has 


demonstrated amazing accuracy in predicting loads on wings of arbitrary shape. 


The nonplanar thick lifting panels are divided into two sections, (1) 
she outboard section which is defined by a locus of chord lines, and (2) the 
soot section which is a transition. region from the outboard section to the 
juncture of the panel and a body. If the panel is attached to another panel 
share is no root section. The perturbation velocity due to panel lift is 
represented by quadrilateral vortices in the root section and skewed horse- 
choe vortices in the outboard section. The perturbation velocity due to 
nanel thickness is represented by a source lattice. The panel aerodynamic 
finite elements were selected because of their numerical efficiency and 
proven accuracy in predicting flow fields over wings of general shape. 


The panel singularities are placed on a mean surface instead of the 
actual external surface of the panel to maintain computing efficiency. This 
will sacrifice surface pressure accuracy at the juncture between two panels 
or a panel and a body, but for this type of general analysis the savings in 
computer time makes the compromise practical. Second order corrections to 
account for the interference between lift and thickness and to account for 
blunt leading edge airfoil sections are included. 


The source and vortex lattice influence equations are formulated in 
avis of the same quantities, which allows the perturbation velocity due to 
-€t and thickness to be computed simultaneously and thereby save computing 
effort. Also, due to the limited range of significant influence of the 
quadrilateral vortex the influence of any quadrilateral vortex is only 
conmuted at those points within a given area of influence. This can save 
considerable computing time in developing the aerodynamic influence martix. 
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Computer time is also saved by reducing the number of unknowns by 
transforming the aerodynamic influence matrix by constraint matrices. The 
-onstraint matrices constrain the body and panel vorticity, thereby reducing 
the number of unknowns from that of the number of vortex elements to the 


number of constraint functions, This is an option in the program and can be 


applied in just the longitudinal direction, lateral direction, both directions, 
or not at all. 





The pressures are integrated by numerical means and the panel section 
drag is computed by means of the Kutta-Joukowsky theorem. The total 
configuration induced drag is computed by means of a Trefftz plane analysis. 
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tH ORETICAL DEVELOPMENT 


 Confisueation Representation 








~The theory discussed in this report is capable of predicting the surface 
pressures and integrated loads on any configuration which can be synthesized 
from lifting bodies and thick lifting panels of arbitrary shape. These two 

















basic elements can be attached along longitudinal subpanel edges in order to 


represent complete airplane confisurations of arbitrary shape. 





_ Lifting Body. - The arbitrary lifting body, as shown in figure 1, can be 
of the solid or flow-through type. The body external surface is divided into a 
grid, which represents the edges of the body subpanels. The body bound vortex 
lines are placed at the quarter chord point of the subpanels and the fixed 


_ trailing vortex lines along the longitudinal edges of the subpanels. The 
_ trailing vortex lines are shed from the aft end of the body and extend in the 





X direction to infinity unless the free wake option of the program is used, in 
which case the locations of the trailing vortices aft of the body are deter- 
mired such that they are force-free. If the body is closed to a point at the. 
aft end, all of the body trailing vortices aft of the body cancel and are, 


therefore, neglected, 
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Figure 1.- Arbitrary lifting body. ‘OF POOR QUALITY: 


The body vortex strengths are assumed constant around closed paths defined 


by the bound vortex lines of two subpanels, adjacent to each other in the longi - 


tudinal direction, and that portion of the longitudinal edges of these adjacent 
panels between the two bound vortex lines. The contribution from the vortex 
loop, or quadrilateral vortex, is defined as positive when the Biot Savart line 
integral is taken in the clockwise direction by an observer looking at the 
external surfaces of the two subpanels. The body subpanels and quadrilateral 
vortices are shown in figure 1 along with the location of the body control 
points. 


The body control points are located at the three-quarter chord of the body 
subpanels. At these points the total flow is summed and forced to be a minimum 
in the direction normal to the body surface. If the body vorticity is not con- 
strained by functions with unknown coefficients, a control point is placed at 
each subpanel and a discrete solution for the unknown body vortex strengths is 
obtained. If the body vorticity is constrained, control points are placed at 
as many subpanels as is necessary to obtain a good representation of the body 
shape and to insure that there is sufficient control of the constraint func- 
tions. In this case the unknown coefficients of the constraint functions are 
determined by the method of least squares. | 


The body subpanels are further subdivided in the lateral direction so that 
the vortex grid is mapped to the body surface more accurately. This also 
allows the use of only one quadrilateral vortex in the lateral direction for 
the case of a body of revolution in a uniform flow at zero angle of attack. In 
this case the two side edges of the subpanel are coincident, and therefore, the 
fixed trailing vortices cancel leaving a longitudinal distribution of ring 
vortices located at the quarter chord of each subpanel. 





Figure 2,- Body subpanel defined by points 1, 4, 5, and 8 with 
Subdivision points 2, 3, 6, and 7. 
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The coordinates cf the body subpanel subdivisions, as shown in figure 2, 


arc used to compute unit vectors tangent to the body surface in both the longi- 
tiidinal and lateral directions, unit vectors normal to the body surface, 
directed subareas to he used in the integration of the body surface pressures, 
and the centroids of the subareas to compute moments about the configuration 
center of gravity. 


> The unit vector tangent to the equivalent incompressible body subpanel in 
the longitudinal direction is given by, | 3 
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-- The unit vector tangent to the equivalent incompressible body subpanel in 


the lateral direction is given by, 


f47 G+T, k (5) 
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(%3-%2 } BL Y5-¥,)*8'( 25-2, ) + (XX) + B(Y6-¥>) x 6 (25-24)° 
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The unit, vector normal to the equivalent incompressible body subpanel is 
given by, 7 


N = Ny i +N, 5 +N, eee (9) 





The directed subareas for the middle lateral subdivision for the actual 
bedy are given by; 
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eee - --~ *** The components of the position vector to the centroid of this subdivision 
—_ 2s a 6 4S 
_ 7: = ey + Y a ee (17) 
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Se Analogous directed subarea and centroidal position vector expressions are 
- obtained for the other lateral subdivisions. 
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— 7s The body description is defined independent of the subpanel grid and can be 
===="input to the program in a number of different ways. Both cartesian and polar 











7 ~ Géordinates” can be used to describe the body cross-sections at a set of chord ata 
| stations. After the basic cross-sections have been described they can be 








scaled and then translated in planes perpendicular to the body chord or mean 
camper line. This provides a means of inputting the correct side and top views 
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of an arbitrary body with a minimum of input data to obtain a preliminary 
estimate of the loads on the body. It also facilitates in inputting the exact 
_ shape of some bodies which can be represented by longitudinal segments, over 
which the cross-sections are mathematically similar, The actual arbitrary 
shape can be input directly without the use of the scaling and translation 
options if it is more convenient. 


If polar coordinates are used to describe the body cross-sections, the 
lateral location of body subpanel side edges or fixed trailing vortex lines 
are also given in terms of angles measured from the local section Zp axis, in 
a plane parallel to the (Yp-Z,) plane, at an independent set of chordwise sta- 
tions. These subpanel lateral edges are either specified at a given set of 
angles or defined to be at equally spaced angle locations. The subpanel longi - 
tudinal edges are specified at given Xp stations, evenly spaced in tems of 
Xp, or evenly spaced in terms of p, where $p = cos-1 (1 - 2 Xp/Cp). 


The following procedure is used to determine the coordinates XBps 
(YBe - AYp.)/ (YB) » (ZB, ~ AZB.)/ (ZB) of the subpanel corners. 
1, Rp. versus 6B, is input at Xp. - 
2. OB, (subpanel lateral edge location) versus Xp g iS given as input 


5. Interpolation on 6p, versus Xp. is done to obtain OBr versus XBp: 
(subpanel longitudinal edge location) and XB. 


4. Interpolation on Rp; versus OB; is done at each Xp, to obtain Rp 
versus Xp; at each 6 Br 


5. Interpolation on Rp versus XB. is done to obtain RBe versus XBe at 
each 6Br. 


6. XBr; (Ypp- AYp,) / (YByp » (2B, A2p,,) / (ZB) are computed from XBr: RBe» 
and OBE: 


If the subpanels are subdivided the coordinates for the corners of these sub- 
division are computed using the same procedure. There is always an odd number 
of subdivisions in both the longitudinal and lateral directions. 
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If cartesian coordinates are used to describe the body cross-sections, 
the lateral location of thc subpane] side edges are defined by the percent 
circumferential length Ni. = SB, Ss PMN, at independent chordwise stations. 
Tables of (YR. 


iS AY). j/Y¥p , and (25. 
= Sp. «/ Suny. arc de re at the input longitudinal stations Xp, - 


.- OZ. ) 2, versus percent circumferential 


The | same " procedure that was used to obtain BE OBE AY'3E)/V By and 
(2 BE” AzZp,)/ 2B W when the body was def ined by polar coordinates is also used 





for this case, “except that opis is replaced by (p,- AY. YB and 
(2p.- Alp, ,)/ 2 By? and ie is roplaced by Tpi SB, i! SBunx, . Also, step 6 


umecessary, The procedure is s eycled through twice, first for 


After XBe py %5,) Nay ana d (2p, 355) ave been computed, Yg,. 


= “is determined by multiplying (BET Yp/ Ypyy by the multiplication factor YBa 











and ‘then adding the translation increment AYRE | ZBr is determined by 























“ wultiplying (ZB, AZp.)/ a by the multiplication factor 2B and then adding 
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| length, 


the translation increment AZRp: 


Thick lifting panel. - The thick lifting panel, as shown in figure 3, can 
be warped in any manner laterally, have an arbitrary distribution of nord 
thickness, twist, and camber. It can be used to represent a wing, 
canard, Fin, pylon, horizontal or vertical tail, or be wrapped around to rep- 
resent a flow through nacelle. The panels can be attached to other panels, 
The panels 


can have plain leading or trailing edge flaps, ailerons, rudders, or elevators. 











such as in the case of a pylon on a wing, or attached to bodies. 


‘These control surfaces can be of the full or partial span type and their hinge 





lines are not restricted to ies een chord lines. 


The panel is siyided 4 into two ; Sections; a the outboard section where 
the subpanel ‘longitudinal edges are assumed to be straight lines in the X 








-- “direction, and (2) the root section which is a transition region from the 
| ‘outboard section to the intersection ‘of the panel ‘and a body. 





24 ‘The line of © 
intersection between the panel and a | body does not have to be a straight — a 
dine. | Therefore the subpanel longitudinal edge lines will change in shape 

















_ = from that of the line of intersection at the side of the body to a straight 





Jine in the X direction at the outboard section. If the panel is not 


| attached to a body, there is no root section. 


“~~~ SOURCE LATTICE 
VORTEX LATTICE 





Figure 3.- Thick lifting panel. 
any distribution of sub- 
are specified as a 


The root section and outboard section can have 
The subpanel lateral edges 
evenly spaced in terms of percent chord, or 


panel lengths and widths. 
list of percent chord stations, 
evenly spaced in terms of ¢, where 
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The subpanel side edges are specified as a list of percent of surface 
percent of surface semi-span, or 
= cos7i%, 


| Semi-span 1, evenly spaced in terms of 
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evenly spaced in terms of 8, where 0 
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The percent of surface semi-span 7 is defined as the percent of length 
of the line projected into the (y-Z) plane by the pane] leading edge. A 
constant 7 line is in general curved in the root section. In the outboard 
section, lines of constant 7 arc straight and in the X direction. The curved 
constant 7 lines associated with the longitudinal edges of subpanels in the © 
root section are computed such that the sweep of the leading and trailing edges 





‘of the subpanels vary linearly from that of the sweep of the panel at the lead- 


ing edge to that of the sweep of the panel at the trailing edge. Also, the 
comer points of the subpanels are equally spaced in the lateral direction | 
along lines of constant percent chord. The chord at a lateral station in the 
root section is the length of the curved constant 7 line at that station. 


If the lateral distance along a constant percent chord line between two corner 


- points is defined by |AS|, the unit vector tangent to the leading or trail- 


ing edge of a subpanel at any percent chord station is given by; 
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and Ky equals the number of subpanels in the lateral direction in the root — 
section, | 


C5 = Gy - i+ ,-y) f+ @,- 2k (21) 


where (Xp, Yr, 2p) and (Xj, Yj, Z j) are the points on the curved constant per- 
cent chord line at the line of intersection of the panel and the body and at 
the juncture of the root and outboard sections, respectively. The equally 
spaced corner points along a constant percent chord line are then computed by 





K=Ky 
R XYZ) qq & R XYZ) + >, ar) | 48, (22) 
| Kel K 


where Kn is the number of subpanels, in the lateral direction, the constant 
7 line is located from the line of intersection of the panel and the body. 
Since constant 1 lines in the outboard section are straight and in the X 
direction the subpanel corner points are defined directly from the input list 
of subpanel edge locations in terms of percent of chord and lateral y station. 


The coordinates of the leading edge point (XL.e.> YL.E.» 2L.8, J associ- 
ated with a constant 1,- line are obtained from the panel perimeter descrip- 
tion by converting the leading edge input points to tables of XL .E.4 versus 
n » YL.E.j versus 7 , and 41,,E,4 versus n , and then interpolating in these 
tables to determine the leading edge coordinates at a desired lateral n. sta- 
tion. Unit vectors tangent and normal to the panel at any Ny station can 
_ also be determined from these tables by evaluating the AYL,E., and AZp pF. 








about the lateral station nx. The unit vectors are then defintd by; : 
Piha ; Piles ; 
Np = - j + k = (23) 
K 2 Z 2 2 
AY + AZ AY + AZ 
L.E., LE. L.E., L.E., 
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- -edterated for such that the wake is force free. 








ie ee Lee: 
T= at ja s i (24) 
, - A 
oe K med AY, . rs A a". AY, Z + A zy 
sea pe ae. wo" K ~K __K 

















he the ese corner ‘points ‘mown, the coordinates 

and seep ‘of the vortex and source lattices used to represent the perturbation 

~ velocities due to lift and thickness, respectively, can be defined. The bound 

vortex lines and control points are placed at the quarter and three-quarter 

chord roe of the subpanels, re The fixed trailing vortices are. 
Ac In the root section the vortex lattice 

is a Gunteiaterai system, where as, in the outboard section the vortex lattice 
is a skewed horseshoe system. The skewed source lines are placed at both the 


_ quarter and three- “quarter chord points of the subpanels. 


With oe ‘weorlinaces 5 of | 





























The” section geometry ae ‘input in terms of a percent thickness, percent 
camber, and twist. Where the percent thickness and camber are based on the 
rie chord and measured in the Np direction. Twist is defined in the 
7 lane described by the unit vectors i and NPK- The X component of the 
unit vector normal to the section mean camber line and in the plane defined 


by — a cous NK is: meee piven oe 

















(25) 


where Zc is the perpendicular distance between the section chord line and 
the mean camber line, e€ is the angle of twist, and 6 the deflection of any 
“= eontrol surface. All three of these quantities are functions of both percent 
chord and ae The thickness is defined in the same manner as camber. 





| =: The “‘teailing vortices aft of the trailing edge of the outboard section 
are straight lines in the X direction going off to infinity. The trailing 
vortices aft of the trailing edge of the root section lie along curved con- 
stant lines to the end of the body and then go off to infinity in the X 
-.direction. These curved constant lines are determined in the same manner 


as th If the free wake option of the 





~ as the constant n lines in the root section, 
prosram is utilized the location of the panel free trailing vortices are 
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Discrete Influence Equations 


The perturbation velocity due to the arbitrary lifting bodies is 
represented by quadrilateral vortices on the external surface of the bodies. 
The perturbation velocity due to lift on the panels is represented by quadri- 
lateral vortices in the root section of the panel and by skewed horseshoe vor- 
tices in the outboard section of the panel. The thickness is represented by 
skewed source lines in both the root and outboard sections of the panel. All 
the panel singularities are placed on the panel chordal surface. The source 
strengths Z are defined by the change in thickness over that portion of the 
subpanel it represents, so that: | 


2BAZ V. 
z t X 
~ = eee (26) 


V1+ (tan ajyp2 Ve 





where <A is the sweep of the source line and Vy is the total onset velocity 
in the X direction. The quadrilateral vortex strengths K and the skewed 
horseshoe vortex strengths I must be solved for utilizing the boundary condi- 
tion that a minimm of flow passes through the external surface of the bodies 
and the chordal surface of the thick lifting panels at a finite numberof contro] 
points. In order to satisfy this boundary condition the total flow due to all 
Singularities and onset flow is summed at each control point and the scalar 
product of this sum and the surface uit normal is minimized. This results in 
a set of linear aerodynamic influence equations which are solved for the 
unknown vortex strengths by means of Householder's method, described in 

' appendix A, 3 


: .th ; 
The influence cquations for the j~ equivalent inconpressible body are 
given by3 ‘| 
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» and N, are the components of the Hed equivalent 


B. | B, a ae body surface unit normal vectors, 
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er aes a 
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= 4s | are the components of the perturbation veloc- 
Se 7 , AY » and Ay ity induced by the jth equivalent ee a aes 
ng is “BB, BB, Bri body unit strength vortices onto the j 
ie equivalent incompressible body, 





e* p » Ay p jane » a p are the components of the perturbation ese 
ane B. K j K induced by the kth equivalent incompressible 
panne — _ -pariel unit strength vortices onto the jt 


[Sear 7 oe eS - equivalent incompressible body, and 
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S ; Sy » and So are the components of the perturbation velocity oo. 
"B, Py BSP B.P, induced by the kth equivalent incompressible { 


panel unit strength sources onto the jth eel 
equivalent incompressible body. 


i) 


Vy Vy Vo 
cra ee i ae I and vy { ate the components of the onset flow divided 
” by the reference velotity V_, 


{ K} are the quadrilateral vortex strengths on the te equivalent 
i incompressible body, and 


strengths, the outboard section skewed vortex 
strengths, and the skewed source strengths on 
the kth equivalent incompressible panel, res- 
pectively. The vector {e}p. is the flow through 
the surface of the jth equivalent incompressible 
body at the control points, for a discrete body 
solution this vector is zero. 


1K}, Ep , and {2}, are the root section quadrilateral vortex 
K K K 


The influence equations for the jth equivalent incompressible panel are 
given by; 


i 





\N A S \N K “— 
» Xp pe "ip 2 * 1 2p hee ce 
j ji j ji j ji ol 
i=] \ 
N K 
~ INS TPA + fy TPA + Pe [4 
r > Py P, PPL Pt) PP), 
\ . 2% 
K=1 V P 
k 
N 
p = 
S \N, S \N S . 
x + Y + Z Z —— 
A Xp. Pp Pp. Pp p PP INV 
2 1 jk j jk j jk | a 
=} | \ \ k 
~ V. N V V 
X Y N z)_ 
% ref 6, ir °3 ie "| 
. j 
(28) 
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are the components of the j~ equivalent incompressible 
panel mean camber surface unit normal vectors, 


“are the components of the perturbation velocity 


ae .  nduced by the ith equivalent incompressible 


body unit strength vortices onto the jth 
equivalent incompressible panel, 





arc the components of the perturbation = 
velocity induced by the k™ equivalent incom- 
pressible panel unit strength vortices onto 
the jth equivalent incompressible panel, and 


are the components of perturbation velocity 
induced by the kth equivalent incompressible 
panel unit strength sources onto the j 
equivalent incompressible panel. The vector 
{el,; is the flow through the mean camber 

PJ ; 
surface of the jth equivalent incompressible 
panel at the control points. This vector 


Boreas is also zero for a discrete panel solution. 


<The onset flow velocity ratios 
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_. =<, > and~—are given by the following expressions, 
Vas Ve Mee 


Ae Gy 2B (Y-Y. ,) 
of Se C.G. _ 


8 
a | 


b 


¥ 2B (2-2. 2 (XX. g 


BBY eae ne 5 (30) 


) 2 (X-X ) 


2B (Y-Y 
C.G." , gs G. 633) 


ees ee 


if 


ic 


Where P*, q*, and y* are the nondimensional roll, pitch, and yaw rates, 
respectively. These are defined as; , | 
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ZV 2V ZV 


oo 


Pe = P/ Ts aR = g/ =", and y* = ¥/ 5 (32) 


The above onset flow equations assume that angle of attack a and angle of 
yaw # are small. The coordinates Xc¢ ¢ , Yc.g,, and Zc.q. define the location 
of the center of gravity. 


Equations (27) and (28) can be combined into a Single matrix aerodynamic 
influence equation. For a completely discrete type solution the influence 
equation is, : 


K 
[Alp g [Als p i>] Is] 
ji Jk ; j 
| 1 
KL ; 
Mie 
19 
‘Alp ap. Alp p Ve {5p on 
i j 90 j 
Py 
where 
B. BL» Boy s+ oy By 
Be = By, Bo, a Ee ee N 
P.=P,P.,...,P 
j 1 ia Np 
PL. = Py Pos So Secciacy Pa 
. P 
NB = Number of bodies 
and 
N,, = Number of panels 
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The matrices 


[A] 


[al 


B.B.’ Bp 


are defined as; 


A (34) 
ay 
7 1 

A (35) 

Z 

N A (36) 
| fea 
j Ja 

; (37) 


The known quantities in the influence equation {¥V S| 2. and IVS}p, are defined 
as; 





Ps] [al DS] 


car Beet ; : 
| \vs| 27 > 
Ke 





tes - = a = RK. (ras * = 
n Wl LR | Jz t Ga 
as V B. V 
J © J ? 
ie 




















Vs 
k-] 
JN j k = Py j Veo 
Jd | Veo Jy | Vo 














| The elements of the above matrices are computed using the influence 

equations derived in Appendices B and C. Each element is associated with the 
influence of a singularity on a control point. The Singularities and control 
points are ordered such that all of the longitudinal stations for the first 
lateral station are cycled through first and then all of the longitudinal sta- 
tions for the second lateral station. This process is continued until all 
Stations have been cycled through. The longitudinal stations start at the 
leading edge of the panel, the nose of a solid body, and the tail end of the 
inside surface of a flow through body. The lateral stations start at the 
inside edge of the panel and go toward the tip of the panel. The lateral 
Stations on the body start at the top of the body and progress in a clockwise 
direction when looking at the body from the tail to the nose. 


All bodies are cycled through first and then the panels. Both the bodies 
and the panels are cycled in the order that they are input. The colums of 
the influence matrices are associated with singularities and the rows with 
control points. The elements of the matrices and the submatrices of the com- 
bined matrix influence equation are sequenced, for both the columns and the 
rows, in the same order that the singularities, points, bodies, and panels 
are cycled. 


The influences of quadrilateral vortices in the panel root section and on 
the body surface are computed by equations derived in Appendix C. The influ- 
ences of skewed horseshoe vortices on the panel outboard section are com- 
puted by equations derived in Appendix B. The influences of skewed source 
lines on the panel are computed by equations derived in Appendix B. 


If the force free wake option of the program is used the influences of 
the free trailing vortices are computed using the influence equations in 
Appendix C. The locations of the force free trailing vortices are computed 
using the following iteration procedure. 
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A solution for the vortex strengths will be obtained first by 
assuming the location of the free vortices to be in the longitudinal 
direction and to be straight except in a thick lifting panel-body 


the total velocities ti; ¥) and Win the X,Y, and = directions, ~~ 
respectively, will i: computed at the midpoint of each free trailing 


vortex division. a 





The actual mean camber surfaces of the thick lifting panels will be 


computed so that thc correct telationship between fixed vortices and 
control points on the thick lifting panels and the free trailing 
vortices is maintained. | nae 


aa Bu ce cae blanca: 
~ = 2 ug TL SD , a % : 


The ratios V/U and W/U are integrated in the X direction from the 
edges of the actual thick lifting panel mean camber surfaces and body 
aft ends where the free vortices are assumed to be shed, in order to 


obtain their new locations. 


The influence of the free vortices on the actual thick lifting panel 
mean camber and body surfaces is computed and the difference between 
this influence and that from the free vortices at their previous 


location is added to the influence matrices for the complete 





- configuration. 








A new solution for the vortex strengths is determined and a new set 
of total velocities U, V, and W along the new free vortex lines is 


__computed using quadrilateral vortices and source lines on the actual 








thick lifting panel mean camber surfaces and quadrilateral vortices 
_. on the body surfaces» 





‘The above procedure is iterated between steps (4) and (6) until the surface 


When the influence matrix of any individual body or 





“pressures converge. 


panel is no longer significantly changed due to a new positioning of the free 


The vortex strengths on any of the bodies or panels c 
representing the vort 
unknowns, equations, 


_ yortices, the calculation of the perturbation to that matrix is terminated. 


Constrained Influence Equations 


an be constrained by 
icity by a finite series. This uncouples the number of 
and vortices used to represent the perturbation velocity. 


The series used must be capable of producing the type of perturbation velocity 


ber of equations or 


needed to satisfy the boundary conditions at the control points to an accept- 
able degree. If such a series can be found for a given body or panel the num- 


control points can be reduced from that of the number of 


Nae 


a 


vortices or subpenels to that. necessary to describe the geometry. This sub- 
stantially reduces the computer time, since in general the number of vortices 
needed to represent the perturbation velocity is far more than the number of 
control points needed to represent the geometry of the body or panel. 


This reduction in the number of required control points results 
primarily in reducing the computational effort necessary to set up the aero- 
dynamic influence equations. A further reduction in computer time can be 
realized in the solution of the influence equations, if constraint functions 
are used, since the number of unknowns is reduced from that of the number 
of vortices to the number of terms in the constraint series. This results in 
a system of equations which is overdetermined and is solved by the method of 
least squares. In this process the vectors {e}p. and te}p, are minimized. 
The body constraint equation is given by;, J 


K _| T a. 
{s]t . 
1 : 


where 


T 7 t, ye _ | 
B. |=] R ad. Cae (41) 


The matrix [R] is a transformation matrix used to obtain the quadrilateral 
vortex strengths K from the bound horseshoe vortex strengths ry. 


R |=] Il, (42) 
% 


where [R]yn is made up of ones along the diagonal and in the lower triangle. 
The upper triangle is filled with zeros. 


The indices m and n refer to longitudinal and lateral subpanel locations 
on the equivalent incompressible body, respectively. The vorticity ratio 
(V/Vio)um is the value of the vorticity series at the mth subpanel from the 
nose of the body and the nth subpanel in the lateral direction around the body, 
The prpen er er length ¢ Mon is given by; 


| ‘J; — 
Hi 2 a co 


(43) 
mn Mom-1)n m M m+1)n 


| w 
|e 
ro| 


20 


CG 


if 

















Iii 
it 
i 


Fa 
ibi 
io 


\ 


ral 


1 | 

| 
Mg 

nf: iii 





| 
it 





I 
} 


! 


| 














— 


| 
I 
| 

1 


wes 
bite 
My fi 





i 
i 





l 
f 





: 
‘RI! 
BUA! | 











wt | 
i 


where L, is the length of the mth subpanel from the nose of the equivalent 


M 
mn. th : 
incompressible body at the n” lateral station. 
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The vorticity series (V/s3,,, is the product of a icngitudinal series and 
a lateral series. The cloweits: of the matrix ap | are the unknown coeffi- 
ene oe ‘ nae PSG. wae, is ; 1 

cients associated with the teius in (We), produced by the product of the 


‘longitudinal and lateral scrics. The terms in (V/Vo),, are ordered such that 


mn 
the products of all of the longitsdinal constraint functions and the first 


lateral constraint function are first, then the products of all of the longi- 
tudinal constraint fimetions and the second lateral constraint Functions are 
second, and etcetera. This process is continued until all combinations of 
lorgitudinal and lateral constraint functions have been cycled through. 


Both the longitudinal and lateral constraint functions are defined over 
seyments of the body. The same functions are used in all segments. The 





origin of the segment is designated by the subscript o and the end by f. 


Tre longitudinal and lateral constraint function segments are given in the 
data input array. 


The longitudinal constraint functions for the body are defined as; 


Vo O/C), = 1 2 
i XB ‘a 


Ty 2 1/2 


Y | 2 2 





vy : 
vy. (pz = COT %/, (46) 
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bs | —— 
7° (m), = COT (7/2 - p79) (47) 


o,- 9» 
of 
(Oe OIN ‘(Gee : =) (48) 
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(59) 


(31) 


(52) 


(53) 


(54) 


(55) 
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above sequence. | 


The fuactions to he used on the body are designated by their number in the 
The last tree functions listed above are special functions 


acid are input at the focetions in the data array, designated by their suh- 
The indepeirlent variasc ty, 14 2 1Vven. By, o, = Con? iL. = 2{X,/C,) | 


Se cipts. 


whe lateral constraint fimetions For the body are defined as; 


es Fad] 73 (87) 


ae Nv 7 |Nz 4 Ny (59) 


= Sin | 7 (60) 
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| 6,- 8 
| Y en ak Bo 
| Vo 6p), = Sin | 2m G 9 (62) 
0 
: 6@.- 86 
| -Y 2 {32 
! Ve, (65). cos |2m a. 3 (63) 
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Here again, the lateral functions to be used on the body are designated by 
their number in the above sequence. If no constraint functions are specified 
the program will do a discrete solution. 


The panel constraint equation is given by; 


(65) 


tl 
~~ 
Ae) 


The | Tpy| transformation matrix condenses all rows of the discrete aerodynamic 
influence matrix by the following procedure. 


1. That portion of the row which is associated with the vortices on 
panel Py is divided into the elements due to each lateral station. A 
hew temporary matrix j|Ay| is developed with all of the elements due 
to the votrices at the first lateral station in row one, all of the 
elements due to the vortices. at the second lateral station in 


row two, and etcetera. 
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2. This temporary matrix is then premultiplied by the spanwise constraint 
function transformation matrix [T(n)] and postmultiplied by the 
chordwise constraint function transformation matrix [T(NC) J. 
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the ‘trangfored Lenporary matelx Tp] = [r¢n) Jag Irov0)1 is che 








~ opened up into a row again and replaces the old row in the original 

- discrete aerodynamic influence matrix. The new row is formed by _ 
placing all of the elements from the first row of {A_] in tha por- 
tion of the row, from the discrete matrix, due to the P, panel first, 
then by placing all of the elements from the second row of [A] 
next to those elements from the first row of [A,], and etcetera. The 
result will be the reduction of the number of elements in a given row 
due to the P, panel from that of the number of vortices or subpanels 
on panel Py to the number of chordwise constant functions times the 
number of spanwise constraint functions. 
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- only the chordwise constraint functions are used, then the number of ele- 
cuits in the row due to the Py panel will be reduced from the number of 

Wit ices on the Py panel to the number of chordwise constraint functions times 
he n uber of subpanels in the lateral direction. On panels with freestream 
































‘pile! t span stations other than the tip, only the chordwise functions should 
ie aes ‘the available spanwise constraint. functions are not sufficient — 


the necessary perturbation velocity to satisfy the boundary condi- 
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Were the indices m, f, k, o, h, and w indicate the number of the subpanel aft 
of the panel leading edge I /V, is defined on, the number of the trailing 
cdge flap, the number of the leading edge flap, the first term of the 

Birnbaum series, the number of the sine term in the Birnbaum series, and 

tte number of the spanwise constraint function, respectively. The quantities 
Nar Ne» Nqy Py, and N, are the number of terms in the Birnbaum series, the 
number of trailing edge flaps with unique hinge line locations, the number 

of leading edge flaps with unique hinge line locations, the number of special 
spanwise constraint functions, and the total number of spanwise constraint 
functions, respectively. 


The columns of the matrix [T(X/C) Jare made up of the (I /Vo),, values, 
where m indicates the number of the element in the column. There is one col- 
umn for each set of (I’ /Va),,'s associated with a unique chordwise constraint 
function. The [T(X/C)] matrix is filled such that the (I’om/Vq) values are 
in column one, then the (T’Nn/V, ) values are in the next set of N,-1 columns, 
then the (T’.,/V..) values are in the next set of N¢ columns, and then the 
(T'iq/V~ ) values are in the next set of Ny columns. 


The [TX/C) | matrix is solved here in the same manner as is shown in 
appendix G. of reference (26). The (I’/Vq)m values are solved for such that 
the same downwash is obtained at the three quarter chord point of each sub- 
panel, except at the last subpanel, due to a distribution of discrete vortices 
as would be obtained by integrating the vorticity distribution in the Blot- 
Savart integral for each chordwise constraint function. The additional 
condition, that the sum of the discrete vortex strengths must be equal to the 
integral of the vorticity distribution is also used. These conditions result 
in the following matrix equation which is solved for [T(x/c) ]. 


fe] [Tex/c)] = [W] (67) 


The matrix [E] is the two-dimensional discrete vortex influence matrix and 
is defined as follows. 


[E] = (68) 


Oe ee i el 
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Where 0 s (x/c) <1.0 and the indices j and i as used here indicate the loca- 
tion of the three quarter and quarter chord points, respectively. 








There is one column in the [w| matrix for each unique chordwise con- 
streint function. In these columns is the downwash due to evaluating the 
vorticity for each chordwise constraint function in the two-dimensional Biot- 
Savert integral. The columns in [w] are ordered such that the column due to a 
given chordwise constraint function is in the same location in [w] as its 
corresponding ( ’/Vco )m column is in the [T(x/c)] matrix. The five basic 
types of columns in [w] are defined as follows. 


Tue to the first term of the Birnbaum series; 
bee Ee oe ts . | ; 
oo cote/2 = (69) 
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Due to the second term of the Birnbaum series; 
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Due to the third and higher order terms. of the Birnbaum series, 


-cos n , 
-cos n d, 
W = (71) 
sin (n¢) ee os ; 
0 
-} 
where °; = COs E - 2 Gc), | 
Due to the trailing edge flap term; 
Wel 1) 
It ; : ad 


sin 1/2( 9+ %¢) 
Sin 1/2(%- %¢) 


—@a @gaewewee @ fF we = 


log 








1/2 sin Of 


where wel, ) = -(t- @-) for os, < Oe, Wel ®, ?.) = O. for %¢ 95s 
is the an coordinate at the Flap leading Sit and, equals the sritnber 
of subpanels per chord. 
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due te the leading edge flap term; 
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a | | (73) 


Sin 1/2( + $1) é. 
sin 1/2(%- $4) | "k w,-1) 
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log 


1/2 SIN ¢ 








ere ow, )=-(r - @) for 0 § @. <9 aM (9 .) ez for 9, as sm, and o 


is the 1 polar coordinate at the Flap lea ang Kedde. 
tee ce the root section the a X/C) ] matrix is trans fonned to if (X/C)], 
“whare each colum of [T X/C)] is transformed as follows; 








Pe ae TWO, = 2 TO), ae 




















whe eas tadicate the number of the element in the colunm in the [T (X/C) ] 
_ and [T OYC)] meric espectively. 

‘The rows of the (r ii matrix are made up of the spanwise constraint 
functions. Each element in a row is equal to the value of the constraint 
=——— =~ ~~ fimction at the lateral subpanel which has the same lateral index as the 
————— number, of the colum in the [T (1)] matrix. There are rows in the [T (n)] | 
SS matrix, with the standard functions described by ae fe - né in the first 


SS —-—™ fe 


ee _ Pi rows, and then the era Py (us functions in the aast Py TOWS . 






























































=== ae The: us spanwise constraint Functions are needed to account for flow 
—— indiced by discontinuities in the sweep of the constant percent chord lines, 











for body-panel juncture induced flow, and for partial semi-span flaps. There 

_-aré two basic special spanwise constraint functions, (1) polygonal functions 

| uhich account for the discontinuities in the sweep of the constant percent 

_ Chord lines and the body- panel junctures, and (2) flap functions which account 
re ‘the partial semi-span flaps. These special functions are derived in 
Appendix F of Reference (26). 
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a 
| 


If the range of influences inboard and outboard of a discontinuity in the 6 
sweep of the constant percent chord line are defined as Anj and Sno, - 
respectively, and the discontinuity station np, then the following expressions | 
define the special spanwise polygonal constraint function. 


For o< < . 
ny! an, 


(An; a an) (1 = ny) 





ec eee z 1 
P (6), (1 1,/4n,MCe,>) + a ey - casaacly ay | as P (6, 8.) 
1- O 
(75) 
47; * Ano 1 -1 : 
| On, an, ; On; iene rao) ny | P (es cos (ny + 4n9)) 


For Sn, < |m| < (1 - 4n,) 


P(@)o = (1-1, + An.) (Q/an,) P (8, cos” (a, - 4n,)) 


An + An. 
| lot An; 


Bo | | 1-14 P (6, 6,) 


= 
+(l- 4) - Ans) (an) P (0, cos (1, + 4n,)) 


| (76) 
And for (1 - An.) < Int <1.0 


P(¢)," (1-4, + 4n,) (1/ Sn.) P 6, cos” (WA - 47.) 
S b i i b 1 


; ji + G- "yd / ans] P Ce, 64) 


(77) 
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a 
sin 1/2 (6 - 6) 


Rr 
(cos 6 - cos 6) LNG, = 
aie Vee 28) 




















* 
cos 1/2 (6 + 6) 


i” | 
+ (cos 6 + cos 6) LOG ; 
cos 1/2 (@ - 6) 








& * ; R . 
+ (4 6 cos 6 - 2 sin @) sine 





i tS) 


a Oe a1 
6 = cos yn 
The special spanwise flap constraint function is given by; 


P(6), = M(8, 6,)-M(e, @,) (79). 


: ~] - : 

Where @; = cos cos} 1 o» ‘4 is the inboard station where the 
control surface begins, and "9 is the outboard station where the control 
surface ends. 


Nis 85 = 


* 


0”) sin 1/2 (6 - 6)’ 


a 
m= = (cos @ - cos @) LOG > 
sin 1/2 (6 + 6) 


M(E, 








*% 
cos 1/2 (6 + 6) 


e 
x +26 siné 
cos 1/2 (6 - 6) 








R 
+ (cos @ + cos6) LOG 


(80) 


The discrete aerodynamic influence equation (33) is transformed by 
substituting equations (40) and (65) into equation (33). — 


fon 


LA I, [A I | "3, | ll. vs| B, ele, 
= + 
[A 3, [A hp 0 P| {a P, Ivs | P ele, 
| (81) 
The constrained aerodynamic influence equation is then given by; 


CA Ia, LA Ia lls, | VS] 5 ele 
+ (82) 


LA Jp, LA Jp al. | 95} P; rel p, 


Where the arrays {a}p. and {a]p, are solved for such that the sum of the 
squares of the elements of the arrays |e Ip. and fe ID, . are a minimun. 
| J 
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Surface Velocities and Pressures 


_ The velocity tangent to the surface of the jth body in the longitudinal 









































direction is given by; 
\ N\ \ ‘ 
j ‘y} = L. T pee oe Wes +i | T K 
(tari 2 e? i “SB i a Yap, : a my. ie 
WN j he IN a B. 
Np ‘ \ | K 
‘ > i. My A 1 ary hy 1 ™, A, Ve 
2 ; aa: B Lars = 
B p B. B_P r 
Kl a ik j B Pk a JX Ve 
Pp 
K 
—i: Np = ae \ 
1 M, aes Hae 1 {7 z 
+ =- a + = 
Ka] : f oan " 'B “, ~ p B ac Ss, 
22 a8 ik x jx a Bk Ve 
Pk 
N \ \ 
T. Vv T V. T V 
X 
: | Ms, | Ve. bee F si : ” — 
jx 2 : ix © BS Vie 
AV 
M ee = = 
+ (2) B Ty + on + T : 
mn Y 
B. 
j 
(83) 


Where 








Vex 


mn Ves 


Man i Mim + 1)n 


{ =x] *b | Ie +2 | 
N: . {—_ Mn f{27 a ee os Mon *M im + 1) 
\ I *Man - ayn * Mon * aM ») | aaa | 


e bons (az in eA | 2B 
‘ Veo / ( 4 {Mon Mon +i)n 4 Mom + 2)n ) ae | 
n 


| * IMan + ayn 
he : (84) 
And 
s K K ) | 
mi _omn _@m- ijn (85) 


Veo Veo | Veo 
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Veo 


Veo 


The onset flow ratios Vy/Veo, Vy/Veo, and V7/Veo are given by the following 
expressions. 


es. pn ae 5 aa a = Tog? re 
Veo : o b (87) 
V, 7) Aiea earn 20k 3) 
Y= C.G. C.G. 
oo ee (88) 
V 20 Ge eae 6a ae | 

C.G. C.G. 
Veet ere eR eee (89) 


ce 


The components of the unit vectors tangent to the actual body subpanels in 
the longitudinal direction are given by; 


T = Pale Wl BT. + T + T. (90) 
My My > a ae 
_ Ds RE a 
T ml BT + T + T (91) 
My My My My M, 
= Das = a 
Te. 8 Dele ey Boe ee, - eT (92) 
M, M, My M, M, 


ae | 
—j 
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~ 
BD 
to 
yy! 
~x< 6B 
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ae 
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—j 
CJ BO 





34 


Bom SAN: 2 Kom + 1)N kN 
~ Veo Ne (86) 


-) 
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' The velocity tangent to the surface of the jth body in the lateral 
— , civection is given by; | 
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= (93) 


Ty Be 7 
sr a res T 1 ~ a T. . 
as (qt) = -4 | mesa) (% wt + Bins ani?) , (7 7 


a Bag” Wht 


== - _ as 
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: | 
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ih 
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If 
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And 
N 


mt K K 2 K 
m(n - An/2) . m(n- 2) _ _mn _ > m(n - An/2) (95) 
Vo Ven Van Vo 
Y = l Bs 
r N 
T (n+ An/2) K K : K 
mm - tn | m(n+1) | S m(n + An/2) (96) 
Ve Veo Vo Ven 
rel 
Py. 
Whe re 
N : N 
= 1Ki (n - An/2) — | Kae /2) 
m(n - &n m(n+ A 
2 V5 and > Ve 
r2z=i1 r=] p 
r 


are the contributions from the trailing legs of the panel vortices along the 
juncture line of the jth body and rth panel to be attached at the (n - An/2) 


and (n+ An/2) lateral stations, respectively. Ng is the total number of 
panels attached at any one point along the line. 


The velocity tangent to the surface of the jth panel in the longitudinal 
direction is given by; 








B 
‘i Se aN 2 ah es K 
(7 - fica + tan‘ A), (dZ/dX * az sax? |? | Pp p | oa | 
. P. \]] i el 4 B 
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B © ae Ven ve 
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P, 
J 
Yy 
+ —_——. 
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Where 


(Z,/C) VOI Vex A 


X ‘4 
genes —S< t ———— _M (98) 
Vn (X/C) (1 - X/C) 1+ tal) 2 
6 


and where (4Vy/Vo)mn is computed using equation (84) and X/C is the local 
percent chord. | 


The velocity tangent to the surface of the jth panel in the lateral 
direction is given by; 
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at . es ] 
ee WG “7 Y 
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where 


r UC" 
Ti (n-A2 T (nm an 
m(1 m(™ 52) 
7 r 
eu 7 1] Veo Veo OV 
7 a i ge oe a 
m4 = e. + e e. + e = mn. 
Tm m(n-1) mN m(n+] ) 
(19) 
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Veo 
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and 


represent the contributions from the trailing legs of the panel vortices along 
the juncture line of the jth panel and the rth panel to be attached at the 


lateral stations, 





(n - =) land (no) , respectivety. 


The surface pressure coefficients at each of the control points on the 
bodies and panels are then computed using the following expression. 


1/ (r-1) 


oe V V 
} See Vl x2 oy (ge, ee ee 
Ce ap [ple SMe] 1 (st) ( = | (101) 
mo Mo a 


where Y is the ratio of specific heats. 

_ Note, in equations (97) and (100) the top sign is used to compute the 
velocity on the upper surface and the bottom sign the lower surface in those 
terits which have a plus and minus in front. 
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Section and Total Loads and Moments Vg 


The section loads are computed in coefficient form on each of the bodies 
by interpolating for the surface pressure coefficients at the centroid of the 
subareas, defined by the corner points of the divisions of the body subpanels, 
and then using the following equations to sum the product of the pressure 
coefficients and directed subareas. 





Cy W ] ; 
pixie ee = ee, , 
‘i : > C,  AAY (102) 
AVG B 1K 1K 
K A. Aaf{— j 
L Ch 
K 
Ch | 
Y = od 
i; : X yy pL ee. (105) 
AVE B : 1K 1K 
K Ro st == 1 
y Ca 
K 
<> .og 
C_W aa 
: ee eee > ch a (104) 
W X P. ae 
AVG B ; 1K 1K 
K A, &{(—_ i 
Z Cc, 
K 


where i is indexed over all subareas in the longitudinals segment A(X,/Ca dy: 


th | 


(105) 
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5S phy, 


(106) 
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Also; 


and 


‘The total foade on a body § are then stained is. summing in fhe yond tad nal 
ction. 
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(107) 
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WAVC | C (108) 
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Davi = a (109) 


ee - ; (110) 


JK K K 
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a Ay Y B 
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: ee 


axe poe. af 38 
A » W C (112) 
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4) 


The total moments on a body are summed about the center of pravity. 


io, 
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1 

Sie Y..-Y a a, me eo C 
My Aye 2 : | (Yin 7 Yo.) Z. (Zix ~ 2.6.) ve | be 
; (113) 
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1 

(114) 

a X.. - X - (Y..-Y C 
My Ae 2 2 (Xx 7 Xeg.) sap ht co) hy hie 
: (115) 


The body center of pressure position vector components divided by C are 
given by; 


x (¥) ae 
= {+ an? Senet See (116) 
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(118) 
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The panel section loads, moments, and centers of pressure relative to the 
leacing edge are obtained by numerically evaluating the following integrals. 
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az = local panel width (123) 
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- 2 2 . panel width at leading 
ey. ee Oe edge. en 


and @$ = cos? [1 2 (x/C) | where (X/C) is the local percent chord. 


Also; 











tana = a ae (125) 
1 - (; ee E+ «tans 


where the top sign is used with the upper surface and the bottom sign is used 
with the lower surface. The control surface angle is equal to - 5, along the 
leading edge flap and equal to b¢ along a trailing control surface. 6 is 
equal to zero at other points on the chord. 


The section normal force coefficient is obtained by taking the following 
scalar product. 





Ci. C Cy G ae & 
AVG AVG e AVG z 


_ The section moment coefficients about the leading edge are computed as 
follows. 








C rT 
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section zero percent suction drag coefficient 


cc « 
wet 
ne Cave 





the section induced drag due to lift coefficient 





CAVC 


the section leading edge thrust coefficient Cr C/Cayg, and the section induced 
drag due to thickness coefficient 


are derived in appendix C. 


The panel total force coefficients are obtained by numerically integrating 
the section force coefficients in the spanwise direction. 


T 
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sin 6 dé (134) 
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6 = cos 4 (136) 
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The panel total moment coefficients about the center of gravity are 
numerically integrated as follows. | 


AR b a | : 
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The panel total zero suction drag, near field induced drag due to lift, 
leading edge thrust, and near field induced drag due to thickness are given by: 
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C = ——— sine de (143) 
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The panel center of pressure position vector components divided by C are 
computed as follows. 
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The section loads, moments, and center of pressure for the control surfaces 


are ccmputed as follows. 
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T HINGE LINE 
x. | 


The equation for the trailing edge of leading edge control surface; 
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*k *k; oy pig 
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~The equation for the leading edge of trailing edge control surface; 


The equation for a control surface hinge line; 


%y : .\ 
as . x, = *h, + re | (n - ns) 


1 / 


49 


(147) 


(148) 


(149) 


The unit vector in the direction of the hinge line is given by; 
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The section loads on a leading edge control surface are obtained by 
numerically evaluating the following integrals 
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where 











and C, is the chord of the leading edge control surface. 


sale 
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percent chord of the panel. 


The section loads on a trailing edge control surface are obtained by 


numerically evaluating the following integrals. 
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(154) 


(X/C) is the local 


(155) 
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The section normal load on the control surface is given by; 
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The section moments about the hinge line for the leading edge control 
surface are obtained by numerically evaluating the following integrals. 
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The section moments about the hinge line for the trailing edge control 
surface are obtained by numerically evaluating the following integrals. 
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section hinge moment is computed by the following equation. 
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section center of pressure, due to the control surface loading, 





to the leading edge is given by the following expressions. 
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The total loads and hinge moment on the control surface are given by; 
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The total loads and moments for the c 


omplete configuration are then given 
the following equations. e e By: 
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The induced drag for the total configuration Cp. is computed in the 
Trefftz plane with equations derived in Appendix F. The center of pressure 
for the total configuration is given by the following equations. 
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Fe. iS a symmetry indicator which is equal to 2.0 in equations (174), 


(176), and (178) and equal to 0.0 in equations (175), (177), and (179) when 
the panel has an image. 
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COMPUTER PROGRAM RESULTS 
Airfoil Section Velocities 


In order to establish the degree of accuracy that can be obtained by the 
use of a source-vortex lattice procedure with second order corrections to 

account for the interference between lift and thickness and to account for the 
fact tnat the boundary conditions and the pertubation velocities are satisfied 
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a 


_and computed, respectively, on the chordal plane, the results from the program 











have been compared with two dimensional exact solutions for a Karman-Trefftz 
airfoil in figures (4) and (5). Also, comparisons are made with data for a 


‘orty-five degree swept wing with an aspect ratio of five and a taper ratio of 








one. ‘The airfoil section on this wing is a twelve percent thick R.A.E, 101, 








_ These comparisons are in figure (6), (7), and (8) for zero angle of attack 














— and in figures (9), (10), and (11) at 4.2 degrees angle of attack, 
- Of the above cases discrete solutions were obtained with twenty subpanels in 


’ the upper and lower surface 


For all 











he chordwise direction and ten subpanels in 








the spanwise direction. Both 
pressures were plotted for the zero angle of 
degree of accuracy involved with the test data. 








attack data to indicate the 
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_= Section Induced Drag Due to Thickness 






















































“= The section induced drag Ca, .C/C, _, or potential form drag due to thick- 








Oe Be See av 
Ness, 15 computed in the present program 


by means of a source lattice with 

















quite well with the exact solution by R. T. Jones for 


equaticn 124) of Appendix F. The source lattice as shown in figure (12) agrees 
a en ceca a sixty degree swept, 
ten percent thick biconvex section, taper ratio one, and aspect ratio six wing 


given in reference (51). It should be noted that Woodward's equations, 























deriyed in Appendix D, for the constant and linearly varying distributed source 
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___density panels can be superimposed to obtain the same solution, for the bi- 
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-Kuchemann and Weber in reference (2). 


by Kuchemann and Weber. 














: convex ‘section, as obtained by Jones. Even though Woodward, in reference (40); 
only computes induced velocities at the 


computes induced © centroid of trapezoidal panels, it 
was shown that the correct solution for the longitudinal perturbation velo- 
City is also obtained at any spanwise location, including the edges. 





— Tf two taper ratio one semi-infinite swept panels are joined at their 
Side edges to form a wing center section or kink the Woodward distributed 
sourcé @quations will give the same center section solution as obtained by 

nar In addition, they will also give the 
correct spanwise variation of the kink effect due to thickness, which must be 
obtained by semi- empirically determined interpolation curves in the solution 
hen The Woodward equations also treat the spanwise vari- 
ation of thickness and other planform induced effects such as taper ratio, 
tip, and crank effects equally well. These same effects are also correctly 
treated by the distributed constant density trapezoidal panel equations 
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derived by Hess and Smith in reference (22), In fact Woodward's constant 
density source panel influence equations are identical to those derived by 


at Sl 
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sit rey ii 





-- As pointed out by Kuchmann and Weber the perturbation velocity due to 
thickness for a taper ratio one finite aspect ratio swept wing can be divided 
into two parts; 1) the two dimensional infinite sheared solution, and 2} the 

__ kink and tip effects. The two dimensional solution does not produce any 
section potential form drag provided the airfoil sections are closed. How- 


ea he 


ever, the kink and tip effects do produce. a section drag and thrust, respec- 
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| |” -s.samescfjvely, which when integrated across the span will give zero drag for the 

= “eamlete wing. Other planform effects such as taper ratio and cranks will 7 
) {=a ase Produce a section drag due to thickness which also integrates to zeTo for 

| [SAF of the above methods’ will give the correct spanwise distribution | 














f section potential form drag due to thickness. However, as also pointed 
ct by Kuchmann and Weber the source lattice does not give the correct edge 
ect right at the tip, kink, or crank for a finite number of source lines 


== ¥ “the chordwise direction. [owever, since this effect is only a function 
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2s -_ of the chordwise component of the thickness distribution gradient and the 
f -. yelue of the inverse gudermannian function, with its argument being the sweep 
oe ==—6; the source line, at the point where the perturbation velocity is being 
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camputed, this effect is easily added. The region on either side of the tip, 
ak, er crank which is not properly handled by the source lattice 1s a | 
function of the number of source lines in the chordwise direction and for 
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=<"“-— Seactical solutions, which require about twenty source lines per chord to 
fapresent the distribution of thickness, this region is of no significance. 
= __ Therefore, due to the superior numerical efficiency associated with the 
== sores lattice influence equations the source lattice appears to be the 
———“hest aerodynamic finite element for predicting the perturbation velocity 
{potential form drag due to thickness. , 
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Section Induced Drag Due to Lift 


"There are two basic approaches that have been tried in the past to 
solve the problem of predicting the spanwise distribution of induced drag 
or section potential form drag due to lift, 1) to accurately solve for the 


thin wing net pressure distribution, including the strength of the leading 
edge singularity, utilizing precise integration techniques to solve the 
aerocynamic influence integral equation, and 2) to utilize a vortex lattice 


«| procedure in conjunction with the Kutta-Joukowsky theorem. Both of these 








_ approaches have failed to predict a spanwise distribution of induced drag due 
to lift which when integrated is equal to the induced drag computed in the 
iT eeta te plane. The reason for this is that in these attempts the assumption 





that the vorticity is constant in the spanwise direction along constant percent 
chord lines, even for only an infinitensimal distance, leads to a nonanalytic 
influence function for which no finite value exists for the induced velocity 

at span stations where the gradients of the constant percent chord lines are 
discontinuous. 


Therefore, at span stations where the constant percent chord lines 
are kinked or cranked the constant vorticity distributed panel procedures, 
such as Woodward's, give a logarithmic singularity in the downwash. The 
lifting surface theories, such as Multhopp's or Wagner's, also produce a 
logarithmic singularity in the downwash which cannot be handled. Wagner's 
theory, as given in reference (50), has been investigated in great detail 
and is commented on in appendix E. 


In the case of the skewed vortex-lattice the downwash at the control 
point is not singular, however, the Cauchy principal value does not exist 
for the downwash on the vortex line at a span station where the vortex lines 
have discontinuous sweeps. Therefore, the Kutta-Joukowsky theorem will give 
an infinite section drag at these stations. 


All of these problems can be eliminated by using an unswept horseshoe 
vortex lattice. It is proven in appendix F that if the bound vortex lines 
are all parallel and the horseshoe lattice is evenly spaced in the lateral 
direction, the integral of the spanwise distribution of induced drag and the 
induced drag computed in the Trefftz plane are identical for all planform 
shapes. This is also true for multiple lifting surfaces, provided they are 
all parallel, and for lifting surfaces with jet flaps. 


It is not proposed that only unswept horseshoe vortex-lattice procedures 
be used to compute the net pressures or loads on wings of arbitrary shape. 
However, this appears to be the only numerical integration procedure known 
at this time which will always give the same value for the induced drag in 
the near and far fields. During studies of the error involved in using a 
skewed lattice it was determined that the error in the section induced 
drag was limited to a vary small region on either side of the discontinuity 
in the sweep of the vortex lines. Also, since most of the wing is represented 
better by skewed vortex lines (because the lines of constant pressure do 
coincide with constant percent chord lines over most of the wing) a good 
overall answer can probably be obtained at less expense, (since fewer skewed 
vortices are needed in general to represent a wing than unswept vortices) 
if a skewed vortex lattice is used to compute the net pressures and the 
unswept vortex lattice is used to compute the drag once the net pressures 
are known. 
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A comparison of the section induced drag divided by the section lift 
four different procedures is shown in figure (13). 


sphere Surface Velocity | 


The velocity over the surface of a sphere is compared to the exact 
solution in figure (14). 


X-15 Wing-Fuselage-Horizontal Tail-Vertical Tail 


Surface velocities and pressure coefficients, section force and 
moment coefficients, and total configuration force and moment coefficients 
are computed for the X-15 wing-fuselage-horizontal tail-vertical tail 
shown in figure (15). The program input for this configuration plus the 
ventral is given in the sample input section. The program output for this 
configuration is given in the sample output section. In the program output 
the fuselage, wing, horizontal tail and vertical tail are designated 
as components 1, 2, 3, and 4, respectively. 


Some results for this configuration are shown in figures (16), (17), 
(18), and (19). The data for these comparisons were obtained from 
References (59), (60), and (61). The force data is at .6 Mach nunber 
and the pressure coefficient data is at .2 Mach number. AlI1 of the 
theoretical results are for zero Mach number. 


The total configuration CL, was determined experimentally to be 
.061. The program predicts .0617. 


70 


é 


2 eS ee 





eS ee 


} 
GOHLAW INASAYd OD 


sal 


19VX3 — 
! 


‘azayuds e punoie oT}2e1 AYTOOTAA -" pl oansTy 











71 











ie 


X-15 Configuration 


Figure 15, 





yoey.y FO a[suy searseq eat 
ye aSelesny ST-X uO UOTINGILYSTG peoT TeULON - "OT s1NsTe] 


3/x 





1 
GOLA. ENaSadd | —_ 
! VING ! cy | 
ee Rann: ee ReneS Creme snoeaeey! Seamer See 








r ! ah 
) 


Prd 





SUTA| §T-X UO UOTINQTIZSIG peoy uedg Tuy -*{T sinsTy 


q/K2 


GOHIS). LNESAadd 











=e 


7 


g 
= 
og 
UF} 
az 
Rot 
< 


1°73" =e 
Meow 1 Be ' | 
7 ee 


[le], TEZUOZTIOH ST-X uO UOTINGTAISTq peoy ueds TU 
G/X2 


— 
a 








Bi) abil " i i 











- "ST 


a 








| 
i 
ii 


oind tj 


ii 
EULA | 
Hi 


ni 


hy 


| 
- 
| : 
fone es Pee Ree needs Oe i : 


i 


; 9°T 


— — 


an 





iit! 


75 


OD 


>’ 


yoeqqy Jo eTsuy seeidaq eATy ye uedstuses 
qusoleg 3°ZT #e BUTM ST-X UO UOTINqTI4STq siNssetg sSTMpIOUD -'6T Inst y 





C4501) GOHIS. LNSSTd | : | 
(udddN) GOHISIN LNASHTdd | 
(YSMOT) VIVG : | | 
(Wdddn) Vivd | 


5 


76 








i 


i) 


iit 


ia 





Ali de | 


ial 

















alt 
il 


: 














| 


l 


—oe ee - 





i 
He 











COMPUTER PROGRAM USAGE 


Program Setup 


There are two types of data decks which can be used in this program 
1) the NASA Langley input format described in reference (52) and 2) the 
NR input format described in the next two sections of this report. Lach of 
these input formats requires a different program setup. 


Program setup for NASA format. - In this form of data setup, routine 
"OUTIN' preprocesses data from Langley's format and creates a card image file 
for use by the rest of the subroutines. Changes or additions to the data 
associated with bodies and panels in the Langley form of input created by 
"MOULIN are made by body and panel "INFO" decks. There should be a body and 
pane] "INFO" deck for cach body and panel in the Langley input array, respec- 
tively. These body and panel "INFO" decks utilize the NR data array format. 


Additiona! configuration bodies and panels not described in the Langley 
data array are added by means of additional body and panel input decks utiliz- 
ing the NR data array format. 


Most of the input from the Langley data array can he directly converted 
to the NR data array format. ‘The only exception is a nacelle which is con- 
verted from a body of revolution to a ring wing. This occurs if the pod 
data has a nonzero value at the nose of the pod. If the radius at the nose 
of the pod is zero, the pod is considered a solid body of revolution. | 


The input data in the NR data array format ahead of the body and panel 
descriptions pertain to the total configuration and must always he input as 
an "TNFO"' deck. Control cards for the C.D.C. 6000 series using Scope 5.1 
are listed in figure 21. Figure 29 contains the entire deck setup. 
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ADDITIONAL 
PANELS 
IF ANY 
ADDITIONAL V 
PANEL : [F ANY 
DATA 


FORTRAN 
SOURCE 






"(DATA 
PREPROCESSOR) 
CONTROL CARD ; 


Figure 20.- Deck setup for use with NASA Langley data array format. 
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IF ANY 
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SOURCE 
DECKS 






os ead 


CONTROL 
CARDS 





Figure 22,- Deck setup for use with NR data array format. 
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Program setup for NR format. - In this form of data sctup all input 
decks utilize the NR data array format. The total configuration or universal 
"INFO" deck is followed by an input deck for each body and then a deck for 
each panel. Control cards for the C.I.C, 6000 series using Scope 3.1 are © 
listed for this form of data input in figure 21. Figure 37? contains the 
entire deck setup. 


Card No. 
Is SEQUENCL CARD 
2; CHARGE CARD 
Le PFSTAD, CM31000, T7777, P6 
4, RUN, S. 
Dis SEL. 0% 
6. LGO, 


Figure 21.- Control card deck sequence 


Input Format 


The \R data array format is given in this section. Data locations of 
all input data with limits on size of inputs are designated. A more detailed 
description of the input is given in the next section. 


80 





es é Agoty 6. DePbl Wee 








ue S Ou (*T-) sumds st Apog ATuo 51 Sure f 
Z-X 3noge Wis Sao suasuai\s Soxon poe pad AT Aoobotpur us | 
=) 23 | 


ade} wOly peor pue aeEeTNOTe Sua ten Wh gto 
Uo OAeS pue xXTIyeEW Meu aynduod (°9) :xtaqeu [y] 10} 1ORBITpU] |. 


(pact) i 
ae Apog) "ey 4 
“OX qusuoduio) } 











arte tle RE ET Rl LY eee OT . 














) _ Baed Apog i 
(ier MEA TeuOTsuaurTpUOU) (a AC} /4t = x 
(ayer ipird Teuo TsuauTpUoU) | (2/%, \c)/b = 








(oI [TAL Teuotsusutpuou) (a0/"Az)/d = ad 





MAE seit 30 oTBue) 


(yeaa jo oTsue) ® 


a enna antaaeaaeanaamnaaenae ame meamanasnamcettinan ec netadamtanie sein iiteie cd memati cnet amie 





(sq Tun) | . eer Pee a isp 
(satun) } = : "9Oy Pe as eb ae Ie 
(s3Tun) oo ——— (u Busy TeLIABI gouatayat) qh) . ae — ac 


te atte rte a ceenerteenimte Fhe Be teh Inere dette tart vite ef | cetieriihitnt tee tether bm meet mains Aapntiainitihnt A RARLA TRARY MEAN — ome 


(Y{SUueT TeUTPNZTSUOT aouatesor) 2|: 





posn joNx 


gs saith ae Neh ial gant en le pe 


_ (eee fenaen ay: 


< UPL, | me ' | | | ue | 
7 stows 0 rox b - Coe | 
ect aoapemte aos Bisa a sh Ne a Dc a eae Go eee eee ap ais ages Mi he ptih Aiea a at mee mientras eer te eee woe , aly a ie 1 ie ae 


| Satpog 30 toy Fe | a 7 
——S Ghai Ee pe Sse Pray os caemcemeted cates pala aimcatbemeestemat ein a yi gad et .. a Sous a aoa Sah 
see TESIOATUN | aa 


a tt ee 8 ee ee ene SaeIeeee 01. MINT :LSORONL OH 181 rem ve 
ear TNS LEEETIEPTONSPREOTEE PTO QTRTTNe: S| “WOTES TOES GE ELISEGLIE PER RRR NEE 8 REET TST SNORE SED cruceumeecimens easecemanseetsenmemm 4 Sons ot see Stet eee See a2 | a — apap! 


HONNd ASH LON OO NOlLd!y9S30 NOILVOISTINGGI | 4MIGANIN 
ee OR GOR TU oe ~~ 59Vd 90 BiVG —- " YSWWYY9CHd oe ee a 4330 


ViVG WWI0S0 LISIG OF OIXis. NVYLYHOS: | : 





Se te ee ee ee ed on Co ee ee ee 





saab 














ea aa 





8] 





n 





S ) 








(Op xeu) | 
= Bat eee es Gee ee 


| (662 xeu) UOTI2IIP TEUTPNITSUOT ut s2UTOd TorzUOD FO “ON! | : és 
DS Cee os ee ae | al. ae | alte Rawle cai Sk ee ee 


re rT Reeere Sr ikcieetinentiimabiaimadiaitedentemndl eee ay ai — kt pe - ~ + 
a -~-—— —-<— ' 


UOTZIOITP Tes19zeT UT squTOd TOIqUOD FO *ON , | ‘9 


Se I ae nell 








. (0z xew) SJUDWZOS ZUTEIZSUOD [R19eT FO “ON 


F i -- -* 2 om 


oe _——s- EY Te ee Se 


' (9S xeul) sjusuiZas JUTeIZSUOD TeUTpNITsUoT FO hel 


(oz xew) SUOTIOUNJ YUTBIYSUOD TRIAIT JO “ON 


: Siow, ae ~~ : ese atiabr gla GRR Sek ood 
| Ss | 
a i a CS i EIT I SSR RON oR ne ar : : } Ng es = ry Be Tengen: |b gawere spn Cage —— pe angen gacnnncnyaaraeie 
l (oS xeu) SUOTJOUNJ JUTEIZSUOD TeUTPNATSUOCT FO “ON | ; fos 


(x03aqur ppo) suOTSTATp Touedgns TexeyeT Jo “on |” fed 








ae ES 4 A S$ “te ee ee 


Seer ee ppo) SUOTSTATP Teuedqns [euTpNyTSUoT jo. “ON : 





er er 





no a ry 


I ee ce, ee ee 


ce i I a —_—— a — ome. & 


BIBI 


UOTIDIELIp TeutpnyT3u0y ut sTouedqns jo "ont 


te A Sere _ - 








i 
{ 
' 
I 
\ 
t 
t 
oJ [=] 


t eo i te re = — + a ke ey . ~w re ee me oe oe eee ee + we — ' - — , : Fs ae 
|  (s9Tun) uTzT10 Apog og “a ep 
oe eae eign at ene en eer een RN = A Ny Seam a ee ne a a ee = 


Ginn) ~_ ° | 


Pe i a lS ANS A I A My + Bh nt a mmm nem ma 


| UOT}IALIP [1IYeT UT BoUNNTJUT FO vary |. 


me ee re eee ME Re aiettter ital itnt SEAL FA RR AR RR AE a 


| UOTIIOAITP TEUTPNITSUOCT UT sdUeNTJUT FO early 
ge ie en alle Seyi, i a ce ne oT RR of nm AE MUPRERITIP cum nee come , 





BG 
Cut LM? 


err ntrerentetititr eh ja + 








— es ee a _——_ ee oe Laan anes wwe ree tee ee 


4 mle ee eae yt ne ae 


: LOTS su 11UtTjap Apog se oues (*T-) $ /S)_10 @_UgAI3 
TP PRES APP ASAP SRS (hy Sah bate ‘fon 8 rie i Be Sess s 








L _typ uaas ("1T-) {7 uora (°0) :40,e2TpuT Jouedgns ~eurpnyTsuo]_ Apog 
“stxe YX Apog 02 FT Jnd uoT}Des Apo | | 
_{*]T) aury_1a gum _o7 T and uot Das. Apoq (*O) 1O1eQ;PUL zaque> Apog 4 c. Se eotee. is 8 
BW AN, Cp Naz 2B 05 +g) OPLOTPUT 10jDeJ UOTIeOTTdTa [NU Apog - 
<T 


Z 
sea ees Secor ee Ee eg Oe ea ee AG, SRT eee” ze 


HONNd ASX LON OO = NOIdId9DS30 


NOILVOIstLNSO! YIGWNN 


Smee A COR Or ———"" 40" 73200 ra 6 Os 16 7 ' ot ae YSAWVHSONd = ‘ON WID30 


ViVG WwWiosG LISIG. O1 G3Xld) NVYLYOS 





Lt) 









82 








pee 5 uUBY 1918013 ST OP & 
Bo XPu) UOTIEOOT UT sequmu JT s,g 10 ‘Ss Has Cinta] £849 Jo ‘ON f 





4 
tea 





uotieiys by _puz 10} Sd JLUTPIOOD [RIAILT 


te 





























(Bop) Zi 
(Sep) lg 
ae tr uty. 190 eet et Qt UOT}EDO] 
AS xXeu ut tequnu Fr $,$ tC ‘s,Ta/(8av-da)] *S,8 Fo °oN 
uoTqEIs ox IST OF SOICUTPIO0D Tes97e] 
sis te a Ne ace e 
(pe32 Tuo aq t uED 4 68-It 30T ut co 
(s31un) Se. awes SY. ST tore 71). 
| (s}Tun) a ee 
Gr _suorgeis 88x Jo vox 
: _paurgap age s,[4./(p-Sy)] 40 5,9 a20yn suoTIwas X 40 3STT_ 
ail 5 
J ai ea aa eg 
bi (rca se a ee ‘ey | 
see ee reac 
y (tt xew) : SUOTZLIS X FO ON | 
: nee ree eee “STHEY STeTERIOGS Apo | 
HONNAd ASH JON OO NOld!dOS30 : 
ONEOLD 30vd LW, a ee 7 ; 


WL: TWoad: Lisia’ OW GK aide 


























ib ' 


Lu a ie 








(s3tun) —— CWg- : 


(s3 tun) "ab 


(s31un) 7 oA : 


S10308J uoTIwITTATI TN 


| - ee 





tA te tm 
rr, enon 
= : ee 
© oe 
= . ws 
ear . . rake 
ier ReMiTe vines 
aoekeey eras 
pee BN orn a 
+ ae ue 
iootects ca oktatep grtye 2 
ee Sebo 
ww ee — higiee 
ga gear ag ae gh, Le: ewig 


(s3Tun) TE 


(stun) 


(6p xeut) : SUOTIEIS X 30 “oN | 





te 


(229 pue fuotyeIs Te19IwT 





puodsas ust{z ‘UOTIVYS TeldozJeT ASATF We : ee 


UdTIDaITP TeuTpNywuoy, suote ynduy) eT v/C 7 - 47) ] 10 
T 

: aie: -_*2)] 29 

T : 


xXxeul ’ . 10 
(008 ) SUOT1EIS q. 7B SONTPA ' 2/ Oz im pal geen 


HONNd ASH LON OG = NOl1d!I¥9S30 









‘ON BOF 40 I9Vd 3170 


ViVG IWWI9S0 LIDSIO § O} 





= — WIAA DON 


Q3axt3d3o NVYLYOS 








NOILVOISILNSG I YIGWNN ) 


‘ON 930 





wm é~ 
ne 


Os 


PB) 


84 


OS £4 ADY CID FliWdAOA 





A 
-—_—T oR ee eee Te 


(satun) (6pIZ - OS8T SOT Se ous FT “4 1u0) Fay : 





ae eam an 





(s}Tun) Th. 


OO LT oT SS LS 








(6p xew) SUOTIEIS X JO ‘ON : 
POeLTFop oe ras/s) 10 chs; OL3yM SUOTIEIS Y 





TNE ere none Ne ey 





S}tun) 
(s3Tun) 
(Coos xew) 
‘ Lf 
a an ee ene ie RR ies ates esa eee eee ale peer nes of li} oO 
Ze 
(sz Tun) | Isa 
sjtun cod | bea} 
(Samm) ee — | fe 
squausoeydstp * 7G nay 00-8 1 ao 
Poaaecneen enere emery eer yt eae Dae seas eee Dee NCA ROMERO a Pe 
oe te ee Bs Pee Mite t 2 8 oe es oe ei 9 4 { eo) 
= | be} 





squauiase eTdSTP A xt 





ate ee Fe eee rere era ia ea ttt he Mlle ot aes = sama Lat + atannne wale Hades ate poerrers he 


| HONNd ABM LON OG NOUdIYOSIG | NOILVOISILNSGI | 43aWnn a ee 
“tnsiniiiniaietncinaaneninanannanasmensiine ene sibs wi Sani tse Ss in eA pees eens cco ete orev os Gero vB oN TOT 3s meee Saas jo ; ; ; i : | l 
eee) ‘ON 3 9 n ei ae a een aks rf wash = =9%d aaa ERE 3 two dane a@lucreenyenste Bos oboe. wctesnoae oe 4 WEN “ of) 50 tied ae eee mie ON) MOGd:;; . ; 


Viva WWI03G LIDIG OF aaXxls iis 























: ba | hy 
ud " i ips" 


=e: : oe ts arate “ 
a. ee ee te. 
LE NTN TITTLE LEE I AT IS “SS | SEO =} hah nmi the “ih mh ete ee seme ss ‘ “ , = : , . ag ee eee EES UNE Po 
s : ge: — 
Sea ee eee A, A (A at er: me re —— —— + . ee ete eee e 8 = ee + od Ta : Me ne eee ee eg rey i eg 
. . wa . 7 
= a 
cielo poeRet Es eke, '; 
Bre tehtye bre See . hen oe nS - " 
LENE fila ARN ha nh a rer eet op er SY OT eo Et hd La . jae Pal —_—e 








ee ee eee 





eee eo 





: 


Sl 

= 
TEES 
ies Mw) fee TY 


posetd aq 07 oper te | ee 
____-21e_S$zUTOd Jo12U09_a1eyM UOTIDeITP [e1a2eT UT sTouwtigns Fo IsTT IAs ee AO ds 








| 





| = I ate et 


g oe : 4 


TL A RH 





pooete oq 0} O1e@he - 190 8 2 3 
~—- -Squyod To. Torauo azaya uotacsatp TeurpmrBuor wr shou GIS JO TSU oa eta oar 


a ttt ter ernie | eet ee 





ee se a eae 
| Bap) eCMNS/s)_ 10 29. 





(Bop) . TeeNs/s) 10 |e | 








i] 
my, , wiih wep. iow 
Mele isi ley le) Lj Le Lele} 1S) Le 


uotae4s =P puosas 1 2 pris X8TIOA Teroqe7y | ; 








ee es we - 


(8ep) 8 RP gyg) a0 29 





| _@e) ee ee T PRU jy io lp er 
uoT}23s dg YSLTJ 3 ps x91W0A [es19Ie7 a 


a ee : dei eee et eS = meee = 4 - —~— —— er ee eee ——— ee re er ee ee 


HINNAd ASH LON OG = NOlid!I8DS 30 NOILVOISILNSGI | Y3IBWNN 


—— A BOF gO” Yd OO ”~*«*SLLG M3WWVY90Nd =~ ~~ ON DRO 


VIVG IWWID3S0 LISIG OF GAaxXIlsd NV4LYOs 








86 


Wh-2 AS Cie Deli WHO 


Bio 





rng te eT TE 


| $98 /¥x) - FS 78x) 1 41° /8x) - (H/)] cot! 
| [(O - 36) /(°% - 9o)uz}soo -6 
. (Ce - Fo) /(% - Tejuzquts +g 
ra 
9 





! | (°C /A - F(%/9x1 /[°(80/9x) - (0/0) - 
[(0o ~ Fe) /(% - Fpjyay soo - 


[(% - Fo) /(%% - U)uluts *¢ 

[(z/#9) - (z/4)]209 “p 

(z/9¢)309 -¢ 

(32 + $x) / zn /2Xwy 82 + (yxy tal + 1} /o't *2 
1° 


“ON OdUenbsas UOTIUNF JUTeIASUOD TeUTPNITSUCT E 





i 
mein ote bh i er 7 i rt A= tesa maple se a Al a te = ote 


need 





SS | 
J 


| Nc ; 
3 Ic ; 


| Sjuauses JUTeI4SUOD TeUTpnzTsucT f 


i re ae i ee re ee ie ec i he he Oe he 1 Oe ee ee 








Oe ee ee ee ee as 








HONNd ABSA LON OO NOM dSSSG NOILVOISILNSOI | YSGWNN 
petein eran AERIAL TIBET ICO ERO ff O Na 0 r Fini Gantt aoe 4o Sl hate pratt 394d MAES ota be 4 4 iva eee ee . is JWAYY9ONG, . Saale. Van ideas, PENIS ON 4330. 


ViVG IWWI0SG LISIO OL G3Xld NVYLYOS - 


87 








ae | 
al 


D 


YS-4 ABae Lt-D- Fil NAG A 








-- ee ds taliaan ceie 


1 mths 





Tt 
uoTIDUNZ JUTeIWSUOD TeUTpNyTsuoT Tetoeds 4sity Eee hE EE ne | 


eT cma ee 


[(Cu - $ujy/(Ou - Guy] 20 [(P - 30)/(% - H6)4] *6 
[(°o - F0)/Co - 46)#2]s00 -g 

[(°o - F0)/(Po - fojuzjuts *Z 

9 


ee Fre Homie | tr attire ett rev emtee nein eerie ae eer enki rarer a lr 


Ou - Fuy/(Pu - Guy} a0 [(% - Fg) /(% - Ye) 4] ‘ 


ra) on 


(0 - Fo) /(% - 4) 4} 503 z 
[(O - 40) /(0¢ - Apjajuts yp 


i nm ei 


_ zt iz +7 2n) PAN “¢ 


_ véttz AN + 2° NI /22n °z | 
OT Ti 
‘ON adUenbas uoTIDUNF quTerySsuCo Te193e7 





wo 





i) td 





(H/o 8 ouse PS ke | 
OAD 1 oozes bee . : . 





ab ee 
terres 


_(/90 "Si __ ouze Fei 


| 

| 

| 

{ 

) 

| 
BIBIEIEIGIFs 





HONNd AZH LON OO = NOlidI¥DS30 NOILVOIAILNIO] YSJSWNN 


‘ON GOP 7 0 39d TN BB SANE90Nd ON D30 


VLVG IWWID3S0 LISIG Ol QG3axXlad NVYLYOS 





88 


ae A Be fb D-Wib wed A 








SUT] PIOUD 03 OATIET OL UOTIDaTJap yuadiad (*]{) { 
B18 JO saysue IOJPLITPUT INOJUOD [Aue 
voly one (°T) 2yu8te13s (°Q) IOJBOTPUT Xa LOA BUTT TeAL 











Seneechecemne= Thathe ene atteaseeaiaa ideaanmn sen-o-apunsantamnenctetatett-ateeradinen scented 


Opts [BWUIOU SATJeESSU (“T) apts Teurcu 

SatytTsod ("Q) [aurd O} poyseqe St xXaj]1OA Butpier3 IBOgUT TSUue, 
DtiqzamMAST Ue (77) Sew ou Ce) Toue ‘gsni 51 Cy he 
Seut ue sary sumouyun pue youed 10D Tput un 
UOTIETNITeD XTTQew sseddAq (*j 


XT1jZeu Mou _oynduod (*n) xt13eu [y] eaes 03 10}e31puy. 














A A eS SRP TE ML Sy a A EE ET NN ATU = OE ea a sreguastaey maaan 


qtwo {*9), SeqeTNITeD (*[) ‘seip p[eT} esau sQeTNDTeD 03 10D TpUy 
(sjt1m tt” | 3 ueds sdualojzay 





TIOLD adUa1azay 


Bale DdUdLOJaY 


peter fae ett Nee Pome A RS 


“ON QuoUodUO) 


eleg Teueg 
tee whey Sai oh ele oe ge OPO Od 
ee _ wn COX 


a a ee ae Peel 
UOT}DUNF qUTEAI SUD Teurpn3 Tsuoy tetoeds pxrtuy | 


Cee ncn ETE iemnsamnsenamsan spearesasse asatneraa susan made londa sain siciniteniaeeeneeemenennenenemrnen ane LA mee ee ae a PORE rE e+ ceRasmeuE csseecemmaes —— et REN AA eM IR Sa Ss 





« 


Ze) f° oe ue 


pin aks md 














| ree th be atest tatessite, cats COO Od ps | 
| see siinince a] _ 
‘0 | f* z a 
UOT IDUNF JuTe1Suo0d Teurporr uo ead Puodes | fH? ¢¢ | oT 

HONNd ABM LON OG NOldidoS3SO | NOIMLWOISILNGG! 4 3QWNN | 





i ‘ON GOR 39Vd BG ate” reel MSWAVH9DOdd, ON: WIG. 


viv WWi0aG Li9IG | Ob G3Xts NVELYOs 


ol 


89 





He-2 ABs ¢1-2-Pli waco 





SOOVJINS [O1UOD aspo Sutter, jo “ON 





sooeyins ee re ident jo 





palpeqqze ST XSLIOA BUTTTE 
pipoqut vied UDTUM O2 XOVIOA Satiieii to 








loyautied Yaued Sututjep squtod Fo ‘on f0 
squtod JOLJUOD TR1ee[ JO *ON = 
“ON 


“On: 


e 


squtod TarjuOD [TeuTpnysuoT Fo 








SUOTIOUNZ YUTEIZSUOD Te1eIeT JO 


-_—_ a=. 


(ST xew) SuOTQOUN} BUTeIZSUOD Teurpn3 Tuy JO “ON | 








(Op xeu) S[auedqns Telaqe[ JO *ON © 
Oy xew S[Tauedqns [eurpnyTsucT Fo “on FE: 


UOTIDAIIP TeISIeT UT aoueMTzFUI FO wary |: 


UOTJIAITp TeUTPNITSUOT UT SdUONTFUT jo vary = 





CO 6 ee ee ee te Re 8 eee 8 ® 


(s3Tun) UIZTIQ TaUeg ee 
| (sztun) 6 a, |. 
| waAts” (“TT 


ee ee 


ov uara (*]-) uy uard ("0) AOJEITPUT xXdzIOA TeI9Ie] Jaued fe ose oe 
me 9) USATS *T ) VY UsAG (*T- pes Bae ae Peas 
¢y usAd (°G9) IOJEDSIPUT X9TIOA [eu D114 TSUO] [Louedg Prec ae oe es 


ne 
ay cen yee ae 
ae Te we 


| 
| ee eer ee 
| 


HONNd ABM LON OO = NOLdINDSIG | NOILVOISILNSC! | 43GWNN 


39Vd 31Vva = YSWWYYS0Odd = ‘ON O30 








‘ON GOFr 40 


VLVG IWWID3G LIDSIO OL G3aXld NVYLHOS 


90 


—e 


U6-24 ABeY 21-D-Ftt wes 














SARA re trea ner rma enn : oo 

| | Nal) 

| — 

| u ER ca even Ne een | ie] : 

: tw} | 

| (Of seu) suoTJe3s ueds-twas juediad jo "ON | ! 
pouTjop dle SIM pue Aoqued aloym SUOTISEUS Terzer | 





(O¢ xew) SUOTJEIS PlOYD Jued19d Fo *oOVN 


pautjap ST 1aqued a1ayM suOT Jes [BUTPNITdu0T 






























area ae 7. =f 
\) eden S ety NH 
| (satun) Z 
2 
|_(s3tun) ee a 
| ay AN) 
|_(s3tun) q, 
(s3 tun) vA 
(s3tun) | T 
(s3Tun) 4 oY 
; (s31un) DUE P AUN POUT Apoy | I 
| on UOT}DES O01 |y} UT 
{| (stun) . UdATS aes DIOYD at aS On ty 
: UOTIAdTADSap LIaqjout1ad [aueg ; 
— HONNd A3yX LON OC = NOWdIYOSRO | NOILVOISTLNO! 4YAGWNN _ 
| “AN ROR 0 29¥d ‘ON O30 


31v0q °° HaWWYHOONd 


VIVO) TWAHOJG LISI Of G3Xt4 NVYLYOS 








“t B 1 
vl Se ‘| l th 








* ) 


eee . . ead 
sir snitch memento ies. empl wae - me. ¢ rr a ee = * “ woh ee 
* pe 
ui ve bee ' Brot Sit 
. ‘ ’ 3 ' 


(08) suoT}@3s Ueds-Tulas uated JO ON] 


ee Meee lee 0: ee AE wee > athe i a ee, Mat et etn aie 
poll foep ST SSouyoOryy alum SUOTQIEIS Te13iey] ae te ae eee es ‘ ee fae eae ‘ 0 9 T t 
aan ero AE : oe ” . 5 Sanaa can a Ne Spa 0°— ‘ea gh Clas eee as fae oe ~~ 


ge nee atten 8 ea i Hiya teenage es ms = igen ae cam 





Z(2/%) 
T(3/x) 


eee 
“he 
oe 
+s 
te 
oe) 


(02) “SUOT}EIS PIOYD qusored FO ‘ON s - 
PaUuTjap ST SSAWTYI a1ayM sUOT}EIS TeUTpNyTZU0T a ar oar ales 


SU} te) Le} Ly LS} os 


KK 
wil 


TE 
— ™ 


— 








- tne oe = 

(peu) EEO Terao eG TOE Be 
ee ee eee eS Be 

(ped) [tu fe(9/x)] 9/72 to [Tu S2(a/x)jo fees 2) 


(peu) (Tu TS] 9P7 wo (Ta TEyT© be: 
OTIIUTJap SsvjAns Isques uray o : 


NOILVOISILNO! | Y3BWNN ! 


39Vd 31vG0 oo BSAWAVEDOdd ‘ON W930 


VIVO IVWIOS0G LISIG Ol Q3xXisd) NVYLYOS 


oa 
ate og 
aa . 


HONNd AZM LON OO NOlidIYDS30 


jo 











‘ON 8Or 


92 


OG-£ “ABM 11-D-11 HO A 








~_ 








or 
OO 





Tey TET TST STs] 


Oe  ) 
ee 





Oa 
et) ee 
Parad rae 

















ae tae it tal cate oa 
= eR ee » 
ease Sa he 











* 
nee he th 6 oe ee 
pi i ee one 





SUOTIEIS Toueda 





oe! oY, 
ee ere ae 





slecead : | 
TES) CIEE CETTE 





2 

















‘ON gor o 


39Vd 31V0 


viva WwWAI03ad LISI) (Oo! 

















f 












WHI 


i 






| 


It 











t 





ae 








oo] 
Oa 
1 | | 










baju 


i 
Hy 


ts 














(s3tun) 











(adBFING [OIJUOD Fo 1aUIOD “g°] 10 “g°7 psRoquy) 


(suotjouny Terseds TeEuolzIppe 1OJF Yeulloz owes asf :370N) 





OIZUO ‘ : r 
SUOTIOUNF BWTeIWZSUOD Te1ezeT TetToIeds 


(suotqed0T Avie O8Zt 


91B SUOTIIUNJ PIEpUueS :330N) 


SUOTIOUNF YUTEIYSUOD Te19RBT FO 3ST] 








39Vd 31V0 


VIVO WWIDS0 LIDSIG § Ol 





vee 
ae ee ee 


Ee Ee es ee 







Sree 
ae 


d4axtsa 


* . 


. 





= 
eee 


oo . ei “nee 
. 
if * ae 7 


MIWAVYED0dd 


NVYLYOS 


Oa-£ 











“ASW 21-3-Fit NHOS 


‘ON %930 











ide we Go UP 21 8s 


“o06b ‘O68h 28 S20S Eep ALeIS “Sad"FANS [O1QUOD [EUCTIIppe 02 anp 
SUOT3> qUTeIySUOD [PIQIsSuoT TeTIOdGS TOF WeuULIOF owes ssf] :920N 


(UOTIDVSIazUL SUT] e8uTH - aspq preoqynp) Clix 


(uOTJDasia4U] auT] a8uTy - o3pq pxreoquy) fx 


panUrtjUuC) SUOTJOUNY ZUTeIZSUQ) [eUTpMAtdUOT] TeTIEdS 


HONNd ASW LON OG = NOILdIdDS3AG 


VIVO “WWIOS0d LISIG Ol 











‘ON BOF 








OG-2 'AReM Ab-3-Pit WHOS 


+ cra eRe I ET ATT TTT TT ETE TITTIES 
as 





ret yt. 
oe ee ee ee ee oe 
ee .s 


NOILVOISILNACI 


a HS AB ot 


~ 4 SAW 90d; 








a) 
re 
a 
Le 




















* 











it 
filet 





95 














Location 
1 


2 


* 4 


10 
1 
la" 
13 
44 


15 


* 16 


17 


Input Description 


Description 


Total number of bodies used to represent the configuration 


Total number of panels used to represent the configuration. 


This number does not include those due to images when 
location 3426 has a zero in it. 


Mach number 


Total configuration reference area. 


' Total configuration longitudinal reference length. 


Total configuration lateral reference length. 

X component of center of gravity position vector. 
Y component of center of gravity position vector. 
Z component of center of gravity position vector. 
Angle of attack. 

Angle of sideslip. 

Nondimensional roll rate P/(2Vw/b). 
Nondimensional pitch rate a/ (2Ve0/E) « 
Nondimensional yaw rate Y/(2Va/b). 


The number of the component. The bodies are numbered 
first and then the panels. 


The body reference area. 


The body chord. 


* Items not needed if data format from reference (52) is useu. 
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location Description 





i 
| 
: 
B 
oo 
oo 


This is an indicator to bypass the calculation of the aero- 
dynamic influence matrices due to the vortices of this body. 
Input a zero to compute new matrices and a minus one if the 
calculation of the matrices is to be bypassed. 








#19 ‘ This is a symmetry indicator used to take advantage of both 
symmetrical geometries and loadings. If both the body geo- 
metry and the vortex strength are symmetrical about the X-Z 
plane input a zero, if only the body geometry is symmetrical 
and not the vortex strengths input a one, and if neither 
symmetry exists input a minus one. 


#20 This is an indicator used to signify whether cartesian or 
polar coordinates are used to input the body cross-sections. 
If polar coordinates are used input a zero and if cartesian 
coordinates are used input a one. 


on This is an indicator used to signify whether the lateral 
scaling factors are equal in the Y and Z directions. Input 
a zero if they are equal anda one if they are not equal. 
The ap array is not input if Bu 2B 


_— | *22 This indicator is used to signify whether the body cross- 
sections are placed perpendicular to the mean camber line 
or the X axis. Input a zero if the cross-section is to be 
placed perpendicular to the mean camber line or a one if it 
is to be placed perpendicular to the X axis. 


23 This indicator signifies the type of subpanel spacing used 
on the body in the longitudinal direction. Input a zero if 
te wee the subpanels are to be spaced at even increments OF 0; 
et ee ~ “where d=cos-1 [1-2(x/c)], input a minus one if the subpanels 
ee are to be spaced at even increments of X, and input a one if 
the’ spacing is specified at a given set of xX stations. The 


es ae cia > 


Xs stations are to be given in locations 1850 - 2149. 











2Hss 5 FES === 22 ee a 


24 : This ae ee the Syne ee subpanel spacing oe 
- on the body in the lateral direction. Input a zero if the 
subpanels are to be spaced at even increments of 8 or "IB = 
(S/Smax) a one if the subpanel side edges are given at a 
edges are at the same lateral stations as the body geometry 
, is defined. The subpanel side edges are input in locations 
ee 2200 - 2799. 








* Items not needed if data format from reference (52) is used. 
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Location 


25 


26 


27 


#28 


*29 


30 


31 


bY. 


Description 


This is the longitudinal distance, in terms of the length : 
of the influencing quadrilateral, at which the contribu- Cea? 
tion of the quadrilateral vortex to the perturbation 

velocity is considered small enough to be neglected. 

Points at a longitudinal distance, away from the centroid 

of the vortex, larger than this value are not cycled 

through the influence equations in order to save computer 

time. 


This is the lateral distance, in terms of the width of the 
influencing quadrilateral, at which the contribution of 
the quadrilateral vortex to the perturbation velocity is 
considered small enough to be neglected. Points at a 
lateral distance, away from the centroid of the vortex, 
larger than this value are not cycled through the influ- 
ence equations in order to save computer time, 


This is the X component of the origin of the body 
coordinate frame. 


This is the Y component of the origin of the body 
coordinate frame. 


This is the Z component of the origin of the body 
coordinate frame. 


This is the number of subpanels in the longitudinal direc- 
tion on the body. 


This is the number of subpanels in the lateral direction 
on the body. 


This is the number of equally spaced divisions the body 
subpanels are divided into in the longitudinal direction. 
These divisions are used to compute the subareas which 
in turn are used to integrate the surface pressures to 
obtain loads and moments. This number must be an odd 
integer. 


* Items not needed if data format from reference (52) is used. 
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si aa is 50. 


Description 


This is the number of equally spaced divisions the body 
subpanels are divided into in the lateral direction. 
These divisions are used to map the vortex grid closer 
to the actual body surface and to compute the subareas 
which in turn are used to integrate the surface pressures 
to obtain loads and moments. This number must be an 

odd integer. 


This is the number of functions used to constrain the body 
surface vorticity in the longitudinal direction. The 

list of functions used is given in locations 3210-3249. 
These same functions are used over each longitudinal 
constraint function segment. The list defining these 
Segments 15 given in locations 3140- 3189. The maximum 





This is the number of functions used to constrain the 
body surface vorticity in the lateral direction. The 

list of functions used is given in locations 3250-3269. 
These same functions are used over each lateral constraint 
function segment. The list defining these segments is 
given in locations 3190-3209. The maximum number of 
lateral constraint functions is 20. 


This is the number of segments in the longitudinal direc- 

tion over which the body longitudinal constraint functions 
are defined. The list of subpanels defining the segments 

is given in locations 3140-3189, The maximum number of 





‘This - is ‘the number of segments around the circumference 
of the body over which the lateral constraint functions 
» are defined. The list of subpanels defining the segments 
is given in locations 3190-3209. The maximum number of 


ee is 20. 


This is the number of control points in the longitudinal 


-. direction on the body. The maximum number of points 


in the longitudinal direction is 299. The list of body 
control points in the longitudinal direction is given — 
in location 28007 3099. 


* Items not needed if data format from reference (52) is used. 
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Location 


39 


*40 


*41-89 


*85 


*86-129 


m130 


*131 


Description 


This is the number of control points in the lateral direc- 
tion on the body. The maximum number of points in the lat- 
eral direction is 40. This list of control points is given 
in locations 3100 - 3139. 


This is the number of longitudinal stations where the body 
cross-sections are defined. The maximum number is 44. 


This is the list of longitudinal stations where the body 
cross-sections are defined. This list starts at the nose 
and ends at the aft end if the body is solid. For a flow 
through body this list begins at the aft end, continues 
along the inner surface to the nose of the body, and then 
along the outer surface to the aft end. 


This is the number of longitudinal stations where the body 
cross-section lateral stations are defined. The maximum 
number of these longitudinal stations is 44, 


This is the list of longitudinal stations where the body 
cross-section lateral stations are defined. The body cross- 
section lateral stations are constant between each longitu- 
dinal station in this list. If this list 1s the same as 
that in locations 41-89, it can be omitted. If the body is 
solid, this list starts at the nose and ends at the aft end. 
If the body is a flow through type, this list begins at the 
aft end, continues along the inner surface to the nose of 
the body, and then along the outer surface to the aft end. 


This is the number of lateral stations at the first longitu- 
dinal station in the list at locations 91-129. The maximm 
number of these lateral stations is 39. 


This is the list of lateral stations at the first longitu- 
dinal station in the list at locations 91-129. These lat- 
eral stations where the body is defined can be given by the 
angle 6, the lateral fraction Lg ADB)» or fraction of 
lateral circumferential length §,? 


Note: Lists of lateral stations at the other longitudinal 
stations, given in the list in locations 91-129, follow 
this data as shown in the input format. 


* Items not needed if data format from reference (S52) is used. 
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Location ena: 

















#300- 1599 In these locations the body. cross-section radii or frac- 
tional distances [(%-A%V2p,,], for each of the Xp; 
longitudinal stations and oe or [Q3,- sy lateral 
stations, are input. 

















*1600 This is the number of longitudinal stations where the 
multiplication factors %, and 2,, and the translation 
increments AY and AZ are given. The maximum number of 
stations that can be used here is 49. 


1601-1649 This is the list of longitudinal stations where the mul- 
tiplication factors Yay, and Zp, and the translation incre- 
ments AY, and AZ,are given. 


» 1650-1699 This is the list of multiplication factors Z,,at the 
ee Jeng! fading! stations given - in locations 1601-1649. 


1750-1799 This is. ‘hie list of translation pre AY, at the 
oe | peng cueinal Stations given in locations 1601-1649. 








* 1800-1849 This is the list of translation ‘increments AZpat the 


Pe stations given in locations 1601-1649, 








* 1850-2149 This is the list of fenieaiaat’ stations where the body 
i wiser edges are defined. Oo 
2150 “This i is the number of longitudinal ‘stations where the = 
lateral subpanel edges are defined. The maximum number 
of these longitudinal stations is 49. 


2151-2199 This is is the list of longitudinal stations where the body 
lateral subpanel edges are defined. This list can be 
omitted - if it is the same as that in locations 1850-2149. 


2200 This is the list of lateral ial edges at the first 
longitudinal station given in the list at locations 
2151-2199. These stations can be defined by 6's 

or. fractions of body circumference "b= (S/Spax) - 


Note ; Lists of lateral subpanel edges at the other longi- 
| ‘tudinal stations, given in the list at locations 2151- 2199, 
follow this ae as shown in the input format. 


-* — Items not needed if data format from reference (52) is used, 
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Los Angeles Division 
North American Rockwell 


Location 


2800-3099 


3100-3139 


3140-3189 


3190-3209 


3210-3249 


3250-3269 


3270-3419 


Description 


+ 


This is the list of the number of the subpanels where the 
control points are located in the longitudinal direction 
on the body. The integer in this list indicates the 
number of the subpanel, aft of the nose of a solid body 
or aft of the tail end on the imer surface and around 

to the outer surface of a flow through body, at which a 
control point is placed. If no longitudinal constraint 
functions are used, this input can be omitted. 


This is the list of control point locations in the lateral 
direction on the body. The integer in this list indicates 
the number of the subpanel in the lateral direction from 
the top of the body at which a control point is placed. 

If no lateral constraint functions are used, this input 
can be omitted. 


This is the list of longitudinal constraint segment bound- 
aries. The integers in this list indicate the subpanels 
between which the longitudinal constraint functions are 
applied. For example, if this list included the long1- 
tudinal number of every tenth subpanel, the longitudinal 
constraint functions would be applied over a range of ten 
subpanels and repeated over the segments defined by every eet 
tenth subpanel. 


This is the list of lateral constraint segment boundaries. 
This list is utilized in the same manner as the longitudinal 
constraint segment boundaries. 


This is the list of longitudinal constraint functions used 
on the body. The functions to be used are signified by 
entering here the sequence number of the function as given 
in the input format. If special functions are to be used 
the data location where the function is described is input 
here. 


This is the list of lateral constraint functions used on the 
body. The same procedure is used here to signify the 
desired functions as is used in the definition or the 
longitudinal constraint functions. 


In these locations the special longitudinal constraint 
functions are described as shown in the input format. 
These are described in tabular form. ( 


* Items not needed if data format from reference (52) is used. 
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3420 This is the component number for the panels. 
#5421 This is the panel reference area. 
3422 This is the panel reference chord. 
5423 This is the panel lateral reference length. 
5425 This is an indicator to bypass the calculation of the 


aerodynamic influence matrices due to the vortices of this 
panel. Input a zero if a new matrix is computed and a 
minus one if the calculation of the matrix is to be 
bypassed. 


*3426 This is a symmetry indicator used to take advantage of 
a configuration's symmetry to save input effort and com- 
puter time. Input a zero if both the panel geometry and the 
vortex strengths have images on the port side of the X-2 
plane, a one if just the panel geometry has an image, a 
minus one if neither have an image, and a two if the vor- 
tex strengths on the panel and its image are antisymmetric. 


3427 This indicator is used to signify whether the panel is 
attached to the positive normal or negative normal side 
of another panel. Input a zero if it is attached to the 
positive normal side and a one if it is attached to the 
negative normal side. 


3428 This indicator is used to signify whether the wake is 
force free or fixed. Input a zero if the wake is fixed 
and a one if it is force free. If any panel has a force 
free wake, the wake or a body it is attached to is also 
assumed to be force free. 


#3429 This indicator signifies whether the mean camber surface 
of the panel is described in terms or local angles or 
attack or as fractions of chord. Input a zero if local 
angles of attack are input and a one if fractions of chord 
are input. 


' +» 


* Items not needed if data format from reference (52) is used. 
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Location 


3430 


3431 


*3432 
*3433 
*3434 
3435 
3436 
3437 
3438 


3439 


3440 


Description et 


This input indicates the type of longitudinal subpanel 
spacing that is used. If the subpanels are spaced at even 
increments of ¢ input a zero, if they are at equal incre- 
ments of percent chord input a minus one, and if they are 
input as a given set Of percent chord use a one. 


This input indicates the type of lateral subpanel spacing 
that is used on a panel. Input a zero if the spacing is at 
equal increments of n, a minus one if it is at equal 
increments of 96, and a one if the spacing 1s specified at 

a given set of n Stations. 


This is the X component of the panel origin point. 
This is the Y component of the panel wien point. 
This is the Z component of the panel origin point. 
This is the area of influence of a quadrilateral vortex 


on a panel in the longitudinal direction. This input is 
in terms of fraction of quadrilateral vortex length. 


a ) 
“ 


This is the area of influence of a vortex on a panel in the 
lateral direction. This input is in terms of fraction of 
vortex width. 


This is the number of subpanels in the longitudinal direc- 
tion on the panel. The maximum number of longitudinal 
subpanels is 40. 


This is the number of subpanels in the lateral direction 
on the panel. The maximum number of lateral subpanels 
is 40. 


This is the number of longitudinal constraint functions on 
the panel. This number does not include special functions 
due to control surfaces. The maximum number is 15. 


This is the number of lateral constraint functions on the 


panel. This number does include special functions due 
to control surfaces. The maximum number is 20. 


* Items not needed if data format from reference (52) is used. 
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Location Description 


S44] This is the number of control points in the longitudinal 
direction on the panel. The maximum number is 40. 


3442 This is the number of control points in the lateral direc- 
tion on the panel, The maximum number is 40. 


"3443 This is the number of points used to define the panel 
perimeter. In the root section the perimeter is defined 
by X, Y, and Z components of points along the leading 
and trailing edges. The outboard section is defined by 
X, Y, and Z components of points along the leading edge 
and the local chord. The root section leading edge is 
input first, then the outboard section, and then the root 
section trailing edge. The maximum number of points is 50. 


5444 This input is the number of the body or panel to which this 
panel's inboard trailing vortex is attached. 


5445 This is the number of the trailing vortex leg or subpanel 
side edge to which this panel's inboard trailing vortex 
is attached. The subpanel side edges are numbered consecu- 
tively starting at the top of a body or at the inboard 
edge of a panel and going in the clockwise direction when 
an observer is looking in the negative X direction. 





3446 This is the number of the panel to which the outboard 
trailing vortex of this panel is attached. 
3447 This is the number of the trailing vortex leg or subpanel 


side edge to which this panel's outboard trailing vortex 
is attached, 





S448 ‘This is the number of leading edge control surfaces on the 


3449" == This is the number of trailing edge control surfaces on the 


_-*3450-3599 In these locations the panel perimeter is described. 


* Items not needed if data format from reference (52) is used. 
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cr 


Location 


* 3600 
*3601- 3629 
* 3630 
#3631-3659 
*3660-4099 
4100-4129 
* 4130 
*4131-4159 
* 4160 
* 4161-4189 
* 4190-4599 
4600-4639 


4640-4679 


4680-4719 


Description 


This is the number of longitudinal stations where the 
panel camber is described. The maximum number is 30. 


This is the list of chord fractions which define the 
chordwise locations where the panel camber is defined. 


This is the number of lateral stations where the panel 
camber and twist is defined. The maximum number is 30. 


This is the list of lateral surface length fraction n 
where the panel camber and twist is defined. 


This is the table of local angle of attack or deflection, 
in terms of fraction of chord, of the panel mean camber 
Surface, The table is input as shown in the input format. 


This is the table of panel twist. 


This is the number of longitudinal stations where the panel 
thickness is described. The maximum number is 30. 


This is the list of chord fractions which define the 
chordwise locations where the panel thickness is defined. 


This is the number of lateral stations where the panel 
thickness is defined. The maximum number is 30. 


This is the list of lateral surface length fraction n 
where the panel thickness is defined. 


This is the table of panel thickness in terms of fraction 
of chord. The table is input as shown in the input format. 


This is the list of longitudinal subpanel edge locations 
in terms of fraction of chord. 


This is the list of lateral subpanel side edge locations 
in terms of n. 


This is the list of panel control points in the longitudinal 
direction. 


* Items not needed if data format from reference (S2) is used. 
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Location 


4720-4759 


4760-4779 


4780-4879 


48380 - ... 


Description 


This is the list of panel control points in the lateral 
direction. 


This is the list of exponents W for the standard lateral 
constraint functions or data locations of the special 
lateral constraint functions. 


In these locations the special panel lateral constraint 
functions are described. The input is as shown in the 
input format. 


In these locations the special panel longitudinal constraint 
functions are described. The input is as shown in the 
input format. | 


* Items not needed if data format from reference (52) is used. 
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CONCLUSTONS 


The theory and program described herein are capable of predicting surface 
velocities and pressure coefficients and section, control surface, and total 
configuration loads and moments for a diverse class of airplanes. In particu- 
lar the source - vortex lattice with second order corrections was shown to 
agree well with the exact solution for a Karman-Trefftz airfoil section. 


_ the quadrilateral vortex element predicted surface ve TOCTE LES which agreed 


Also, 
well with the exact solution for a sphere. 


Because of the options built into the program simple configurations can 
be input with a minimum of effort while the capability exists to input com- 
plicated configurations in significant detail. The program is general enough 
to analyze complete configurations at angle of attack and in side slip, 
pitching motion, rolling motion, yawing motion, and with control surfaces 
deflected. The configuration can be run in each of these modes to obtain 
static and rotary stability derivatives or in a combination of the modes to 
predict the loads and moments on the configuration while in a quasi-steady 


maneuver. The program also has the capability to account for a free wake 
while operating in any of the above modes. 


4 


Within the scope of this study a significant effort was devoted toward 
understanding the limitations of and the relationships between the different 
types of aerodynamic finite elements and lifting surface theories available 
at present. It has been concluded that those used in this analysis are 
presently the most efficient and generally available. It was also shown that 
the spanwise variation of potential form drag due to both thickness and lift 
are correctly computed by the program described in this report 
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Appendix A 
NUMERICAL PROCEDURES 


Discussions of the prime numerical procedures used within the program are 
given in this appendix. There are essentially three such procedures; (1) 
straight line interpolation and extrapolation, (2) controlled deviation inter- 
polation, and (3) Householder’s simultaneous equation solution. 


For straight line interpolation and extrapolation about two given points 
(X1, Yy) and (X2, Y2); 





Y= @Y> + (1-a@) Y1 (1) 
where 
X-X] 
Qos 
x (2) 
_dY , : 
The slope aX for this case is given by; 
dY ' , 
where 
a’ = 1/(X,-X,) (4) 


In the case of the controlled deviation interpolation method (CODIM) para- 
bolae are used to curve fit a set of four given points (Xy-1, Yn-1), Qy, Yn), 
(Xyt1, Yytl), and OG+2, YytZ) to obtain interpolated Y and dY/dX values for 
Xy sx Xyt1. Only that information, relative to this method, which is necessary 
to judiciously pick input points will be discussed here. A complete derivation 
is given in reference (20). 


a* 
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. One parabola P is fit through (Xy-1, YN-1), (iy, Yn), and (XNe1, YN+]) = 
=.) The other parabola P2 is fit through (Wy, Yn), One | Yn+]), and (XN+2., YN+2) « 

This curve fitting process involves the solution of’two sets of simultaneous 
equations. If, 


Py) = AyX* + ByX + Cy (5) 
and 
P2 = A?X* + BoX + Cy (6) 
Then 
2 
Ay yer Wa} YN-1 
B) S$ = | X4y Xy 1 Yn (7) 
Cy Xap XN YN+] 
and 
z= 
A2 Xn Xy I Yn 
Bo = XN] XN+1 I YN+] (8) 
C2 X’ne2 = -Xne2— 1 YN+2 


The interpolated values of Y and ay/ax iB aicen Xy and Xy+1 are defined sie 
either Pj, P2, or a linear combination of P,; and P27. The amount of P) and P2 
used in the linear combination is determined by comparing both of the parabolae 
to the straight line. 


S= a Yye] + (1- @) Yy (9) 


” 
“io ix 
bail 
one 
ee, 2 ee 
—__ 
——_—aee 
rr eee eee ~ 
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: 


Therefore; 


dD _ da dE; _ da _q) 22 
ax Tiax 2) Pa axe eo ay oH) 
P) for X = Xn 
@E)P2 + (1- @)E2P 
: <4 ee E, + (1 ayEy for Xy X Xyyy (19) 
P» | for X = XN+] 
and 
dPy 
dx for X= Xy 
dY dN d)).. 22 
Fz] ge /D Nay /D° for My OX Aye (20) 
dee —— 
— dx for X = Xy4] Pe 
In the case of an end interval Xy-1*x $Xy; Pj is set equal to; 
P, = S + K(P2 - S) (21) 
where 
In - [M2 | 
S413 Ss (22) 
(Mp + [Me | 


(- 
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where | 


(10) 


The parabola which has the least deviation from the straight line is given 
the greatest weight. The weighting factors E) and E2 are determined as follows: 


Hl 


dH 














e weighted expression for Y in t 


<= ———— 
a Ee. ES 


e range Xy<X< Xy4] is then; 


















































The derivative dY/dX in the range X, = X <Xy47 is then; 


where 


And then 


dN 
ax 





dN dD /2 
as ~ aio 0) 
x) Nx D° for DF 





P. + (1- @) EOP, 


aE, + (l-a) E. 


dE @2 4 dB p yp, Bi 
SEP? + oF P> + a ELS E2P, + (1- a )ay P,tE2 aX. 


(11) 


(12) 


(17) 


and 


1 Xy > Xn-1 
6 oe (24) 
Xy - XN+] 


A similar procedure is followed for the other end interval Xn+1* x sXye2- 


Householder's method for solving simultaneous equations is used in the solu- 
tion of the aerodynamic influence equations. The method is applicable to both 
square and rectangular influence matrices. In the case of rectangular matrices 
it is not necessary to least square the equations first, since Householder's 
procedure least squares and triangularizes simultaneously. Also, the influence 
matrix is triangularized by means of orthogonal transformation matrices, which 
preserve the conditioning of the matrix. The combination of these two advan- 
tages, along with a reduction in the number of required computer operations, 
greatly improves the numerical accuracy and stability of the solution over that © 
of the standard Gaussian reduction method. 


A complete, but rather abstract, derivation of the method is given in 
reference (21). The method in the subroutine has been altered from the original 
to allow the operation on a single row of the matrix at a time. This reduces 
the required core allocation necessary to triangularize the matrix. 


_ A derivation of the method, developed by the writer, will be given here in 
order to describe the basic philosophy of the method. 


If [A] is the rectangular influence matrix, the upper triangle is given by; 


[R] = [WwW] (Al (25) 


where [W] is the combined orthogonal transformation matrix used by Householder 
to triangularize [A]. 
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The relationship between Householder's triangularized matrix [R] and that 
obtained by Gaussian elimination of the least squared influence matrix [A]![A], 
1s 


[G] = [T] fA]? [A] = [T]IRJ![R] = [DI [RI (26) 


where [G] is the triangular matrix obtained by Gaussian elimination of [A]T[A]. 
The matrix [T] is the Gaussian transformation matrix used to triangularize 
fA] T [A]. And, [D] is a diagonal matrix with the same diagonal as [R]. It can 
be seen from equation (26) that a nonsquare matrix must be least squared first, 
before applying the Gaussian transformation [T]. Whereas, the Householder 
transformation matrix [W] can be applied directly. The least squared matrix 
(A]![A] is usually more illconditioned than [A], and therefore, less accurate 
results are obtained. 


In the Householder method [W] is equal to the product of N+] individual 
orthogonal transformation matrices, where n equals the number of unknowns. 
There are N+1 transformations because the augmented influence matrix, made up 
of the influence matrix itself plus the boundary condition matrix, added on as 
the last colum, has N+1 columns. Each transformation results in reducing all 
elements below the diagonal to zero for one colum. The colums are reduced 
from left to right. 


The individual transformation matrices [W],, are defined by; 
(W),, = CIT] - 2 {ulm {uli (27) 


whére [I] is a unit diagonal matrix and [u], is a column matrix defined by the 
unit vector Up = (am -@pVy)/ By. The vector 4, is defined by the mth colum 
of [A] where the elements on rows less than m are replaced by zeros. The unit 


vector V,, is defined by a colum matrix [v],, with all zeros except for the mth. 
row, which is equal to one. The constants a, and uy, are defined as, 


=| faa) | (28) 


By = V20n(G_ - Vigt am) (29) 
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It can be shown that [a,] is reduced to |{am}|{vm} if {am} is premultiplied by 





({1J -2{un| lu T), Also, that the first m-1 rows of ex 
(Wn-1](%n-2] +++ [Wy] TA] (30) 
remain unchanged by the mth transformation. The result after m transformations 
is then zeros below the diagonal for the first m columns and |{aj}],|{a}2], ... , 
l{a}m-2]. | {a}m-1]>{a}m| on the diagonal. The elements above the diagonal have 
been defined by the m preceding transformations and will remain unchanged for 
the N+l-m remaining transformations. 
((1] - 2fum Ham} {@m) = [amt lt yal (31) 
Un = (qn - &_ Vm)/ Hm Or fU_) = Caml -Om{Ypl/ em (32) 
Om = |14nl| (33) 
) aed 
py V2am(@m- Van) OF bg = Y2enCa- - fat lv_l) (34) 
remain to be proved. It is helpful in the derivation of equation (31) if the 
vector identity 
JaIV + 2G yy) thy = By (35) 
is observed from the following vector diagram. 
2, * ty 
om Lan Yip 
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Then from equation (35); 


[al ¥in + 2 UnGm © Ay) = 4, (36) 

anl¥n + 2 Ud, ° am = am (37) 

Therefore; 

Lan lVin = Gt ~ 2 Unt) 8m (38) 

A A ee ; IN baie ok , - 

where 1,1, and u,u, are dyadics. The unit vector i, is in the direction of a,. 

Equation (38) can then be written in matrix notation as follows; 

((1] - 2 fulm(uj?m {alm = Halmi{vim (39) 

~<— which is equal to equation (31), In matrix or tensor notation it becomes —_ 


evident that the dimensions of /a},, {v/,,, and ju}, are not limited to three. | _ 
}4 pms [¥ jm m : 
On = [ant (40) 


and . 


tn = 2 {aml fag! (41) 





Then if equation (31) is premultiplied by | 2m}? 
(og) {a} - 2(4q)7 Un) tn) = Hag! Hand | (42) a 
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And substituting a, and Ln into equation (42). 


ay T 
a2, - se mn {4m} {Yn 


Or 


oy 


Hy = 2am - (aj {va} 
In vector notation equation (44) is seen to be equal to; 


Um = FOC a Vege) 


Also, if equation (44) is substituted back into equation (31) 


Fo es aaa Se ; 
{an} - “2a_(Om - fant t{val) {am} = Om iV} 


Therefore; 


{a} i. antYnt 
{u_| Ss 
Y er a ee oa 
Zanlam - {amt i {vm}) 
or in vector notation 

4 am - 8m Vin 

u = 

m 


V2Za,(am - V_‘4,) 


™m 
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(43) 


(44) 


(45) 


(46) 


(47) 


(48) 


( 


JT 
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Appendix B 


SKEWED SOURCE-VORTEX LATTICE INFLUENCE EQUATIONS 


£ 


The general form of the skewed source and vortex lattice influence 
equations was derived in appendices B and C of reference (26), respectively. 
These equations will be specialized to the case where the inboard and outboard 
sweep angles for a given horseshoe are equal. The perturbation velocity due to 
a skewed vortex is then given by the following equations. 








bs ae (1) 
Ve Von | 
V Wy 
vy AY, Ve AL, Ve 
<= = : (2) 
(a8) 
Z 2 2 ? 
AY, + Aa, AY.” * Ad, 
V. W 
a aad 
ae aa: ( oO oy, (AY 
” + (3) 
Z Z 2 2 
AY, + AZ, AY, + AZ, 


Where AY., and AZ, are the changes in Y and Z across the width of the 
horseshoe vortex, respectively. 








U (T/Vq )Ey 
VL Vv 
Voo - Ar (4) 
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V (T/Vo JEy, 








2 5 
Vo an (5) 
| Vo )E 
My 7 CT/Veo Ew P 
—_— == (6) 
Vn An 
Where 
=, a (7) 
mu 7 5 (1, + 15) 
Vv R 
2 
I,+1 U,+)D = 
Be ee es ee, a) (8) 
V 2 Rn ¢ Rr ¢ 2 3° 
i Ri 3 2 
W-Ryy +) Ore wIG*YD SYN" 
fie ee eee ee a2 +: 
: Ry 1 2 
Let 
(lect Te) 
ee, ee eee (10) 
pie 
2 
Les 
term 2 = (11) 
2 
Ry 
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—— ae 
ae 
——————————— 
ee 
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1 


| 


H 


| 





rrr ‘: ——= 
—— 





J (12) 





Then 


Eu, Z (term 1) (13) 


I 


{| 


Wy z {-T(term 1) + (term 2) - (term 3) | (13) 


Buy -(X - TY)(term 1) - (V+ B yy) (term 2) + (Y-B Yp )(term3) (24) 


Where 


= X+ Typ (15) 











ee 2 os 7 7 
Y+TX-8y O+T’) oT Xo eet 
_Y+TX-6y GQ+T) 2. Y+#TX-h Ww OT) 


— — —2| 1/2 40 
lK&-Ty,)?+ W-2 yee 


ui 


t:, 2: tees 


bc pce te _ X-TYy (18) 
_ — —21 1/2 4 
| - 1 yo se 8 yy)2+7| / 
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and 


o<| 


—< | 


Al 


2.3 ee 


= (Y + By). 22 


24 Wye 7 ee 


- - By, +2 


roe! 
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— 2 = 7 Me ac? 
= &-Ty,)*@-By, ) +2 
— 2 - 2 =2 
SOP iys) 2 ety 2 
X - X 
q Vv 
AY. (Y - Y. AZAL..<— 6 
P vq Fs vq" WW 
Z 2 2 2 
AY * Ad Ax + AZ, 
AY. (Y - ¥ joy ae 
- B ae: V) + B v3 
Z Z Z 
A 
Ny + AZ, AY, + AZ, 
2 Z 
1/2 AY, + AZ, 
tan A 
B 


) 


Ny 


Z 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


where (Xy, Yy, Zy) and (Xq L q) are the locations of the influencing point 
and the point being seth reepectively. 


ft 
YX 


‘ 
The elements of [Ay ] , [Ay] , and [Az] are computed for a unit strength 
of YT/V,.- , 








Similarly, the perturbation velocity due to a ‘skewed source line is given 
by; H 
Se ek lie 
Veo Ven 
ay Pe: 
Wee ye L, Veo 
—— 3 
7. (30) 
—_ Vg Ns 
W Vv oe 
ee eo ge (31) 
| 2 2 | 2 2 
AY + Ai, AY, + Az, | 
where 








(2 Ve Eu : 3 
Veo An 
. (x /V,, )Ev — 
16 an 
We (z /V., Ew 
Ve . 4m (34) 
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and 





Ey. = d (term 4) + ire see 
ict \ 1+ 1 


Ey, So ee ee Cem 4) 
Ji+? (i+? 
—? _ 
Ey = 1+T° Z (term 1) 


and where 


1 1 


"R, R, 


4 2 


term 4 = 


1 rT 


(X - TY) (term 1) 


(X - TY) (term 1) 


(35) 


(36) 


(37) 


(38) 


The elements of [ Sx | ; [ Sy | » and [ Sz | are computed for a unit strength 
of 2£/V,,. The influences of both the vortices and sources at the quarter 
chord of the subpanel are computed simultaneously due to the similarity in the 


vortex and source influence equations. 


For the case where |Z | s2yy , set Z = 


0 and if both 


&-mwy/[a+ty, 2+ F+ ey, 07] 


and 


&-?/[a-Ty, + &- 
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B yy )?| 


(39) 


(40) 


are less than (0.08716}2, and lY| >By,, , then use; 


c 





ee ee d : 
= 2 zs - = 2 
21 + F K-Ty,)2+%-Byy)? &+Ty,)° + T+ By,y) 





(41) 


in place of (term 1) 


If ¥| sBy, and | X - TY | < By/4 set (term 1) equal to zero. 
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Appendix C 


QUADRILATERAL VORTEX INFLUENCE EQUATIONS 


The Biot-Savart law can be used to calculate the influence of a finite 
vortex segment on a point in three-dimensional space. The incremental change 
in induced velocity at a point in space due to an incremental change in length 
of a finite vortex is given by the following expression. 


where; 
K 


h 


KCcOS? d?@ 
47h oY) 
Vortex strength 


Perpendicular distance from the vortex segment to the point in 
space. 


Angle between the line formed by h and a line from the field point 
to a point on the vortex segment. 


Velocity induced by the finite vortex segment perpendicular to the 
plane formed by h and the vortex segment. 


A vector expression for q can be determined from figure (D-1): 


=. 





Figure C-1.- Velocity induced by finite vortex segment. 
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| The magnitude of the velocity q induced at (Xq, Yq» 2 ) by the vortex segment 


S is given by the following equation after equation (1) has been integrated 
— from 9; to $f. 


= k 
a | = amp (605 B - cos a) (2) 
where 

Sak 

cos B= TE 
5 Re 
3.3, 

COS @ * 7 WZ 
S R, 

i 











The vector h 


( 


is determined such as it satisfies the conditions of being 





perpendicular to S and,equal to the vector sum h = R¢ - as where "a!" 
defines the length of h. 
ae Gs _ -Since, 

h = Re -as (3) 
and 

h-s = 0 (4) 
then 

h+s = R, s-as*s = QO (5) 
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lk 


therefore; 


b 


os) 
nf 





wf 
wy 


— 
After substituting ''a" into equation (3), h is defined as; 





R,-* S$ 
> > + — 
h = Re- Zs S 

5° 5 


Also, a unit vector q in the direction of q is seen to be equal to, 


Re x s 


ik, x 5| 





q 


The magnitude and direction of q are then expressed in terms of the 


(6) 


(7) 


(8) 


coordinates of the control point (Xq> Yq> Zq) and the end points of the vortex 
segment (Xj, Yi» Zi) and (Xf, Yf, Ze\e Eh dar “Js and k are defined as unit 


vectors inthe X, Y, and 2 directions respectively, then; 


— 


wh 
T 


S.° 5, = Oe X)i* pi, -Y,) 5 + 8@,-2,)k 
R, = 5,-Q = O-X)i * ris wa* p (2, - 24) k 

and 
Re = Sp-Q = Kp X) Et BMe- Yq i + 8 Me” 2) ® 
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(9) 


(10) 


The value of "'a'’ is then expressed as; 











—_ ae | eee ? , 
ut Re S ; Ke - Xq) Xe = 4? +BY, - Y J, - ¥5) * BZ, - 2.) (2 - 23) 

—_ — ; ; D 9 : 3 5 ; 

ss (X- Xs) + BY, Y.) + p (2, z.) 


(11) 


and the components of h by, 


2 2 
(X. - XIN. - X.) + We - YO, -Y,) 8 + Zp - 20 - 2.) 8B 
Pgh fog bd ee ae (Xp = X) 


h = (X,-X)- 
£~ “q ee: oor epaee £ 
X, - Xj) + (Y, - Y;) po + (2, a B 





2 | 2 
Keo Xq) Xp) + Me = Yh Ces Y) BT ply) p28 2 
rete es 
ee et (@, - 2.) 
Z "tT 2 £ of 


2 Z 2 Z | 
(X 7 x3) + B (Y, ‘ Y3) + p (2, 7 Z.) 


The magnitude of q is found by substituting 


on et che | 
x 4 Zz grin ee 





and the following expressions for cosa and cos@ into equation (2). 
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q 
le 
i 
f 
ie 


-Z) 











(X, - XQ, - Xt Biv, - YJ, - +B (2, - JG: q) 
cos 00 8 
a ee ae ee eee le ee oe ae: 2 
(Xe - X)0 + BW, - es BY(zZ, - 2.) (X, x. oP, © ue + BZ. - oa 
(13) 
2), > 24) 





ee stasis < & ae i 


coop = —— 
7 | | aaa: Oa tae ae aes 
(MND FB Ope Ne lee 2)) 











q _are then given by the multiplication of the 


The components of the vector 
components of equation (8) by |q |. 


a || 0 Yo) @_ > 24) - Ze - 2M > 5 | B° 


f | 
ee eee eg eee eee 


q = ae 
x | RXS | 
- 2 
af ee ees EN) 
x | RXS | 


ail 


q, = Seas 
| RXS | 
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where | 





ia 
WY 


ay ~ 2 
1 RXS | = 2 + x . _ 6s S 
p | o% Xp 2) > Ue - 2 %X))| 6 


ns VN. Se ; : 
bY, ¥ je 2) - Zp - 2M, 2] 





(15) 
He 


+ : ae, Y;) - (Y, : tae . x) |e" | 


i 


The velocity induced at a control point by a vortex segment is then given by 
equation (14). Since a curved vortex can be represented by a number of 
straight segments, this equation can be used to compute the induced flow pro- 


duced by a vortex of arbitrary shape. 


The components of velocity induced by a quadrilateral vortex can be 
written as ratios computed by the product of influence matrices and the — 
vortex strengths. 


























= cm | Ay Vix (16) 
V K 
— = A cancer 
2 = 2 and . 
ae ie W _ K 
Vo; ~ | ee Veo i 














where the elements of Ay, Ay, and A,, are computed from the following 
equations. | | 
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i 


7 p*| cos B - cosa | |g - Y) - Zz.) ae Ft 2), - YD) | rae 
Ay = ». Se me ge a rie eg ee es 
47 | h|| RXS | 
(19) 


. 7 B | cos B - cos a [Rep Set) - Ze2 XX] 
ae » —— a ee (20) 
Y 3 4n | h|] RXS| 


and 


B| cos B - cos oe | Xe - Xe soa Olea. * YQ) & - Xi) 





Ay aa >: a © — —- 
4n | h |] RXS | 
(21) 


The Z sign indicates that the contributions from all of the sides of the nee ot 
quadrilateral vortex are summed. 
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Appendix D 


WOODWARD'S DISTRIBUTED PANEL 
INFLUENCE EQUATIONS 


The equations will be derived for supersonic flow first, then for 
subsonic flow in subappendix A. 


Preliminaries 


Generalized potential function. If 


A N 








B og “S Qn, M55 
Then 
ee Veen): sn 2 |X oe + =e ds 
os & | {ce oO) 
where 
o = cosh” ——— 
p iCy-7) + AZ=E) 
and 


ie (-6 n,n), 75) > _%! 


with n= (nj), 7M, Nz) the unit normal to surface S 
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Vv 





Region of integration. - In the (£,7,¢) coordinate system the plane of Ne 
the semi-infinite triangular surface is determined by ¢ = a¢. The lines 7 = 0 
and n = mg are the projections in the ¢, n plane of the triangle edges. The 
area of integration in equation (1) lies on the semi-infinite triangle and is 
within the Mach forecone from the point (x, y, z) given by &<x and 
(x-£)4 > B“(y-n)4 + p4(z-t)2. The surface integral is carried out by inte- 
grating first over & and then over”. 


The ¢ integration goes from the leading edge & = §] = 7/m to the inter- 
section of the Mach forecone with the semi-infinite triangle, &€ = &€2(m), where 


Since {= a€, (x-&9)4 = 8 Z(y-n)2 + Bé(z-ak>) . The Nintegration goes from 
n = o to the intersection of the Mach forecone the leading edge where 7 = 3. 


Thus, 
N3\¢ 2 2 2 aa 
(x -—) a (y- ns) + £ (; tea oe 


Looking down on ¢ axis 







0 7] = ”, ewe 
(x,y,z) 
n 
* 7 = 
_ a MACH CONE 
\ 
(eh ep Gan & Ses: 
g 
Thus 
mee 2 
by are dee (2) 
2 7 : 
(x-e,f = 87 -m? +p” ( - aby)? (3) 


‘4154 


[ie INE 





| 


( 





2g 


(15 -ml = aiming =v + 8%(ang - my 2 


These relations may be manipulated into other forms. 
2 i 2 2 2 | 2 
Cig Ye (mee) |". Big eye Bala sy) (ay - mz)| 
OT 


. 2 2 
(ng -yP [2 - B2(a? +m] = 2[tm - y) + Belay - m2) ]z - ¥) 


2 2 Z 
+ (m - y)” - Bo (ay - mz)° = 0 


~ “Thus on the forecone 


2 a 
1 care (a° +m‘) eye (mx - y) + p‘a(ay - mz) 


rr Nth Nein ne OS 


: \/ Lem sae aes pa(ay - al - F - ror + m°)| | (mx : oe 2B Cay - ay 


The following relation may be used to manipulate this further. 


Z 
m(x - £B 45) - y( - p 2a’) = | (mx - y) + B afay - mz)| 


Z 
= p ’m (2 - sie ¢ (1 - B Es | (mx - a - p (ay 7 mz)* 


(5) 
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Thus Rey 


1 - ae + n\n, = m(x - Ba) - Boney 





(nx - y)7 + (2 - ax)? - B “(ay - mz)‘ 


(6) 
oA seer Zz 
Or starting from (4) again after multiplying by (1 - gp a) 
Zed, 2 22 Ze 2 
(1 - Boa )(n, - mx) = Bo m(1- Ba )(n, - y) 


4 


2 ie 
OBS Ll eB a-){a(m, ~ mx) - m(z - ax) 


which becomes 





22.4 3 2 22 
(1 - Boa’) (nz - mx)” + 28 am(z - ax)(1 - B a J(n, - mx) 
2 2 2 a0 29 a ti 
+ [pan] (2 - ax)? = pm} - Ba), - vy) + (2 - a) 
Or 
2 2 | 
(1- Boa), - mx) + B am(z - ax) 
? 
= £B male = p*a")(n, ; y) © 2 = ri 
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OF 


«2 ; 


n.(2 - B a -m(x - £B pe = one 


p m|(1- g “a*) (n, y) 


Surface Distribution of Sources 


Boundary conditions. - For a distribution of sources on the plane 
¢ = ag we will assume 











{t 
rN 
ow 
I 
C 
ad 
< 


= ee ta -u') + (w- 9] 


fee 


const 


— 
— 


WwW + Pas = 


fee 


The statement Woodward makes on the bottom of page 17, of reference (55), 
u' = -u and w' = -w is not true and not necessary. 


“57 





In equation (1), if we set Q =¢, then (a) says the second integral — 
vanishes and (b) means the quantity 00/av + 99'/av' may be removed from the 
integral. These assumptions may be checked after the integration is performed. 


+ (2 - ax)‘ 


(7) 


(8) 














Evaluation of the intesral over &. - The integral which results for a 
surface distribution of sources 15, 





; a3 59 
g(x, y, 2) * _wt+B au d§ dn 
TT 
(x Ey p(y -)° : B°(2 5 at) 
‘@) -] 
(9) 
Rut 
—_ de 
pe 2 
fe 8)" p(y <n)? - 87 (2 - ab) 


= p(X 3) = Bia (2-08) Hl -Bo") ((x- 8)°- B*(4-3)"- A C2-a 39] 
2 Y'-B*o* | 


(K- 8 )- lal e-a3) Vl Ae) L&- 89° at (q-g)*- Aaa 8) 
Which may be verified by differentiating 


<n: ase, OX 


A = (%-3) - Bra (2-23) 

Bt = B* (G~pra*Wq-9)" + a0] 

che Urptat)[ Ge ar*- attra - At (2-43) ] 
dt. Sa ae, 

ne BP 
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Se vetfay Ale Sdefiithe ttgueD we meet B50 Je. 


Q A+c Ab+c' A-c 
os eae Are Ac 
7 2(ac'-Ac) 


A*- a= 








a B* (4-7) - 2 (e-a8>° 





/ix-3)° 


QE. D. 





‘iw Since from Eq (3) + 


2 ( 
«ei? = Br yen) + Bo Gab © 


and from Eq (2) 





b> 
at 
ae aes ak eee 
ae)? - a? em)" - Bo @- ag) 
= | 
J 





(nax-F) - Bra (met-a7) + (-pa*){ cmr- 1) - pe (4-3)° - 3 (mit-a7) | 
ee ee 


| 
27i-pa* Awe, (amx-74) ~ fa (me - -a7) ~yi- ~ pat) | (mx: -4y* - pow (4-4) © B'(m2- a7)" | 








(10) 
ae 2 
Therefore from Eq (9) 
¢(x,y,2) = 
2 ~ : 
- Bi! : pra [ Sox (-mu- 4) - ~ppaltma- -o%)+ (i- pia) {(mx- 1)- fon (7- 3) -p “(me a7) J a7] 
on alma -a9)-¥ Ui pia) (ms eset q) eGo gy) 
pig AY at Bae , sa rye (-p'e*) (Car seg) 58 " 


(11) 
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Evaluation of the Intecral over 1. - To evaluate equation (11) we must 


a ee 


first integrate by parts. If we let, 








f =  (mx-7) - Ala (2-27) 














(12) 
gh = Btont [ U-ptat itary + (2-00) ] 
poe (~ ptat) [maT - Biet (y-a)*- pt (meraty J 
ae fis aie ht 
LA. oy ara Av = am 
rae - Ke a ray _ eae oo 
dw = eee es EN = ee 
= fw 1 h (4?-h*) “7 
a a CE 
hg Sa: 
rca ik Wa Cs a 


ha” 


iL 





At gg = fiw Cy -ptat) (qa) 


fo = = (1-p'a*) 





16] 





Now we can write 


> 


-n [gf - e'gf] 
naa) ea les oie Ge pe) = ye Pe) Cee) 


= pon? [m (x - paz) sae ae Oe pa’) (= pa”) (n- y)? (2 + ax) “] 


+ Bom y GQ - Ba’) | [m (x - p7az) - y (1- B’a")] ( n-y) | 


+ Ban? (z - ax)? In (1 - p 74%) - [m (x - p*az) OY Che pa’) | 


sf AB ne ly Gee sP ae) ya Pa) = Boa Cay # aa) 


+ p’n? a - 6202) | y me - o%a2) -y A - pa + Ce - ax)? | (4- y) 


2 222) 


22 
a’) - [me - Braz) -y GQ - Bad) | 


| + nr - ax)? ly a - 8 


And since from the integration by parts 








7S 
j £ th teh ~ f - gorf 
mi a = bee + ee os ED, 
. | Xoo, c-h 7 2] ree i, : d7 
BV aoe «2 
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ic 








ty * 
2 (-mx-7) - lal ma-a7) - 


qa pote ap os a (- ye (t- eer rae eee | 
eta Saks Es ceed oo 
-yeG- ax)*),/ In (x - @ 702) - nt - 6724)}7 - 6 a Bay 0- ne G- ax)'] . 





t yl mec SS 








= Bi (me2 09) ] 





(1-@a? )(enn-1)*- ate? (4-4) 
(i= pat) (Cr)? 384-9)? - at (m2-27)*] 








o*n [a p42) On ys + ax)?] 








2 
-y (l- Boa)} ie 


iaaa aid ee 








(A+B) +fAL AT +207 +¢] 





a - : eS ne 
2A a -(Aq+B) - fAL Aq  +2a7+C] 


(1S) 
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In our case expanding the denominator gives 


Ae (1- n?4?) {i - 2? (a? : mn’) ) 
ee 


: tz 
B= [-m (x - 0 7a:) + p?m?y] l- toa) 


2 
Ce wx: p7az) - on ty" (1 - aye G- ax)"] 


-(A4+8)= [ (aax-) - Ba (02-07) +P (1-3) ] (1-f'o*) 
A{atts2aqec]= (peti p latent] | [aGe-sar)- 1-6] = pie [(1-601-3) ae-40"T} 


= [atcatemnt)] (i-miet)” [ Gauep*- Blan? (t-a)'- Pi (ma 01) 








Therefore 
dn 
2 22% 22 22 2 2] 
m (x - ptaz)- 10- 620%] - pnt? [a- 67a) Cn-y)?+ @- wd ae, 


rat -t) - fb (ott-0.9) om (TS) fli-prtem (G02 1) Atle pep (meg) 


_ - | ¢ 
2 (1 -f Cateam™) } (i- pa?) Chimn-7) -$o.(m¥-a7)+pin'(4- 4) - fp <p (atone )}[ (wr 2 ‘ant (9-5) - ptm aa n} 





The last integral in equation (14) may be evaluated using the method 
described in Appendix A of reference (55). 


From Appendix A of reference 65); 


If 


22 


4 : 
Y + (c- ae)" - b%e = 0 (17) 
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ee 


then 


(A + By) dv 


5 
ig +e) a + 2bv +c 


Ab Y + BY> =4 Y Vay" + ¢by + 


a tan 
4 27 
Y +be re 


3 


ee 2. be? J 
eee ie Lag wy" + be’) + ve Viv’ + abv + « 
oe ay eye) -(v¥" 4 bet) = Ye av" 2bvee 


(18) 
































From equation (14) we can write 


3 | 
ee p7az)- 1 (1- pad) - p’n? [1 - Ba) (n- y+ (2 - ax) “] 
2-67) 1-67 @ em] Cn- 39° 


-2fm (x - B2az)-y - p7aQ) Cn - y) Ge Bn”) 


| _ 
+ [m (x - paz) -y O-- B27) - Baa (z - ae 


ao 





: 2 2 
Therefore if we factor cut a (1 - 8a’) we can write the denominator as 





, ee ? 2 2 
ra p2aty Cn yy? (2 ax) yf in (ane) - 9 88D] _ pm ta - Ba) Cn y)! © (2 - ax)*) 


3/2 mi 
s (1 - pea’) (v? + °°) \ jay" + 2bv +c 
(19) 


where 


ede: Wave Bo fer 


b = - [m {x - paz) ay a pa’)) 
2 -] 
c = [m (x - Baz) -y (1 - p°a°)) - ane (z - aor (i - ore 
-1/2 
e = (z- ax) (1- gia‘) 


Vo n- y 
Referring to equation (17) 


é ? 
a0! . im {x - paz aie ea 0 ee ewe) : Bn (Zz. ax)? - {i - p? fa +m) (2 - ax)” 





ct 
Qa - Ba”) 
2 
_ b a a1 ga) 
(1 - Ba) 
2 
a Ve Be 2 bree = et lone re (1 - let = 1g 
(1- Ba) 
Therefore we can choose 
y2 o ae a - p*a*) 
or 
Y == (2z- ax) 


(20) 
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Using this the numerator of the second integral in (12) becomes 





@ an? tly a 67s) - me - 6222) - Pit 28 ayy) Bethea ety | patmte a peg atta ga) eo anh? | py 5) 

































































= eat AE Aedes. geee a Lee. 2.2) | Keen 
— oS Se Therefor Ye. fr ee Lo the coefficient of tan” is 
; ghey kh oa bey 
: 74°) + b| * (1 - p*a*) ¥? - yb + y?| y(1 fa ) | Y . tt 
_ | 1- B°a’ 

















2 











z - ex 


























1 - 27,2 























ae - Ba’) 
“Using (19) and (18) ee the coefficient of =z = ; 
Z 2 2 
: “ty = p22) + bp + be stant 
: — SaCeecaanaae 372 
wake yg 0 “, 4 
a ar * ie (Zo ax) ( a + bé e ne ~ Ba’) 
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And from (19) 7 


(1 - pat ray? + 2bv +c] 


2 
= [m (x - Bap) - (vty) (1 - p*a°)] - pn [1 - rare + (Z - ny) 


ye - by = (2 - ax)? + vim (x - 622) -y (1 - B70] - (vy? + be?) 


| Pe 
-(va* + be) = | [m (x - ae. a eC ae ora - v(l- g 2a?) | Zax)” 


Therefore 


-1 ¥ Viv? + 2bv +c -} 
tan —— eer SF TAN 


Im (x = a2) - ety eee ie a 
Fa 


; (x - A%az} - nO - 6%s')) - 8*e [0 ye + (2 - ax)? 
- b yj 
Y +4 ges” “Nee: sak ax)’ + (my) [m (x - 8702) - y {I - 6 a’) 











~(v¥i4 bet) + YeVav +s 2bvece NY 
~(v¥"4 be*) - VYeyav'42bve+1e 


: Sos om(x-Bot)-Wi-Bot) + ea tey -P'aa)-BWi-A's*)) - - Bm? [Bat 1-3) ea") 


ea (Ro pee) A(t Bet )oe 7 [nor-ptan) -7(-p'0)]"- A [O-2e I 3) 18] 


= 0 (wr x- 7) - P'a(-m2-07) + ( i=Bia*) [ (-ex-1)"- Bn? (4-4) - Bt(~ma-a7) 
(mx) ~ 2a (-me-a7) — ff (1-ptet) L (-maeg)t~ Ae q-4)"- Bi (mee 7 )* I 
5 es oan (12) 
2 : 
[acx- tak) -4(1-prat)] - Bim® [G-p%a*) (4-0) +(2-ax¥ | 


= (1- Brat) Gau-4)"- A atty-g) - A (ma-ag)* J ae 
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Therefore combining terms 


c 





(ax-1) ~ Bio (2-07) + fU-B'a?) [Gn x-7)*- Bint (4-9)*= Bn a-0)? 


| Jy - 


(-mx-97) ~ %a (an2 “a4 ) — (1-p'a2) [ (mag) 4 - ln? (4-4)*- 3" (an2-a9)*] 


=! bean) - Bia (ma-o1) +7 U-Boat) [Orx-1)? - Bit (1-3)°- pm a-07)*) 
=3 (1-9) beta Raa sata ioe saad AS ttt ane latin dae it ll Gl Dae 


(om x-1) = B'a (om 2-09) - ¥ (1-4h') { (ma-4)? = Baat(q-a)*- 2) (— 2-07)"J 


4 2UiAlot) =n poe) 9 lnaty-palmzat)epel(ts)s y [1-Aat ne] Lone) Amt Crg) = Onze) 
2 G-2 ar) { L> p (at+ont) | (omx-4}) -Ba.(m 8-0) +Bm'(J-4)- ( 1- B*(a2enmt) }[ (onx-4)'-pimty.4)'- 2? (a12-09)*, 





Ce 9 \ : | ea 


, MEN ng [n(x - A “az) =n 8's) = 8 n't 284) Cn- y)? + - ax)?} 
22 23 3 _ 
ye es yi: Boa \(z - ax)? + (m- y) fm e-Bay -y 1 - p2a2))f 
(21) 


( 


Differentiation of the indefinite integral over 7 . 
ST Tatra ; 


- Equation (21) ma 
be checked by differentiation : a 


De fete Ade fat Leg Tee wre con nifin eck ts 
Ade Lnbagpatin Jeg gears oA ade Lktyal oor J. Ga Aefoe 


s = (am x-"7 ) - Bia (me-a7) 


re) 


SHH aiid AL 


a 


a= Pom f G- p'o* )Oy-x y" 4+(#-axy* ] 


< - Bat) } Can a-7)” = prow CA-4) - Gta | 


tH 
HAW TE 


Xan fo = gt ah 


a +h | | 
Ond — = 
| | zany $-h ha* 
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| 
re" 
os 
op 
oe 
eae 
Cc» 
CP 
J 
a 
|v. 
r 
~ on 
ac 





Ap ae 
f+ - (1-B'at ) 


43° 
(-ay- A lsa-aq). 


2 (m= 1) Ba (m2-a7) Hll-fa?) [re Bm 
i 2 (1-4) Jeg, a ek eee 
ae (mx-77) — Ba (ms-o7) -/(1-A'a?) [onn-q)- Bom (9-9) - A (m2 07) 
Gan) - Pia (mao) +/U-atat) [oan 1) fit (4-y)" -B Gne-4 7)" J 
(1- f'n?) (rn)? ten? (q-gy'- Atm ta)? | 


Pow 3 - pat) (7-3) 








i 
= 2 
(-nx-1) - Bia (am 8-07) - 
re ae 
- (7 y) (1 - p*a*y ze = 
vin cs - pay 9 (1- 8704)" - Bn ta - 84a) Cn y)? + (2 - and) 1a - 879") C0 yy? « (2 - ax)?) 
(22) 


— 


ee cee eee 


\ [m (x - paz) - a - oly: - gem (Ql - 8 2°) (1 - y) + (2 ax) ”] 


De Afhwtik Ale areond oy Tease x, TY, i 


fork sparen re fancinn, rnendegtaion 
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2 


fom (x= aia e) - omy [1- rat omtd 17 } 





2 
= | mix-plaa) - gii-pia*) . [1- Barend] 7-2) f 


= Bom (2-ax) + ()- Brat) { (mx-y) - ~ (an-me)™ | 


ae ag 4) - Bia (-z- -ay)|{ 1 ifs Cate) | (4-4) 








want =p’) [1 pet] (1-3) 








megs gee ieee low 


rea = 























ae = — “po ae . “ac Gm x ron ax + ee -oO x =e - - p Ke a x rm 2 ,; | 











































































































a fe “p ‘at +m a (an x=)” - 20-4) (%-q) + (4-4) 

















vp ee ee ee seyitt 3) * Oo = Cr 4)" | 


— foom* CI - 4)” + 


oe pint [macys Gant ptlag-m at] | 
+ [lpia] | ore p= pia-gr- A eee | 
Neue 22 ve Qt 


+ = wm (x~ Bia B) - pix 4 ~ [1-A* Co ai. 


—< 
=e 
a5 


Kee) at ANIL i] LP MN i 


Hees 


u 


a” = Bm [ Gax-4)* + (2- ax)” - Arcag- me" ] 
h 


c [1- p° Cats mt)| [ Gax- 7)" -pa (4-4) - ‘(me-agy® | 


171 


un s*s ateh” | _ (23) 








{ 


ee eee eS ese = 2 foe, oe (nang 139) 


day * f-h h 


fs ~ [ I- A Carsmt?s | 
ane 


j 


4d 
ea 


Gren 1) = Bolom 2-09) + Bnn'C7-3) +Vote) ome Og)" Bm tra) ] 


Come) ~ Ba Com a-09) 4pina! (pny) -V [1 Plate?) (eT Bimtlq-qy piim? 07) ] 





al ee p* (a? + n’)) 


| m (x B az) - pin’y - n [1 - pé (a? +m’) | ° : pin [(mx - y)‘ (2 - ax)? - a’ (ay - mx)*] 


: ~ Vn? ate mys cr - 870%) 


[m (x - p az) = O- p7a’y) , p °m? ((l-- p’a’) (n- y)? + (2 - ax) *] 





(24) 


er ifl 
( = 
Mg, 


This is because, using (5) 


? 
22 ; 
a9 thew p49 - ey a ee a? + ny | - p'm? ((mx -y)? + (2 ax)? - B (ay - mi)™] 


2 2 ” Z 4 + 
Sa we aah ae) et ah AL ee” Ge Se lb hw Bee ea 


2 
ez ae a ae 44 : 
- Bm [nm (x - Baz) -y tl - Bay] * Bom ot: - ayy” 


3 ? * 2 4 3 % 
2 |i - Gavel eee od 1 f¥-@aty) 2 fi ear) [n. te Fae) ey Cr p 7a") ee o. 


2 P 3 ? 
iS Pease a ea we ee a ee 


Bea? 38 ? 2 ee ; : 
- 1 88 gate ary) fim ix o7az) ya. B7a%)- a 82a) 0m - yy’) - oe CC - a’) On ye (2 - any’ 





3 
2 2 2 m4 22 ? ? > 
= (2+ p* (a* +n] ffm (x - a 7a2) - 9 G+ p7a4}) - 82m? 1» pe) Ca 9 * te an] | 


Loe com alae wute, 


". 

' 
ie 

rh 
NS) 
C 


( 





{z - axe 


—- 


2 
Z iz a: 2 2 2 
f = (z - ax) [Im (x - Baz) - 1 (1 - B a)] - Bm [(1 - Boa (Cn - y)” + (2 - ax)] 
l 2 22 2 22 ay | 
ff = - (z - ax) (i = pia ys eS a = OS BL ee CS 
/ 2) | 2 2 22 
¢ = ar 3 a (Zak) oto Mey) im (X= SB azy ey = g7021 | 
1 22 Z 22 
rate eae a> [mi (ee Baz) ey l= 83) ] 
‘ 2 
gff' - p'f 22 2 22 22 2 : ee - pa’) | 
= (1 - Boa) [m (x - Baz} - 9 (L- Ba) + Am (0- y)] | (2 - ax) + Cm y) [m (x - Boaz) -y O- Bad) 
(2 = ax)’ ven p74? , . 
im (x - B22) -y 1 - 070) | cx - #202) - vcr - 8202) - Bet #204 Cn- > C - 4y 
oe tea? eee ee a ee ae a es ee 
+ fm (x - fae} n (Ul - pa*)| Im (x - pan) -y (Ql - p 4%) ere - p *32) me ee ee pat) (n-y) d - p24) | 
Pas cee eee | eee ew OO Se eee Geen 
2.2 ? 2 Oe ok a. 2 | 
f +e" =. (2 <2) [fm (x - ®7azy - nn C01 Ray Caan | ae ea) fale (2 any ] 
Se. cA 2 Z ae Rea ae 
#41 Pa Te ES NO esas - AIO Pts. to SN Re ee Es 
[m (x - Baz) -y (1 B ay) ie ees (1 “f"5*) fl - pf” (a bane) - atm ee roe a 
= | : 22 | 4 2 zZ 2.2 2 | — 
+ (1 - Ba (z- ax) +On-y) [mf - Baz)-y Q BR “a’)] ; 
, - mh ; aa > : i. 7 Sak. ‘ = 
= (z - ax) mc (x - B “az) -y (1-8 “a)] + [1 - 8B (a  +m)] (z - ax) | Sets 
22 2 2 22.4 a ee: 2 
+ (1 - Ba} (7 - y) Briere. aces meee es: (a +m )} (z - ax) | a 
— : | a 7 rs ; ; > ; ; “ “ . = . 
y}) + (z - ax) | [m (x - Boaz) - y (1 - Bva)) + [1 - 8B (a +m)) (z - ax) | a 
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Therefore 


ee 








5 $$ Pel ~~ tritiation 
ae a os 79 7 ’ 

d ein (z-: ax) \) Ee paz) - nfl - A737 | - ion fe ea y ie y)" + (2- ax)? ] 

dyn i. 85 Fe ee a TCT pee pee ee ee Pe 


2. 


a. wr 7 
va oe ae ea oa *({n- yt Fe Braz) - yf) - pa?) | 





a < 


a ae aaa ce any aa? [ote — fan ae p72) | 


te pre ia a Sere eee eee eee 
adin- y) Ble. ae Ne Prefs Braz} -A(L - Boa)y + Bom |(1- B’a")(n- y) + (2+ ax)?| 








(25) 
Combining (22), (24), and (25) 









: 2 2 
2 201-3) dey me (i= Bat) [enx-) = Bim’ (J-3)" - am z-a 7) ] 
q amt-F) - Balma-a7 )~ ( 1 Blo?) Cam x-1) = Bimt (4-4)? 2 (an 2-e7)*] 






, me: ps!) nlite fia (K-71) ~ Ba (n2-01) +BtT-3) + [1B am I [Gna Bet(4-3)- Blom B09)" ] 
a (1-2 a2) fy ~ £2 Catan] van-1) ~Palma-a7)spoe'(4-3)- iC 


~ A '(atone)] ( Coee- 1) <pek7-3)= pm 2-a7)'] 





> j le ae 
2 e ax) _ tan”! {2- ax} V [atx -_ pant -"n() - p7°) | - pom t - eee yy" + (z- ax)" | 


rr rp a. 


yi Ba 1- Bra [az - ax)? + (ne y) jmtx - Bar) - y= @72)|| 





a ae Jog SD Balms-e Df -Bet Lene TT 2ie -4)"- 2 
2 


(nz-a07)" J 


Comme F) = Be Com a-at) — of C- pe) [ Cmmgy® Bimt(q-4)*- B'(mg-09)*] 


. ~C1-ayi- B'et) [Ci- Bia) -4)" + (z-0n)] 


8m 0? [mtx Bad - yr 872) Cm 90 7973] + [mez - 9asy - yer - 0724] [ar - 822) (m- 2» an? | 


1 
See 
@ = a? [mx - 64a) - 00 -8747)| | x ic Ba yn- yy? + Cr a] [mee - 8222) - nc - 0702)] - oe? [Cr - 82e2)(n- yy? @ Ct ax)’] | 


(ee x1) - Bla (mito) + V1 Ba) [Qnn- 7) Bem) Bln 8-07)" ] 
‘ (mrz- 1) - B'a Cm 9-09) - 4 Ci at) { Camxet)*- gla(4-9)*- 6? (on a-09 )*] 


i. 
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Evaluation of the definite integral. - From (6) and (7) 


m (x - paz) 4 Hig (ls pa’) = gm V(l - pra*)(n3 - y)? o{2 ax)? (6) 


eer | : 
m (x - plaz) - pemy - [L- B2(a2 + m)} nz = Bm (me - yj? + (2 - ax)? - pe(ay - mz)° (7) 
are zero when 7 = 7 
7} = 


Therefore wrens (0),(7) and (23) Lah of Ye Are, terms in equation (21) 
3" However the tan” 
73, depending upon }\ 


term can be either zero or +m when 
snether the denominator is greater or less than zero. 
If it is greater than zero the limit 1s zero aS 4 —- Nz. 
= N- 4 6 
Uy Ms watng (6) ond ARon (5) aonk(7), 
(2 - any? (mg - y) [mlx 


pea) - yl - p*a)) 
tr - eta? ¢ my) (z - ax)? + [ale 


Qy (21) when 


a + 


' 


. plaz) - yO - pra?) 
en ¥ (a " : 32 . . an 





Vim - yi? ¢ 


—d 


Le 


g2(a? « 





(2 - ax? -p2(ay - me)*y (mtx -_ptax) - ya - %a4y] 


Bz)" Smimtx: ear} -yit- a225}} 


(26) 


Now if we examine the numerator of (26) and use (5) 
(1 - pa 


2y2 [(mx - y)? + (2 : ax)? - B° (ay - mz)*} - Bene {im(x - paz) -y (1 
a 


: pea’))* 
e(1- peat) (1 g2(a2 + m°)) (rx - y)? - pay - mey2) + (Q- p2a2y? - p4n®) cz 





=f] - 2 ee | : 
fl bra + m*)) (1 - pea’) (ave a . G > ax) - p° (ay - mz)*} + Bene (z - ax)" 
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(27). 


because 
(1 2 grae)* = Ban’ = (1 7 pea“ z g‘m*) (1 a p2a“ + 2m?) 


But if 1- g(a’ +m) > 0 then from (27) 


(1 - pta2) \fom - y)? + (2 - ax)* - gt(ay - mz)? > pm | mx - gaz) - yl - g2a%)| 





and (26) shows that the denominator of the tan7! term is > 0. But if 
1 - B2(ae + m¢) < 0 then (27) shows 


(1 - B*a°) \ f (mx - y)é + (z : ax)" - B“(ay - mz)? < pm m{x - g*az) ae © ae Bra’) 


and the denominator of the tan-! term is > 0 if and only if 
m(x - Baz) -y{1 = Ba“) > 0 


But 


mx - Baz) - yQ - pa’) = mix - pVYy? +27] + pm Vy" + 22 - p2ame - ya - p%a*) 


But we can also write 


Bem (y2 + 24) ee [ B¢amz + y(1 = B%at)] 


7 Dvd 2 2 22 : 
[6° (aoe mm) = Ijy* 4 B¢m Zz - gem? pa2z2 : aaa ae Ba*) anyz + B%a* (1 . Bea2)y? 
= (A? (a? +m) - iy + pec - p2atycm?z? - 2 anys + ay) 


z Zp.2 Z : 
= [ B“(a* + m*°)- lly? + B°(1 - Beat) (mz - ay)? > 0 if 1 - B&(av +m) < 0 
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And therefore 


pm Vy? + 22 - ptamz - yl - p%a%) > 0 


and 
m(x - p4az) - y(1 - B¢a*) > 0 
a & 
te gp2(aé + m°) < 0 


X > B(y? + 24) 


Therefore inside the Mach cone from the origin x4 > B“(y? + 2°) we have 


(z - ax)? + (nz - y) [m(x - paz) - y(1- B*a*)] > 0 (28) 


and therefore because of (6) and (28) 


(ax) Vinx - paz) - m3 - pta®yy - pena - BY) (m5 - y)7 + G - ae] 


-i ge ee 
tan - _ et ag ie a = {) 


V1 - Béeaé | (2 - ax}? + C15 - yd} {m(x - ptuz) - yi = pra’) | 





If we change the sign of the denominator we would get a for z>ax and 
-n for z<ax. However when we evaluate the term when 7 = 0 and subtract, the 
result is the same. 
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Therefore if we evaluate (21) for n = 0 we get for (11) 





1. ad 
eae : os 
9 (x, y, 2) = pacer (ona 1)-Balon to] ek Bet) [CDP tgp Creo) | J 

2 , a - 2 a a oe aero 

anil - peat 1-1) Balm a-a9)-V-Ba")[ (mx) -A nly yep nto) ] 
; | 
= w+ Beat 2 - ax) tan m(z - ax) fx? - Be(y? + 2°) 
i L- p “at y[(y - mx) - pa (ay =z) |} {2 = ax)¢ 





F y(i - p2a*) - m(x - p2az) Slow x-Alazemg) + ¥[xt- Bat +2] (1-2 ae] 
2(1 - peat) V1 - B2(a® + m2) K-Blatsmgy - VfB sea Tevet] 


y Don So Pat sot piston eo 
2\/1 - p2aé X-Boa foe CE) -B'8*) 


(29) 


Evaluation of the velocity components. - The velocity components may now be 
obtained by differentiation of the velocity potential. However, whoever attenpts 
to take partial derivatives of (29) is in for a long days work. But the partial 
derivatives of (29) may be obtained quite simply if it is noted that all tems 

in the expressions for the velocity components must have coefficients involving 
either or tan7t. This is easily shown by differentiating the expression 


for g given by (11). 


178 


Dé Ws Bad 


= aa Ox “oar l-aa 





Ws pan 971s ) eT) ~ Bol ma-aty) +7 (1- Brat) [OrK- 1g) - Bee (Aye) -BiCr2-2 Ts). | 


















an YRARE | 3K Ray) polment)-firaD amylases] 
Is 
. | = Bow, - fla lm 2-07) +f{Q- Bo?) (emx- = Bom C13)" 6 (2-07) } 3 
: OK (amn-4)~ palma -07)—fC-plat) [Cary ple (q-3)'- at (a 2-07)" 
w+ pad 3 ae fia me ot _ : 
| cesses "yf. peat Nite eae: stan) Toa. Baye peel - pea(n - y)’ + ~ axle} — 
ae oe, 6 | oe 









































































































































/ a=, However. this Semel expression for u involves integrals of the same forn 
ee ae 1 013): (16), or (18). A closer examination shows that the resulting 
: : r_tan- 1 terms will have arguments which are identical with those for ¢, 








a aisle the ee may be different. Therefore all expressions must 
or tan™+ terms. The same arguments can be made for the 
































: Ly ro 





























see a a “With ies in earl Ae partial derivatives of @ can be obtained from (29) 
by ‘differentiating only the coefficients of the Lex, and tan"! terms. Differ- 
entiation of the tan-! or Log terms will not yield other tan71} or Joa, 


terms and therefore the sum of these terns must give zero. Therefore 
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WYN ULI 
AL ALTL | 
HO 
iy E 
: : 
elie 


| 
| 


: os ye “1 ___m(z = ax) Vx? - Baty? + 27) 


= ieee OO a Ee a tan 
= ax n(1 - p2a*) y{y - m) - ies na ss ae se 
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(| 
| 
| 
| 
i'l 
‘| Ie 
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= é\/l_- B26 92 + mi?) X~ AN ad tmAy-y er Ie -A* C9? = 
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| 
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Constant Pressure Surface 
Boundary Conditions. - In (1) we can set Q =u 


Assume that on S 





(a) | USF =" const 
du dul | _ OD) 
(0) ev api | © 
Now since 
Z 
OU) 8. es (OU ge. 
OV 5 88 > oO 
lta lta 
2 
dul _ _ = Pal ju, o-l ou 
ap! ; af 7 ar 
lt+a l+a 
and from (b) 
re) du’ : a 1 
SS + Spt ey -u') ae og ea) Ss. ap 
Ov ov Of af 
2 2 
lta lta 
Now from (a), u - u' = const on S implies that on S$ 
A: oe re; ae 
aT (u - u') 0 
and 
_o_ = ' _9o_ a t = 
ap u') taper tu u') 0 
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jel” 


This means that since 


Bas—(u -u') +> p Cu -u') = Q 




















ieee Herre? 








Op agit) ake eo) = Qons 


F 


= 


Note that assumptions (a) ‘and ion in + (30) are anieoen ene Scpeions. 


(a) Says that Gi yas constant on S 


(b) Effectively says that the "derivative of u - u' in a direction 
normal to S 1s zero. | 


Woodward assumes, in reference 65), that (a) implies (b), which is 
not true. _ 


If we examine the normal velocity on S 


o 
= 
tl 


gu. 


7" . Ww 
then since due to irrotationality 7 = aE 


moe = “= du, OW = epee ou 
d—g n a€ oF 3 d& do€ 
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4 i 
HB 

















and 





2g _ 2 9 hee Gs Oe da. 
lta je U, t 4’) = ae (u oe)  ) 0 


and therefore the source strength, if any, is constant. 


Evaluation of the integral over &. - The integral to be evaluated 1s, 





@(X,y,Z)= 


J iw 
(2 : at)é + fy -n) 7]. (z - ag)* + {z - ag) {c - ax) 

. a pe ee dgdn 
a+ yy =m" | (x -¢)°- @ Sts nee tae are 8 


"4 3 














Taq o ce Se) eS 
[cz - ag)? + ty -n)?| Cee ae p* (yn) - am (= - ag)? 
0 - . 
13 g, 
; (ee rar VAP ” dé F 
4maq 4 
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Except for the coefficient, the second integral is identical to the 
integral obtained for a surface distribution of sources. It has already 
been evaluated [see (9)]. The first integral over — 1s easy to evaluate, 


[tz - ax)(z- at) + (y- 1) Jag 


Lo - 09? + @- any] e892 - 8? 0)? - 82-08)? 


ae oi + ¥ (x3) B'[@-08)% (9-17) 
é &-3) = V Cx 5)* at lte-a3) "05-994 








This can be verified by differentiating, Sak , 


a = B [e-a3) Hly-4)7] 





Ns Oe pie awn 




















gee. ks ba Peat : 
ee oe Ss Jeg fh kg [ #'g - faq'] pie, C12) 
ee Slee ell ial ota _ 
A { | aacer 


29 ws» vee leaner 
(K-35) ~Y(x-3)*= B°f C2-a3) 44-4)" 1 


ten aay +72) ee 





(2 -ax)(z-at)+(y-ny? 


7? [a-aty?+y-2y] Va- ey? - 2? @-agy- rye’ 
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bee | AI sil 
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Since from (3) 





aa 
, ; ; ; aed 
w- b= Pl ae)? +9) 
(x-3,) +/@-3)- A @-a8ot+-t ‘ 
8 &- 3.) -(x-3,)*- ALC 2-03,) eg J 
and f 2) ae oe 
rom (2) Ey . 
b 
2 
(z - ax)(2- ag) + (y- 1) ii 
2 2 2 
ase ee eny Cae ae 
bi 
- = dey SED Venn) Aetna m9)" ed 
(amex =) = of Coma) pflaq- m2) aan (y-9)] 
therefore 


(32) 
$04,3,2) > ps fe EES ET 
ia : aE nce ) ea ga woes 

T) = Vkowx*- Bhetq-ay- AM me-o) 


a [ (omk-4) = Ba (maa) + fO-Ba?) (Cornet) Ban (4-9) “A (reo) re 
ora to ) QD brn) -Ba(-m2-07) - VOnpia!) [ome MB 7-g)™ Bone _ay)) 


The second integral was evaluated previously for a surface distribution of 
sources. 
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ta. Evaluation of the indefinite intceral over _» . - The first integral in 


(32) must first be integrated by parts. 


Qk 44> 


rrr - 7] 


t 2 u u 2 
om (a-7 + B Cmte aq) 


as 
rT! 


a . (on x-q) — Bomt(y-7) - A Cmn2-a7) 
rr. a 


























B° - (an - aes (n - ies 

















{e 2° - gat] = p? {tar - mz)? +m’ Cn - y)? + (mx - 1) [2 (an =m) em (1 -y}] 


| 
I 
UT ET 
il H EL sf) aa itl FOAM AME 4 THe UH! 


2 
a pol om (am - mz) (2 - ax) +m’ (0 - y) (mx - y)} 


aE 






i 
! 







iM 
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ey)" 


; | (a mn - mz) [ca n- mz) +a fm -7 )| em (n - ” les : : ae 2 3 











am. 


Mayet thee 
( ? 
nae 


7 ete” - ze't| . p° 





m (mx - y) - a (z- ax) | ne + fie | 2 (z - ax) - y (m- "| n 











4 > 
eae [« 1 - a" + ( 71 - v7] 
) ¥ 
2 2 
am? EF {mx - y) - a (z - ai (az + my) +™ (a? +m) E (= sax) = ¥ tie: | n 
+ a ———————— 
(a? + nm 


at Z 
aw E (mx - y} - a (2 - ae | in (y + =} : 


But 





‘ (mx - y) - a (z - ax)| (az + my) + ta + a" E (z - ax) - y (mm | 


[n (az + my) - y (2° + | (mx - y) + E oS + a - a faz + my) | (2 - 


ax} 
= -a (ay - mz) (me - y) - m (ay - mz) (2 - ax) 
= fay - es 


| 
the 


C 
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Therefore 


z/ Ler hon no) + Venn 1) = Pow" (I-9)"- Biman) A 
2 (-mx-) = (nx-11)*= plan! (4-y)*- (an e-07)* 





) + fom-1)*- Pon CT-g)- Bima-e)® 


Te Pee ea 
(omx-7) “eer e ia ~ 2" (ma-07 





? 
rl - 
—3 








a ae ora t {a ro Ah] ed 
| = Gat en’) aye e Aigh ons eB at Ua a 


wy f He hn - mz)” ¢ (my ¢ 32) [, (mx - y) - a (z - Gan] ls m ce + 24) i (ux - y) - a (2 - ax) | 
+ ae se ee dy 
: vow ian - nz)? om oe y)?| vc [en omy! + me in 1] 7 - : 


The first integral in (34) may be evaluated using (16) because, | a 


| [n a p’az) - 9 (1 - p72") | - pm’ fa - pa) (n-y)’ + - ax) | 


we ae pte’) | a 0)? - p? fan - mm? - 8 'rr (n - y)?| 


(35) 


Therefore using (16) and (38) _ =: = 


dm 











Yi 0)? ane? we ce yn? 


7 | 


2 fi-pcatenn*) (omx- 1) “Acme -3) ey [ina Co a es eer -04)?) 











(36) 


‘167i 


Se 
zi 


The last term in equation (34) must be changed to the form of equation (18) 
in order to be evaluated. We can write 


(a 7 Sale + a (n - vA = (ar + a - 2m (az + my)” + ii oe + oa 


2 a a ee ee 2) 
hi 4 a mieten) 7 eee 
= (a +m) 7} - + 
d Bs ae 
a *m a +m 
: 2 2 
2 2 m (az + my) m (ay - mz 
ey (ae FI ||) + 
2 2 
a +m (a +m) 
Now we introduce a change of variables. Let 
u = yn - Bie my) 
a +m 
Therefore 
2 2 Z Dt Oy ssi Gayest 

(a y - mz) + Th (7 - y) = (a +m ) u + (37) 

2 27 

(a +m) 
m (az + my) - mx (a +m’) 
n- mx = ut (38) 
2 2 
(a +m) 
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(7 - mx)* - p° (an - mz)° : an” (7 yy 


E : Be kei ea? (az + my) - m (a +m) " 


| eG + a 


2 





a aon (ay - mz)” (a 35 mn”) # n° a + my) - X (a + ny | 
Dna aes hen eects , ae ee ee Ee, 


| ae * a) 








(39) 


































































































(a* + m?) : (a’ : m) 
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ear (ay - mz) > Ca + n”) + - ea + my) - x (a bg m2) | 


2 
(a° + mn”) 





























(40) 
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and referring to (17) 








| uns 
2 nt Z 
2 2 2 | aay 
a, -m (ay - mz) +m |m (mx - y) - a (z - ax) 
c-ae = $$$ 
2 2 
(a +m) 
Z 
22 a Im (mx - y) - a (z - ax) | ne (ay - aie 
(aS +m) 
+ (Cc - ae) y? - Bre. 
Z 
2 
2 2 a ln Gx ey ae ad 
= y 2 _ may = mz) v2 @ lame) -a@- ao] = 0 
ee ie, os 
(a +m) (a +m) 


Therefore we can choose 


m (ay - mz) 


e (41) 
(a” + mn”) 
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The numerator of the second integral in equation (34) can be written, 


—_. 
—-. -—— 
——. .—— 


i +myy" Le a2 + mb oar roe nee ty" : 2*yp 


(a +m) 


2 


2 2 
to + nm’) y° - m (az + =) u+m (az + my) y- + Be tay EL 
(a +m) 


I 


| acer - m (az an u + vn (az + my) + b Ge ia | 


Now referring to (18) we can write 


we F (az + my) + b a + nm”) | 


Z Z 
a +m 


2 L 
a +m 


The fae + i in the denominator of A and B comes from (37) which is in 
- the denominator of (34). From (18) and (41) 





| RRP 






mM 


i a Be 


3 bY E (az + my) + b (a + mn?) | + | (a + a) y" - m (az + ny) | Y 





| 
& 
~? 
+ 
\o9) 
~< 





Sot a. ese ee — 


—= y) + be (at + m)(y* +b’) 


i 
hy r 4 


tt 
WK TUE Th 
1 


y(a +m) sm (ay - mz) 


(a? + m”) (a’ a m°) 


ii 


i 
il 
i 


Nae 


iol 


and la 
: 2 by, oa De 22.2 - 
ae + BbeY -¥ fy feed ey * sass +m)Y Serer te 


vy" P b“e” ee . m-)(-¥~ ‘ b°) 





-m (az + my) 
2 2 
a +m 


From (39) and (40) 


au’ o-ohurae = Ge = pe (a-n or ae pn (7 -y)" 


2 


2 2 
{a° +m) atm 


mim (mx - y) - a (z - ax) 47 4, tet 
y? - bu = v? | n GP = ao | P ia pe mg 


m [mn (m - y) - a (z - a n ae lage ~ 2Zamyz + n?2*) + ee la (z - ax) - m (mx | (az + my} 
= ne sree pa ASE VE = Sef rey} errr 


2 2 2 > 2 
{a +m} 





alae, a 5 | F (m - y) - a (z - 6) n+ m™ E (z - ax) - y (m - a} 


u+b= n - mx from (38) and (40). 
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Therefore using (18), and (41) we get from (34) 

















































































ae} 
2 fay my? + (ry + a2) [n (Mix> #37) aces erate rae aes es <2) (yo 2} 
mo 1 ee te nn at unas oe nnn din 
2 2 
. Pe 2 2 2 ee 2] 
[tan -mz) +m ner (am -mz)° +m (9 - y) 
ee ee ee ee eee ee ees mee cae yy? 
age oe : me Peres edi Santee teen eter ee Des 
oe ee ied Ga) ew ea yey) | 
= omg) (ORT) lomn = Bis Cay piemeap™ (42) 
= * al DET) = Cm) Blom? (peg) Bim are)” 
oT Ther refore Ube UE ne (36), ond (42) in (34) e 





1- dios ncaltodasast~ Lee mre t + Venn-7)" = mt (7-3)" = A (me -07)* 
_ ee ED = Cone) Blom? (> g)" = B Can a-09)® 











am ((atamy) -(arsandx] (wAx~“f) = Aa (nde) Pan (-3)+ 
2 (aren) bona) Ah am aa} )4Banl4-3)- ae tnt) [met peei(-9) ery 





. pe [tan - mz)? eine (1 - 7] | 


(mz - ay) 
+ BOE 8) gan”? (43) 
a+m erayion Gea yeas ee er) 


HH aeailebii 
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| 


Differentiation of the indefinite inteyral over 7. - From (33) 
Gant) = Boia) - Bi lmaeatye 


d (wax-1) + 
dn Comex 7) = (med) — aint (-g) 67 lm 3-07)” 


Cte +s ah ee a ee ee 


~~ 





(nx -nyé = pé | (an : mz)? + m (n - y)* 


(an - aye + a (n - y)? | 
(44) 


From (24) using (35) 


{ - Bita*eun? 5} [nx ~ 4)- Bin (I-4)'- A'lna-at)? J 


1d 9 Gan-2)-Balma-at)+Amlt-3) + ¥ LB ott Mnnet) = Bon (1-9) PR bonte 
(me R= 9)~ fia CB -0°7) + Bem (4-4) = of [1 BiCad samt) ] [Cen 9)"= Pond (4g) Bil z-0-7)] 


2dn 


; rs. g° Cart a 


Vm -1)?- 6 [Can - mz)? +m? (n - y)?| 


(45) 
ond 
figit 
ee, ee eff! - g'f 
an tan ry = ¥} = 7 
(=) + ] f (f£ ae ) 
g 
warn 
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om ee 
( a 





2 2 
(nx - 9 } p iehat aun ieee asus | 


pad 


g = [mim-y) - a 


- 


(46) 
ae) ear [entecan on Gee | 
ff' = (mz - ay)’ len mx) - pea (an mz) - pn? (a - y)| 
gi = laa tea ten | 
fo sg? = (mz - ay) Lon =m)? 8? [Can - may? +m Co oul 
| 2 (47) 
+ | [n (mx = y) =a z- ax} ems ( - ax) ~ y (ox y] | 
i aoe (mz - ay)’ {tan - mz) +m (7 yy? 


Z 


+ Im (2 = ax) + 8 Cm ayy (9 < 


+ [mm (me y) (9 -y) - (2 - ax) (an - m3) 


s ap? (mz - ay)? ke : mx)? +m (Cn - y)?] 


+ [a2 (mx - yy? + 2 am (m - y) (2 - ax) +m Pee mx)? 


+ m? (mx ~y)? Cn - y)” = 2m (me ~ y) (z - ax) (n - y) (an 


mz) + (2 - anys (an ae 
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7 


« -p! (mz - ay)? [tan - mz)? +n (n i 


2 2 
+ (mx - yy? [tan - mz) +m (Zz ax) * n? {z= ax)? fe - yy > (me - y) | 
- 2m (mm - y) (2 - ax} [cr > y} - {mx yy | [tan mx) +m (2 ax) | 


cmttm syle tne yy? = Be (mx oy) 7 ax) Oy) Gam me) 6 fz = a) (a7 - mz)? 
= [im - yy? + 2 - any’ - PMiscaeng eee Ce 
[iis . Ferme ere ee as lern eh ee A ee a ace (oleae |) 
. Ae rr re ae ptlacn sce el acenelel 
: {[. fx -y)-a(z- ax} Cn Jae eee Sees | [1 pat ew) ] On me rae (2 ay mow - id | 
. [a tea mtx) (eo mnt fa 0? ew] s zera[a aay eneea) cme Be Le a ee 


Z 
= penetenqD: ata? sm] [im - y)2 + (2 ahs a? late ax) eee 


=om(n - mm} Lim «y+ @- an? - p?[m? (2 - ax)? + 2 am (z - ax! re aren 


= m( 7 - mm) Lim - yy? + (2 - ax)’ - p? ay - mx)’ | 
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} 


| 


l 
| 


therefore 


qo -y *Y 7 mz) 
aa tan ? 
dn 


2 
Grxcmy +o? lanes ee f4e97 | 





- (NE se) esa = ax) | i mn [2 (2 - ax) - y (mx - y) | 








| (an _ mz)? +m (nm - yy? | C= mx)” - pe (an m2) +m? (1 - yy? | 


(48) 














Therefore combining terms in (43), (46), (47), and (48) 


el aoe | amet) + PORK 1): Boe (I-9)"- Bn2-at)” 


> (net) = Coe) = Bhat (Tg) Cm Beet)" 








A” = 





ees ae f “an 1)4 flomax)" - Bigg) - Bi(me-09)* 
ee, + Comn= 9) = Cmca = ple (ye gy = OPC a0)? 


| mi faicsine) lyon -2e-aolen|s oe . -y} 
Peer We ke etree) te aes 


é fo Soo ee ee 
f= ae) ene { rey | Vin mo? - ag San men Ce gia 





aemnm 


2 2 
aflfaz em) - fale mn) x} 


(a? + m4) Yj (nx “9 Y - po [ta WARE Sk Me én 


a (ay - nz)? (1 - mJ 


a? ea) [G0 gone} ont (4% YI" \y (ex ay - a? liane my em ( o-¥)"] 
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Arranging all terms over a common denominator 


a Cant- 1) + Ylimx-4) = m'(q- 34) ~ 8 lem2-any” ‘ iis | 
la%s 


CO Ble) (me Bln 4-4)*- 8 (m2-ay)” wm) [04-0 8) 3 9)" ] Yew 1)= Pm 1-3) B tea 7)" 





x {ateant) J mlademg) }] om (1 lmx-3) - 2(2-ox)] + [a (B-ax) --m(onx-4)] 1} 


~ [ A(2-ax) - an(omrx-4)] (a4--mz) + an (4-5) | +n (ay-n2) (7 -x) 


The numerator can be written 
ROME eae a derrag: eel a alee oe ha ee a 
+ on{ -m (az + my) a (2 - ax) - mom - yy] + Ga? emt) mly me -y) - 2 @ - a0)] 
» 2fac-axy- mex -y)] [rz + mm em Cay - me)?| 
mn’ (az + my) Ly (x - y) - 2 @- ax) - [2 (ax) - meme - yy] G+ 22) - mex (ay - 2)? 
2 ain) eee ee eae ee lan a aaa nz) | 
+ nx - y) |m (az + my) + (2 + m2y - 2m [a2 + my] - 8 (ay - mz) | | 


+ ac | (z - ax) E (az + my) - a ty? + 7*) + mx (ay - nz) | 


+ (mx - S15 (az + my) +m y? ‘ 2*) + ax (ay - mx) | 
Tae eo ee ee 


ee | eee eae aly = 0 


Therefore (43) is verified. 
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Evaluation of the definite integral over_n. - Because of (4) and (6) and 
analogous to the section in the evaluation of the definite integral for the 
surface distribution of suurces, all of the temns in (43) will be zero 
when 9 = 72 provided the denominator of the tan! term is greater than zero. 
At n= 2, using {f), snd then (5), 


(n (me - y) - a (2 ax) |, . mh (z - ax} - y (mx 29 





ee er ee rm re i a me ee er re 





Tt [n (im - y) - a (2: 2% | [x - a *az) : pny : ay (mx - ic + (z - ax)? a (ay - mx)” | 





2 | E ix * p*az) -y 0 p2a°)| (mx - y} + la : pin?)2 - a (x - pny) | (z - ax) 


es | | E a (a )| 


the a-rim-o] 


ax) | Yiu ys (2 ax)’ - p? (ay - m)" If Be (a's *)| 


Ye ee se 





t4 


mfim ys Ga? BF tay my? ¢§ p———-——5- > 5 | 
ce, a ee | (mx - y} + {z : ax) - B (ay * mz) - B [s (mx ~ y) ~ a (= 7 | | 
1- B° (a +m } 


| But 


fe ; a se ax)* - p° (ay - Aye . a? le {mx - y) - a (z - a 
22 2 


= [3 - p? fa” + i | (rae - y)° + (z - ie + gia’ (mx - y) + gem’ (z - a 


(50) 


a, © 


+ 2 ae a (mx - y) (z - ax) - a? (ay - og 


2 


. I - Bp? (a's a | [mx - yy’ ae ; ax)? | . a’ E Mets eae a (z - ax) | - pg? (ay - az)? 


2 E - p? (a° + ny te - ic + (2 - ty 
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Therefore if l - ae (ar + me) > 0 


(mx - yy? + (2 - ax)? - B? (ay - mz)’ > Bm (m - y) - a (2 - ax)| 


and from (49) 
[nm - y)-a(@- adn, + m[z (@ - ax) -y (me -y)| > 0 
However also, 
2 De. 
m (mx - y) - a (z - ax) = f(a +m )x - (my + az) 
and 


2 2 4 
2 oe! 2 2 
(a +m)x - (my +az) = a + mo) x - fae + a) ty’ + 2°) © clays mz)‘ 


= x tat [f 92 (a” + 12) -1 foe + - - p? (y+ 2) | | + (ay - mz)° > 0 
p 


if 
1- p? ame < 0 
and 


xo o> pt yrs zy 
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Therefore if x > 0 also 
fa" : nx) > (my © az] 
and 
im (mx - y) "a (2 - ax) > 0 
if 
2g 2%) 


ei ae et | Gis 3 


— Therefore from (50) if 1 - B* (a” 7 mn’) < 0 


| 2? Gen <i | [m (mx = y) Bee an) ee eels ae faye ne) | >o 


iat Hie ik I! 


and from (49), if - > pe (y" + Z 


‘) 


AAV PAU TAL 


il PH 


[m (mx - y) - a (2 - ax) ng +m iy ee Ge eee > 0 


and therefore, using (4) and (54), the an term in (45) is zero when 


n= 7. 


201 





Therefore substituting n = 0 into (43), and from (32), 














7] 
oP > 
$(x%,4,2) = tee hee (mx--) +f (nx) = Bim’ (4-3) (a7-m2)™ | 
8 Tra Com n- 4) = (mrx- 7)" = phn? (7-3 ye ~ (07-2) 
© 
[1 ~Biloten?)} [ (max -4)- Am'(1-3)- al a J--m 2 | Af 
[1- Baten )){ axa pontlg-g)® 6a q-ant) 
se tls, - SO os eee) xX + x*- Cys 2?) 
47a 44 2 (atta?) KV Kt BC yte2?) 
_ an {taza gam) - (ateant)x | p X= A (oz+mg) + ¥ [!-B\carent)) | x*~ B'(gte2?)] 
2 (0. * 4-0? ) I- B* latent) X- Ba Z+my4) - V [1-° Caton) ( x% B°(yre2'y] 
Fe ~~". Cay-ma) om =| (an 2-and) ¥ x *. B's 22) 7 


O02 +an? 


Zl2@-ax) - A (oan x- 9) 


+- (Z-ax) y ee, om (B-axd yx & (47427) 
3 [ lye) -Ba(ag-m2d) *a-a4)” 


_ mlx. Baa) - 30-44") x- B(azimy)s V (1- Parr) x2 (ge 8)] 
| 2 fl ~ 0% lat 4m") X- 2° (a@ten4)- y (1- BCom?) {x 24742") ] 


-~ +4 yi- pot Lew CC YS A ed ed Le 
: (x-@lo®) - (1- pla?) [x*- ai(yte2') J 
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Therefore combining terms 


lL [AP 
O.X5yse) me a) 


Aen Jog — + +x? -f “Cy? + 32) 
2 Labo Ve gee) 


2 (a*ennt 


A a = [1-8* (avant) | (x" aliyt+27)J 
=] (mz - ay) Vx" - Rept tye Ba) 2, 22 


_ (mz _-_ay) 
~ 2(z - ax) - y(mx - ice we 


2 ys 
a +m 


tan 


| oa c ) cae m(z - ax) om: p* ty’ +2) 


y|or - mx} - B eas - ~ na)] 5 a a ss 


a Hota cane : a : ae Boa k- Baz + Vf (i- Bat) [x*- A cyte2®) | 
ry een nn ene. Aan AUT ianoonve nO ma S  eSS 
= | TO x pla - So-aby fat pti] 





Evaluation of the velocity components, - Using the arguments presented in 
the section covering the evaluation of the velocity components for the source 
distribution, the partial] derivatives of (51) may ne obtained by differen- 
tiating only the coefficients of cach term. - = 





Modifications and Regions of Validity of the Velocity Potential Functions 


Region ‘of validity Se ¢. - All of the integrations were performed without 
regard to the existence of the limits of integration, negative square rcots, 
etc. Therefore the functions in (29) and (51) must be examined and possibly 


modified for some regions of (x, y, 2) space, 
Orn DR ade tome oct, YESS he ol 
her $(x,3,%) Qadir Dey (51) - os eld sess x > B'(ate2), 
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alaryeme) Vil- “| te? tan’) f X- A lattmy) + Vv [)-B'Ca%san?)) [ x= ArCsten)] 


LY |= 3 Cotem*) < oOo C prgperronce Sooding adee ) 3 
Abe Nog Hine whch cxcttare ADe sper rvot of he 
qurentthig mast he rod fad . Los eS 
ee Se 
7 Rta t 


one Peado 23 1 B*Cat+m?) 


(52) 


< O 


ZVI A emt) Dog 


X-8'(a 24mg) + V[1- Bate?) ] [ x*- A" 3*+2?)] 

x - 2'(az +™3) -~yl I= *(a*4-mt) ] ( x? - Bi Cate 2°) J 

. -1 [A Cast) -1) Cx*- A*(at+ ard] : (53) 
far, 2. pe Se 

a or | X~ "(a B+my) 


2} l= B'Catem?) = 0 Beg na a. Pica Leeda adaa. 
ph pole CoA ke b= "Car ean?) —- © WAR Can WwAAg, 


7 X- Slee tmyh 4 VL 1-2 (o.tymt)| { x?- BC yt+ 27) ) 
s 3 z = 
2 a = ee Xow X ~ "(ad emy) - (1-A*Catemt)) [ x7- 2°( 4242 ) 
= Blend) $0 
[xt Bt(yt+ 2?) (54) 
3 { 1B (arent) Te ee a 
X- 3° (a2 +-m4) 


C 
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(iss ita ait : 








ua 


Ht) 


| 


T 


UL Sa 


PN || TUM 


ua 


TT 


mT 


PA 


fee | 


TERA TT AP al 


Pe ee 
- 


RT EAA 


I ea 


ay 


Al 


[meee 


ch Da og Pane agen Ae Sly oon 





se nde BALL Boag oe A Sle mngeeect onl 
wane genant aaeah gtr od aafoe amnct Ao 
Jaen Donan gettine A Dany agit 

Bae tan facie oaenmyacidl, ve a Be 


ei). a en m(z - ax) a EBC O. +z) ; 


y[(y - mx) - p‘a(ay some) F(z ax) 








(\- ats A (aqtea®)) 
2) ]( x? ~ p*( 4t+ 2*)] 





F2 = X- B*Camyeak) + 


2V1 - ?(a’ + a) d —— 





{[1-*la* Sant 





(56) 


“1 (mz - ay) Vx" - p°y* foes) 


oe. ee “3(z - ax) - y(mx - y 


Vo-paa Lee acatea)y | 


C1-@*o ) txt = ne Cy" +2) 





7 
= Vi- pa = X-foa = 


aos of the tan”! functions Fl and F3. - The terms Fl or F3 are 
always real inside the Mach cone from the origin, : 


x" - B’ ty? +24) ty 


| However as the argument of these functions go to zero the ‘Guetta may take 
on different velues depending upon how zee is approached. Corresponding to 
' each of the four quadrants, if 


then 


f>0O ¢g>0 means 0s 6 <= 


f20 g< 0 means 


f <0 g>0Q0O means 02865 “5 


= 


2<6<T 


T 


v 


f <0 g<0Q- means "3? @2 -7 


Therefore from (52) and (57) as z—~ ax 


0 
lim Fl = 

Ze AX 
-T 

and aS mMzZ—-eay 

0 
Lam BS_t=* 0 

mzZ-« ay 
=i. 
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y <0 or y > m 


O<y-< mx, Z > ax 


O< y< mx, Z > ax 


y <0 or y > m 


O<y<mx mz > ay 


O<y<mx m < ay 


(57) 


(58) 


(59) 


Vsbeok. apkankee’) oe hs * (at+2") Ag 1- Bi Cate ant) >O, Rafarra 

Re (52) FU and FO on ator ad x*. Care 2t) = 0. Sd. Junctions Fl ond 
F3 a TA aero 4 ord mly =A Abe Aanenrnitena 0h Aber Ss OAQ, 
Opnaken Aen seroL ree (6797. Sok we ate och pQDew ns 


x ppt img eazy}? = fx? - 87? 27g) + 82g? + 27) - aimty? + dams + 272”) 


ee a a ee Bh 
a © xt = 0G? 6 27) + a2 - 87 (a? + my) (y? 2) eB fay? 2am + 0727) 


: © xt 82? 6 2*y) oe 2t1 - p2ta? omy}? © 27) © Bi fay - mz)? 


[X- AC-my soe] f x + Play+ea 








seen 2 : = Q ° Xx > | 2 ¢ an ata?) 


ee — | ah x *- i r “+ z*) 20 1- 8° Cat sant) 7° 


=———_- OT 
=—__——— 


For Fl the denominator is 
| 2 2 
y ([(y- mx) - BYa Cay - mz)]) + (2 - ax) 


a y fie 6" (a” + m9) - my [x - B* (az + my)) + (2 - ax)’ 


> Oif y < 0 from (60) 


207 


y C0y- mx) - 67a (ay - mz)] + (2 - ax}? 
= (mx - y)? + (2 - ax)’ - p? (ay - mz)? + mx (y - mx) - B mz ( ay - mz) 
= (mx - y)” + (z- mi pe (ay - mz)? + mxy - mé (x? = B 2%) B “amyz 
= (mx - y)" + (2 - ax)? - 6? (ay - mz)? + myx - B2 (az + my) 
- mex2- p2 gy + 2°91 Lf 


if x? - p? (y’ + 4) = 0, y > 0 due to (50) and (60) 


Therefore 
y (Cy - mx) - Ba (ay - mz)) + (z - ax)” > 0 
on 
) a es a. 
x" = BY (y* + 2%) (61) 
cat 


1 RG ay eG 
Therefore due to (57) and (61) | 
lim Fl = 0 
as x2 - g2 (y2 + 72) +9 
if 


se p? fa” +m’) > 0 
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ade oe 








or FS 


(2 - ax)-y (m-y) = y2 + 27 - x (az + my) (62) 


= - rx! = p 2 (y? + 2 W355 + ry (x - Bp? (az + my)] > 0 


ane - g2 (y2 + 22) = 0, and 1 - g2 (a2 + m2 >0 due to (S4c) if x > 0 


oo 











t ear a surface distribution of 




















aed Se = 4 he e r e fo Ore "ex ani ; n in a f 7 9) A (Se - tae 10WS “a hel 











sources « ‘or a 2 constant peeeeer surface 














if 


‘ iat ¢@ will be continuous, as it must be, 





oe a 


as x2 = Be (y2 + z*) — {) 


1 - pt (a2 +m) —0 | (63) 


if we define @ to be zero 


outside the Mach cone from the origin. 


if 
p= B2 (ate wt) > 0 and x’ < BG +2) (64) 
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Supersonic leading edge and the mach cone envelope. With a supersonic 
leading edge, 1 - Be (a2 + m“) < 0, all functions in the equations for 9% will 
be shown to go to zero for points on the Mach cone from the origin, except for 
a region on this Mach cone which borders the envelope of Mach cones from the 
supersonic leading edge. Inside this envelope of Mach cones the functions Fl, 
-F2 and F3 will be shown to have constant values. However on the outer 
boundary of this envelope of Mach cones 9% will go to zero, and therefore all 
fimctions may be defined to be zero outside this envelope. 


The envelope of Mach cones from the leading edge is illustrated on p-29 
of Reference (55). The Mach cone from any point Xo? Yor 2% can be written, 


«-x) = B ty-y)? + @- 2973 


di J = d = 
on the leading el’ge MX, Y, and ax, Zo 
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C 





Fe pe 





tT Al 





UM 











{ 


Tey 
ah 


Therefore 
2 
(mx - y 7 = B°Cm (yy)? + (mz - ay.)*} (65) 


The Mach cone envelope is determined by the maximum values for z, at a 
given x and y, obtained by a variation of y.. Therefore differentiating (65) 
with respect to Yo? holding x and y constant, and setting dz/dy = 0 


gives 
es 
(m-y) = Bm (y-y,) + Boa mz - ay) 
or 


y,{1- 6? (a? +m)] = m[x- 8° (my + az) (66) 


Pe eo 0 (65) and (66) give 


eG = p? ty? + 2°) and x - B° (my + az) = 0O (67) 


or 
Z 


2 
ie pris | ae | 
B“a 


22 
Boa mG 


ptaty? + x2 - 2% my + Bony” 


at (a2 + m2) Pn kath a- ptt) = 0 
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mta \/B° Peabo] (68) 


where 


y = =~ mt+a B° (a + m) = oh ZU 
B° (a +m) 


For y, > 0 we can eliminate se from (65) and (66). 


From (65) 
F ape (a* + m*) |»2 an| x - 8 carta |y, +02]? - p*(y* + 2°) = 0 
OT 


Yo 


fa- B- a? +n’) | - nx p° an my) | 


+ ™ & ar bern ois B° a? + my | ies pe cy? + 225] 


Therefore from (66) 
x - B° (az + my)|* - 1 - p? ( a’ + m’)| x! - p° Cy" + 2°)| = Q 
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or 


ge (at +m) x? -2x 8° (az + my) - Bp! (ay - mz)? + Be Cy + 2 


i 
© 


and solving for x 


en ee re ee ee Ce eee | 


(a-+m) 








(az + my) + (ay - mee | p° a + a = i| 
é 
(a° +m) 








ee _ (az + my) +|mz - ay] Ve (av +m) - 1 


ence Co + nm”) 











wil If we write 


pee (a” - - 














Og Se. —- OT 
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(a + m*) E = Bn? | 


= (- a + m \/ 3? (a +m) - 1) eer \/ 87 Cat a 1) 
2 @ ea Fs Bam, a Nias Be egos 1 


Therefore use. 


+ for mz > ay 


- for mz < ay 


or 
eo Lat ny) + [mz + ay | Vp? @2 +n - 1 
(a° +m’) 
This is presented on page 29 of reference (55). 
If we set 
zZ = y tan 6 on x? = BG ees) 
then 
x? = gy" (1 + tan’ 6) 
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(69) 





or 


cosé@ 


ee #56 oe By = xX cosé (70) 


Therefore (67) gives for the two points where the envelope of Mach cones 
from the leading edge is on the Mach cone from the origin, 


yy? + 27) = 8? (my + az)" 
Or 
E + tan26 | = pm + 2 pam tan 6 + B*at tan’ 6 


which means 


a1 - pa?) tan = Bam+ \/pia’m?- CG - Ba’) A - g “m*) 























or 
ae pte?) tan o = ane V8? @? om!) - 1 cm) 
but 
, eg dee ; eee 
x - f° (my + az) = 0 
means 
ee ; = 
Taal " 





Therefore the two points corresponding to (68) are 


= pa : ey 2 Deo 
oan oT cos @ pew) ; 




















l 
(73) 
z Ba oe Zit 3 
tan 9, cos 8, + VB“ (a +m) - 1 
(assume a > 0 and 05 > 64) 
From (56) and (57) the value of Fl or F3 as we approach the Mach cone 
from the origin depends on the sign of the denominator of its argument. 
For 
F 1 on <5 = Be ty" + 2*) we can write 
1 Z 2 
Dl = >} y & - m) - Bia (ay - mz) + (2 ~ ax) 
y 
Z 
“ ~ gen) . Bm M _Ba_ 
= (1 Ba’) BB + p7an tan + [tan ¢ ba 
2 
reas _ pe 7,2 2 pa ; : 
1 B™ (a +n?) + [tan @ 2 an|an+ pa tan @ | 


= OQif @=6,,or 6 6. from (72) and (73) 
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If we differentiate this with respect to @ we get 


an 2 [tan re Be. | | sec”e : fase + Bm pa sec’ 9 yet | 


Z 
cos @ cos 9 





- 1 _ a _ - : 
= ; a G ae | | 2 (1 ga sin @) - gm cos @| 
cos 86 








= ! ta - oe | 2 (1 - Ba sin @ - gmcos 6) + gm cos6| 


é 
cos 9 cos 6 


From (72) 


1- gasing - pgmcos@é = 0 at 6 O4> 6» 








7 


and since 


cos @), >0, cos@ , >O and gm> 0 


2 


and using (73) we get 


aDl : 

cv} < 0 6 = @ 1 

dg. 2 Se 
217 


and therefore using (57) if 
ree, 
1 - p’ a” +n’) < 0 


: 0 Q@>85 8 
lim ? <0] 


ee ae a ee Sa (74) 
x—~-B (f(y +2) Tm Z<2xX %)<@< 85 
For 
F3 on x? = pe (y? + 2°) 
Let 


D3 = z (z - ax) - y (mx -- y) 


= y (1 + tan’ 6) - yx (m+ a tan @) 


2 
1p cose m+atme)| 


Z 
= ehh. (gmcos@ + ga sin 25] 


= 0if 9 04 or 6= 8, from (72) 
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_ me 74 D3 = 3 |m sin 6 - a cos@ 


Ze 
= 2 wit. |, (1 - Bia’) tang -a(l- p*."| 
pl- pa) 


Z 
2 eee 2 F “n° -a (1 - pa) +m \ |B? (a” +m) 1 


pl - Ba’) 


= _ x c0s 8 ; ali- o? @t+my Jem Ve? @2+ 2 3] 


pa- Ba) | 


, Bcose Vel on’) «1 |: Va? Gal wat 120 | 


as 64) 











a 3 >0 for 4, (+m) [see (73)] 





: = <0 for 6, (- m) (assuming a > 0) 
because 
nm >a? |p? @? +m -1| 


since 


44 phe ba j djetecd al 
PIM | 


(lie H a) 


(a” + mn”) (1 - p72”) > 0 


id! 


' 


ehet 
LUI 
i: 
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“1 Ride 


Therefore since 


Re D3 >0 6 62 
d@ <0 6 6] 
and therefore 
: >0 8@>6>5 or 0<86 
D3 <0 2:6 ep oo) 
1 2 
We can say 
0 @>0, 6<6 
| 2 1 
) ; Pee y AS me MmZ>ay 6,<0< 6 (76) 
| x ~—-B (y' +z) 0 MZ<ay 6,;<6@<65 
Also on the plane 
mz = ay 


2 2 
D3 = 2 (2- ax) -y (my) = y* fa + 4 sy [ao 


= y Q-m) (a +m) 
a 


and therefore 


0 y<0 y >mx 
F3 = wf OQ<y<mx mz>ay (77) 
=m QO<y<mx mz <ay 


lim 
(ay-mz)->0 
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2 Z ea Z 
To find the limit of F2 as aR (y +z) when 1- B° (a’ +m <2 
C [see (53) for F2] we note that . 


xX - p? (my + az) 


" 
Jnod 


lim 
— p? Cy? + 2) B \ ny)? + (z-ax)” 3 ee (ay-mz)* 


Areca 


4 22 2 422 
B° {erm = y)? + (2 - ax)? - p?fay ~ mz) | = Br (az « myx? + B27 « 22) - axn*tmy + az) © 2A amyz - B’a“y - pms 


= tp? (a* ? n’) - 1J x? ° od Co ¢ 2*)} + x? - 2xp° (my + az) +2 B  anyz + pim*y” + p*a*2* 


© ptta’ «my - 3) [xt - B2ty? © 249) + [x - a2 Gey ¢ a2)? 













































































on fl - ar ~ | = 82. 
x [1 - gmcos 6 - ga sin 6] ~~ D3 


Therefore from (75) 


5 
La 9 
@ 
‘os ae 
77 
pol 
wn 
Ou 
@ 
: 
© 
Qu 
© 
cs 
“1 
vt 
; 


ee rr ee ee eee Oe qe ee 








ee 2 | - (78 
, ae 2 F2 ot ous = ate iets Shs < 8 >0) ( ) 
~~ + Z cd | thee a 

- pee Eh pe 01< 6 <6, 
i | ip (a + mj)-l 

| ee eae Seamer . ee se ee. A ee Paha tetwisat aes ane aoe 


a ee ee ee a a Ar a ie en cea A cg ir ye ee 6 aia se me a gg ne ey. on) 


herefore c combining these results 


Fl tan’? 2 es) Vx? - 6? g? 2?) Vx - PO +2) x* > a? y? + 2°) 
a 7, mx) - p° a (ay-mz)] + (o-ax)" 


6 gee eee oe ee ee ee eee 


= ——_—e= SP ee ee i ee te 
*. 
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If 1- gp? (av +m)>0 


() 


2 
Fl = 0 x? #8" (y* + 2°) (79) 
: 0 y<0 y >mx 
ne Fl = 7 O<y<m Z>ax 
ia T O0<y<m Z< ax 


Then from (69) and (74) if 


and 
x2 ‘ B° Ge m m°) 
Fl = 0 8 >6. 6<64 (80) a 
Fl = 7 sgn (z-ax) 
if 6,< @< @ 
and 


(az + my) + |mz + ay| Ve? (a% +m) pale of 


2 


x > 
, (a + m’) 
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= and | Fl = 0 








= x> Getm)+ may p? (a? +m’) -1 


(a” +m’) 


i il 


a 


= and where from (71) 


Ly 


mM 





a. a _ _ 


(i - pa’) tan @, * p’am bs 8 (a’ ag mn’) “1 


ULM 


wT 


For 


MI 


= Oop of +2!) 


em 


I 


. | 
2 wap | cent 
1 


from (53), (54), and (56) 


Hi 


be x-Biatems) +7 L-ptaa Tern — 
a mae i- (o's? 2?o 
2 \-fiotem*) x2 Cat A 'latema)- Ue la Syant)] {x* a Z 


ee 20% F2 «- V x*- wie | : oe ees 


j~ Zoatyant) = o 
X~ Pi (adems) 7 


| 





Hy 


| 





ee TE! 





[Aiatem) -i){ x. plates?) 


ee ere | X- A'la.d+nms) _ 


(81) 





| 


uD ms 


PPE T TTTE RR | 


NTT 


: 


and from (78) 


2 


F2 = 0 x2 = g% (y2+22) 
= 0 x2 = g2 (y2+z22) 0>67, 0>0; 
= 1 x2 « B¢ (y%+z2) 6,< 6 < @ 
or outside envelope of Mach cones. 
For 


x? > pt (y*#2*) 


From (56) and (58) 


F3 = tan 
and 
lim F3 5 
(mz-ay) + 0 7 » ; : 
and for 


x" os p* (y*+2 ) 
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“< oO 
A A 


1 - p% (a2+m4) > 0 


RR 


1 - B- (a24+m2) < 0 


1 (mz-ay) 7 - p* (y*+2°) 
2 (z-ax) - y (mx-y) 


y > me 
MZ > ay 
MZ < ay 


(82) 


(83) 


C 


rar] 
ie ~~ 


from (76) 


if 
Oo 
© 
A 

D 
a 
V 
D 


F3 


= 1 mz > ay a, <O< a, 


= -7 mZ < ay 6,< 9 < 4, 


= 0 outside envelope of Mach cones. (84) 


Value of g on the envelope of mach cones. For a supersonic leading edge 
[i - ge (a¢+m2) < 0], in the region inside the envelope of mach cones from the 
leading edge and outside the mach cone from the origin, we have for a surface 
distribution of sources [from (29), (79), (80), and (81)] 


4 


i ; ses ie : 
\je? ar Ga peace pea’) a+ fe ge = (z- ax) sgn (z - ax) B? Go ans Lee oe pea’) - m (x - fl az) 


= -x|m+a sgn (z- ax) B° (a? + my a + z sgn (z - ax) \/ a (a” + n°) a ay a - #7a%) + Bane 

~ xX (a? + a) {1 - g’a’) + {rm - a sgn (z - ax) \B° fa® + nm?) ~ 1} [z sgn (2 - ax) \/ 8 fa on") -Ley(l- g’a’) + ganz} 

oes a ie ie ae ea Na Se 
m- aspgn (z - ax} \/ a" (a + m”) - 1 


+ 
(1 - pra’) - X (a? + my + (my + az) + sgn (2 - ax) (mz - ay) \ 8° ta’ + a”) - 1| 
MT 
m- a sgn (z - ax) 8? (a? + n°) - } 
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On the lines 


xX = 


(my + az) + |mz- ay| \/ p° (a +m‘) - 1 


(a° +m’) 
sgn (z - ax) (mz - ay) = |mz - ay| [See figure in Woodward] 


Therefore on the envelope of mach cones from the leading edge 


In this same region for a constant pressure surface [from (71)] and 


7% 
1 AP é 2 = 
9, fe (22) | (ai+m) 7 


= - a (ay - mz) \/ p? (a? +m?) - 14m (ay - mz) sen (mz - ay) + (a° + m’) (z - ax) sgn (z - ax) 
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C 


and in the region where 


: 
F 
F 
f 
a 


Sen {2 - ax) = sgn (mz - ay) 


+almz- ay] \ Ba’ + m’) - 14m (ay - mz) + (a% +m’) (z - ax) 


sen (z - ax) 


= spn (2 - ax} aj- x fa 4 nm) + (my + az) + |[mz- ay| \/ Be (a° . m) ae ! 
(my + az) + |mz - ay| Ve? fac + mn”) - ] 


; eee _ =* “es _.. = @Q on the lines x = 0 a ee 
ot Sr | aa 7 | 4 ao+m 








| “Therefore, fs ~= 0 on the envelope of Mach cones from the leading edge. 


Verification of the imposed boundary conditions. Now (8) can be verified 


for the case of a surface distribution of sources. From (29) and (79), on 
. ee ax, since the cosh"! terms are continuous 
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For (8b) using (29) and (79) and the results of 2.6 3 ) 
eo! 


2—- 
u ® - MAES [ra + mF2| O<y<mxk 2><ax 
£2 
m(l1-pa) 


2- 
- AE |-ma + nF? | O<y<m 2<ax 
n(l-Bpa) 


= i | 

- +6 au | + p°anF2| O<y<mx 2Z>ax 
22 

n(l- pa) 


w+ Beau 


7) | + pam? | O<y<mx Z<ax 
n(1 -B a) 


_ Therefore 


— 2 | = oe 
(w- w') + p’a(u -u') = 7s ri[l - pa") = 2[w + ret = const 


ml-pa) 
which agrees with (8b) 
For | 
y < 0 
or 
y > mx 
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where 
f = m(z - ax) \x" 7 pty + * 
Br, p*alay -mz)j + (z ; ax)” 
and since 
f = Qons 
when 
Z = ax 


canVx2 (1 - p02) - By 


Ou _ a Sa 
8x yy - mx) (1 - Aa’) 


and likewise on 


Z = ax 
du 
——— = f 
ay 


\/ 2 
—aVeq - 6222) - 2 


eu 
Bz yyy - ml - B’a*) 
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Therefore 


ou Ou aaa” 
Q@x’ dy 82 


are continuous on S (except possibly at the edges) and therefore 


_o- Peat - 2 - F ee: - 


and (30b) is verified. 


u') = 0 


Surface velocities on a constant pressure surface. - Flow normal to 





surface U 
n 

ui = a (-au + w) ™ 
1+ aé 

t 

usié : (au’ - w') 

ie 2 
l+t+a | 


th [4 


on 


(Zz = ax) utu' = 0 
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a 


|-a(u + u’) + (w+ w') 


« 


from (51) and (79) > 





- AP\}_ an = Les, X + Px *(1- prot) ~ Brat 
(a“ + ni) é 


= + w') ee ey 





x = x? ( )~ Ba?) - p'a* 


aya Na oe a Dy Sn 
ACB o)-Bimy - VL parses] 





9 (a° Sun) 


-m tan 1 (mz - ay) V2 - 2a’) - py? 
RE ie -y (mx - y) 


eo ee eee 
= 5 l+a (u, u) 


1 | 


“ee , n n 
, ; ) 
\/1 +a | 


net = ee —_ 
7 ewe 


i. ee , | (86) 
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Velocity discontinuity in the leading edge wake. - 


Constant Pressure Surface 2 





O< y < Mx 


= maha at [ =) + ae] (a + ra tie 
2 


4may /1 + | 


= —~4n__(¥) Source Distribution! [see (86) ] 
ee 
 4aVao +m 
: 2 a 2 2.-1 
Vp > Vp -m' =| +m) 0 


2 AP alan 
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1 ace Sac 


ibe i 





| ee fea , 
Alternate sign choice of tan ~ denominator. - 


may - mz) 1 (mz - ay) (mx - 9) - 8° ECan - mz)* + m°(n- y)*) 


(a i m2) a(z - ax) - m(m - y)Jn + mfy(m - y) - z(z - ax)] 


[see (45) and note sign change of denominator and since _ 


[a(z - ax) - minx - y)Jng + mly(mx - y) - 2(2 - ax)] < 0 


The above term when evaluated at 


is | . 


m(ay - mz) msgn(mz - ay) 
2 2 
(a' +m) 


Evaluated at: 7 = 0 the term becomes 


m(ay - mz) pe (mz - ay) V ee 5 rig + a 


ah Seer [y(mx - y) - z(z - ax 
| ca 
2 2 
“ » yim - y)ia Fm 
m= ay  —-y(m- y) - 2(2- ax) = YW Bd 
m 
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which means that on mz = ay and n = 0 


C} 
1 msgn(mz - ay) Vo ee 
tan 7 
0 0<y<m 
ae Ate 0 y<0 y>mx 
tan = 
0. mspn(mz - ay) O<y<mx 
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Subappendix A - Woodward's Subsonic Equations 


If Q(x, y, z) satisfies 


ae Se bare 
OX oy 82 





Then we can write the solution for Q(x, y, z) as 


soer-fffow n,¢) oo a | = 
‘S a4 | 


- where 


Q(x, y, z) = 0 B° = 1-M 


dS 


x = \/ (x -£)*+ pe (y -n)* + ae ple 


and 


Ge, Nw dir... iets 28 
av | 8 eRe ttn Me BE | 


mo = (Mp. y Ny) is the unit normal on S 


We will assume S is the surface { = aéand 
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(Al) 





surface distribution of sources. - Let 2 = @ and as for supersonic 
flow assume that on S V1 + a2 P= (- 82a, 0, 1) 


a) =! 
b) oot Soe Z (W- 8% a) = const. Wew-wi 
| lt+a 


the primed quantities refer to z< axor {< aé 


Therefore since on S £= aé we get 


. Wis B au 
@(x, y, 2) a 





(AZ) 
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This is the eae as (9) for the supersonic cuse except for oe linits, 
a factor of 1/2, and the fact that Bis replaced by -g. ‘Therefore the 
integral over & ave be performed by using the equation below (9) for the 


supersonic case, and replacing 6 by - B 






ee ae | (a3) 
oe ers ar ae ee + B ee Di + 6 Cae Se 
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e(x, y, 2) = + Ga-1) + Bia (mn 8-07) ay 





41 = (mmx) = fia (m2-07) 








= |. | | (A4) spice 
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J (iegha2){ (anh -71) 4 Palas Bnb-a9) 4 bye) +a (ma-a7) ie 
| 22 
| =e Vee? a 


SOB) (Laney s ply -gya fh CmB-og)”  — nk) - Bales -07) 


we J . - ran : ; . A 
5. phere Me xox, ond we z-2, = Z-ax, 

















aie The wean party is “the | same as the first ‘met ‘that > x, 2z ee m are 
| ; “replaced by x-x,, 2-2, and m These integrals are the same as for supersonic 
flow, given by equation (21), ‘if B* is replaced by - B and a factor of 1/2 is 
as 2 “added. | | 
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When the limit 7 = Y, is calculated it can be seen 
that it will be the same as when yn = O if x,y and z are replaced by 
X-X2 y-y¥2 and 2 - 2. 


Therefore . 


w- pau : : ; V2 ie ae 
6, ey, 2m) e wf au —2- ox gt mz - ax) Vx" + p°(y" + 2°) 


Zs 22 
yL GY - mx) * p -a(ay - mz) 14 (2 - ax)? 


‘ “ + p*a*) - m(x + § 242) - log a [ Te 2a" i n)] [x * 6 ty? * 2)| . [x * B 2(az . ry) | 
e “at VA ee “ea! : ") V [a + g(a? +n?) | x ¢ p?¢y’ + 2°) | - [x + p? (a + ay) | 





ye +p tet [x2 +p 2y7 + 24)) + + 62a) 


: y | 
V 22 = 
ei a+ 2%) [x2 +n 2y? 6 24) - xe 0 2a) | (AS) 
- Now we can write the result of (A4) as tes 


fo 
po | 1 | 
{ ' i 


’s (x, y, z) = ¢,&, Ys 2, m) - o,& 2 Xo» > is Yq» z- Zo» m) ee ete : 

| 
| —— 2 
| 9 & a Xq0 y i Yz> zZ - or ay m, ) +o & 7 X43 y 7 Yq? z- Za? | 
: = 
{This means that (AS) may be interpreted as the velocity potential for an infi- 


‘nite panel with leading edge slope y = mx, and at angle of attack a. 
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Constant pressure surfice, - In (AL) we will usc u(x, y, Z) 
= Q(x, y, 2) and 


’ 
gu, ou 


a) oy a 0 on § 


b) u-u’ = const = Au on > 


and where S is the same as previously defined. 


refore (Al) becomes 
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(y -7) + (z- até) 


————————————SESEST didn (AG) 
Vix -@)° + p*ty - a8 p*(- ae)’ 








—S Se ee ets 
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Now we cah write, using AP/2q a -Au 














| = 
ad 
x r 
@{x, y, 2) 7) u(x’, y, 2) dx 
Y2 ae 
- ps- (x -& )(z- ag) : 
: aa | 2 aac haere rte 
. i - &- E+ Bo - 1) + Az - ag) 
n 
(A7) 
+x 
% pm 
_. AP z- ag) a 
- a tdn 
ae f J (y-n)?+(@- ag) 
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aes 
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_. The first term, except for the limits of integration and an additional factor 


of 1/2, is the same as the result for supersonic flow, which is given by the | 
first equation of the section covering the evaluation of the integral over & ' 


for the constant pressure surface, with - -B in place of ®. Therefore, as in 


' the case of a constant distribution of sources, we can use the supersonic 


result for subsonic flow if we substitute. =p for s. 


The second integral is new and does not occur for stmersonic flow. 


This integral will now be evaluated. 
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Garde -S= log v-n)+@- at)" 
7 ny +(@-ae) , 
i | (A8) 
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Therefore © : 
Tox 
Y2 my 3 Y2 
(y- 1) + @- ag; : 

o 
uu 
no 

‘ = log so “7 i + (, : 25 : “| dn 
ae a, 
Pees 4 —_ 
note that: 


| 


‘mz -an = (m - ay) - at n- y) oe ax) + a(m - y) - af n- y) (AI) 


; 2 
i . . 


ea : _ 

| Fen eae ee | —_ 

| log | Sony) * ten d : 
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| 2 2 2 (A10) 

. - Boy > 2 = | aie See ee Se )_ + (a 9 - mz) | zl 

|. a +8 | a 4 
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| Dae - ey) oa? m(mz - ay) : 

| at en (n- y)(a’ om’) - afm - ay) 
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F a eee cae momen eee = ee ee | 

a casa ciimigaias. arabe ease ——— eee Fa 

241 | 








When this is evaluated at 1 = O we get from (A7), (A%) , and (Al10) 


“ AP -m(az_+ my} a 2 by 
¥ 09 Ms y, Z, m) Srq a = 9 7 E log(y +z) 1| 
80 (a +m) 


. . m@nz - ay) 1 (mz - ay)m 


(a 


-tan- 
+m) .  ra(mz - ay) - (@ +m)y 
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Since all of the terms in (A10) may he written using only the variables, 
(m - y), (z - ax), m, and (1 - y), the additional term which must be added 
for subsonic flow is, 
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If we integrate this over S we obtain (A'7). Therefore (A7) corresponds 
to a volume of doublets which means there wil] be a surface of sinks on one 
side and sources on the other. 
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Subappendices E and F consider the case where g(t,7, 7') 1s considered 
to be a function of ¢, as in (3), but the expression (26a) is still 
infinite. 


If tan @ is not continuous at some 7 = 7' then, from (23) and (18) 
or (26), it can be seen that the integral of 
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will give a logarithmic infinity. If tan 9 is continuous the first term of 
(23) will give an odd function of (7'-1) which will give a finite value in 
the Cauchy principle value sense. Therefore kinks or cranks in the planform 
cannot be permitted because discontinuities in the angles of the constant per- 
cent chord lines will result. In fact subappendix G of this appendix shows 
that a logarithmic infinity occurs if the leading edge kink is considered as 
the limit of a hyperbola whose radius of curvature goes to zero. If the term 
in (18) involving tan¢ is neglected for small (#'-7), this has the effect 

of rounding the constant percent chord lines. 


ww 


Now let 7' = cos 6' and 


f(m') = > £5, (@) - (27) 


m=1 
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where 





M 6' = @ 
_ 9 1 mM 
YY. : : t = 
Sf 6') ved Sin p 0. sinp 6 mei .- . 
(-1) sin 6@_ sin(M+l1) @' 
B=] 1 m a £8 
M+] cos @6' .- cos 0 m 
_ mt ; 
O eal discrete points 


Now the first part of (18) may be performed in closed form 








1 an ue M 
i. mee: ae oe 
on oe i > Fam Pym 25) 
-l m=1 
where 
M+1 1 
4 sin6, | on om 
Oi a +7 in 6 
De aoe y sin _ 3 
M+] ? Vv On 


[cos 0, - cos “i - (29) 
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Referring to (18), (19}, (26), and (27) we define 





oo Li i bg: 


1,46 


- kK (é > 1s 1) 
aoe 7+, 4 p #0 (30) 


- q 
ny, =6 Cm, "? 


= -1 
f vt 
KC ay? ny 0 ) Sa cos } | 











and letting 


1 +n, +6 1-7 


EE a a Oe a a 
: . 2 1 


i} K E ono” co) 3 cos 7 ca) 
Peseitie ths Tye i aaeMs : sin @ dé 





Mh x aE »n i» 14 ( 0) S_ feos”? ay a) 


tel | my - 14 (oT 





sin so 
ae 


; ; : or 


where 


Wiehe YL sesiay 





@ oe AE. £12, cee M 


| 
pe 
ous 
+ 
+ 


| 


ow 


| ' ' 
| i ikieddqe WIMP. [biser ie al 
: ' i oa | ria HO Peek Ce Te 
ii Wiel to: sp ha 1 
i Bil: Lit ce gy Ladle tea by 
|: Blt H i Vl esate i Ia. | | 
4 Pe | PIS if $file vei] i H ia | 
: 1 ig Bil ith diits il dle Hd ab dT Ube effet} ad dh NTA title 
} iH 5 wid Bde 1a BIA Al! ate q 1 ‘| ik ‘ 
H be alta) LATA La WEE od aia! Bel Raa 4 MDH Is ea Bal Lae nly at 
re Liab Ue Te UTE Fee, sede ut" 1 Hy eae | Rs A yey ea a |i eye 
ot ‘ bord nye . 








— Wad ae ie | ae 
npym 2 2 
-j (Ny - n') 
and letting 
A Mg 8 1+7,, -6§ 
ae n(6 ) = oo eee 
uJ oC * 
1-"ny -6 K [ Epemyem2(8)] S,Ce) sing 
C = ——_—_—_ is saa ee eerie Ree P 
npvm An 
17 Vale 
} [ ™- 726 )] 
M a e 
tatvce sg Tlerter tet a pad. 
~ 4M, + 1) is SET i” Mel (32) 
1=] [ 2 ?| 


Therefore referring to (18), (19), (26), (28), (29), (30), (31), (32) 


N M 
a(t p’ ony) 7 S > to == Son apy . B ppm 7 Cocina 


which is a set of linear equations to be solved for f 


Pe SO sects: aN 
Mx (N+1) 
v = 2 eae @#eosa 9 M 
Hippy = Ae Cb Nps Ny) 


all of the H's and K's are computed numerically 
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Subappendix A 








Consider 
a 
= lim J (nt-n) FC 7') ae 
= Zz 0 ees 
= -1 ( n Me n ) + 2 
5 7 fn) (n'-n)[FC 9’) - F¢n)| 
2 = lim F(7 ) a n! + pag eg dy’ 
s 2-0 4 (7 -4) + 7 fas “1 (nt-nj) +: 
, a... 
3 . a 2. 2) =m) [FCW ) = FCn) | 
= lin 5 FC") log |(n'-n) +253 4 a A Ge n!' 
2 —» () 2) (n'-7) 
1 
Te Peat ae ee 
1- ; (n'-n)[FCn') -F(n) ] 
= F(n) < log _+ log = + lim ——— ee n! 
Z =e 0 i ( nt-n) 


If F( n'} is differentiable at n' =n then F(m') - F(n) = O( n'-n ) and 
we can write 


Dae 7) [FC 1') - F(n) | 


Lim —— ot. n' 
eee (n'- 7) 
i n-€ 
ean | f FIRM age fray - FEM) gy 
E==() nee (n'- 4) “1 (n'- 1) 
(AL) 
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Therefore 
1 
j F F | 
Lin _(n'-n)F(m') mF Cr) AP: Zi 
z= 0 (mn)? 
-] 
1 n-e 
a we F( 7') F(7') ; 
_— Cnt nyo * Cn-ny 
N+e 1 
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Subappendix B 


He (p,q) = Ho (P.s G) - Vi (p,u) «de 


es , 

| On DP, q) aut (nm, q) 
Hegel os 

uv O4 


efter 











da 


ti 
— 
3 
a 
ed 
+ 
Sa Se 


(B1) 


where C is a contour integral in the p, ; nlane fron (po Oto ip,q) 


l : 
H (p,q) = he ee eee ae 
nh : n —— 


| 
Hei if i Fae ki 


Ii lik ; 
i Ltn thie sii : 
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|.‘ Hil 
1240) ara tL 


Hl 
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aH, (py q) 2 
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rere t atqtglafyp pete pelle Pe ah yp Pet Madea GET Cpe pede git a crptee boy S| a lL | | CO UO | pe flap lltte pays PF {1 een ‘ 


iti] it 


t 
Le 
RO 
| | 
oS | 
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J 


if 


Hits el fet cod Wl 
Mery dale fesfep i} befe-delfpe 


at 


Pp #0, 1 


we can write 


5 


h¢ E') = h. (P) - h’, (pP)(p- é&') + = ee (p)(P - Bry? em ae ae 8 Cae 
| (B2) 


sl Wail na i 


I ll 
i lt | | 


ela LAU) | fd vt 


sk (loli 15) da 
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[ib 
(USHA PSIG ALLS ol 
GL 


lll | 





| 
mi 


f 
mil 


and since 


x 1 x” xr a° 
a eae —z ees + a7 log (x+r) 
r r 
and 
2 
oe ep ie 
| 3 Yr 
r 
where r # x" + a* 
g'=1] 
an are ee eee 2 
3 H,(P.a) = h_(P) eer + 0(q ) 
VO= 49 4 
£'=0 
or 
bog! Be HP.a) = 2h, (P) P # 0,1 (B3) 
1 
ao 7 = - 7) (p- Ex) ' 
a Hy fP q ho E') - | awe 5/2 d & 
[we 2)? + a2] 
O ; 


now write (Pp #¢ 0,1) h_( g') = h_(P) - he (P)(P - &') + rf E') (p- E14 
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Then as q — 0 


(P - E') <= See, Cee oe = 
ee a d & O(q) 


| @- gy? + i] | (- E14 + q’ - 
oO | 0 


and . 





ee cae eee y r CE) 8° 

tt nce, gh a te ee d &' = Q(q) 
a ee ee Z fa Spe 
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(p- gE‘) + log fcr + \/(p- gry? + a | 


2 2 
(P- &'}) +#4q 0 


eee a-) : ia, Maat wie tae eras 
eee Seer 2q9(- —- + log ti Pp) + (=P) t:q | log p+ Vip +q 
= sae as faa-»)° i a Ip" 4 a 




































































2: 2 
> swe tia bans Pp +#q 


2 
= aereEe ae 


q 


= 0(q) - 2q log {q| 





277 


if eae 
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Therefore as q -- 0 


H, (Pa) = - 2h* (P) qin |ql + 0(q) (B4) 


3 
|} 


and since h,(P) and h',(P) are slowly varying for p # 0,1 referring to (B1) 
and using (B3) and (B4) for (P-Pp) << 1 and q <1 


P 
SH (pq) d = f 2 fi pq) a= 2h ()Op-p) 
apn ‘F? ap “ns? no’ ™ ‘oO 
PS 
and 
q 
3 2 
— = : - = 
f 5q fg (P29) da 2h' (P,) q in q dq hn’, @ja Ing 


or to lowest order in P=e. and q (for Pp. # 0,1) 


~ _ ws n | 2 ; 2 
H (pq) = alPy?0) 2h, (P,P eae a(P.) q ing (BS) 
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Subappendix C 


LeFrpeiile tb dpe ds 


(P - E') . d f! 


Ha) = | h_( Et) | 1 + oo: 
(pP- —&') +q 


Leading Edge Let 


Bb AEAGLE TR ae 


een 


or 


S2 

il 

ct : 

g 
als 


and introduce a change of variables 


gl s sin’¢ 


dé'- = 2 sino coso do 


i | q 


Hil 


I) iy 

uth i [Heil 
yD ele Lae bpd! 
ern MARTH 
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re 


1 eae 
yt 
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NORE EEC GL Lee 
ML 


i 


| 


ee 7 1 _ ae 3 , 4 cos éng + cos 2(n+ljo 
h, (sino) h, E (1 - cos 20) | ho Qo) == Siar oo . 
2 cos ([2n+1] - 1) o + cos ( (2n+1] + l}o _2 cos (2ntl)¢ 
™ 2 sino cose 7 sing 
Pear 6° es sin @ + sin’ a) 
0 ee ET een anc 
\ (P- gy? + q? Vein’ o + 2 6° sin ry sin’ o + 5° 
and 
H (P,q) = Hf 6, ?) 
. 
4 - (8° sin $ + sin’ o) 
sae 1 + . cos (2ntl)ccos ad. 
. 4 2 <a ° 2 4 ara = 
te \jsin'o+ 26 sing sin a+6 rd 
us 
afi 2 _ 
eo & - 2 sin @ cos [(2n+l) ¢] cosa da d 
pe \f.4 2 2 rae 
A sino +26 sin @sin“o + 6 
* 
4 : 26 | + sind sin’o + 32 | (87 sin @ + sin’ a7) cos [(2n+l) o] cosa do 
4 
is 3/2 
(0) | sin! o+?2 6 sin @ sin o + 64] 
us 
ee ee 2 
8 8° cos @ cos [(2n+l1) a] cosa sin’ a oe 
" 2 
- isin’ + 2 8° sin ge + 34 


as 6 — 0 the only portion of the integrand which is important is in the 
region of o = 0 
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Therefore let x = sino dx = cosa do and expand 


| 3 
cos [(2n+1) a] = cos = x= =| 


| 2 
21-4 cane)? x? 21-2 1+ tn + An’) 004) 


: 2 
Therefore to first order in x 


| 1 2 21 2 
|1- (1 + 4n + 4n*)| x 
oor 5.3) = 8 oso 8 : — dx 
86 n i - 4 22 nm 3/ 
0 x + 26x sin 3 + 6°| 
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. 8 cos’ é x7 dx ee JU An 4n°} y dx 
, 4 oe 3/2 2 1 a 3/2 
O |x + 2x sing +1 olx + 2x sing+1 
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i 
fe) 
oo 
_ x? dx x? dx a dx 
| 3/2 3/2 °% 4 
0 ie «Be sin$ +1 [xt + 2x sin 3 +1 : 
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5 
oO 
o 2 
eee <* x dx 3 


Jo [x4 + ax? sin $41 
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Near 6= 0 we can Say 








00 Zz 2 
7. fy (1 + 4n + 4n*) 
da ~.  § cos ¢ 2 3 
96 1,65 %) = TT 3/2 dx + 0(6 ) 
. 4 2 a rw) 
O - + 2x sin ¢ + 1| 
integrating from 0 to 6 .we get 
on wl 8 cos“ ¢ ~ 1 2 3 . 
HL (6,9) = H (0,9) + “po 1(9)6-— (1 + 4n+ dn) 1,(9) 6 (C1) 
or 
H (p,q) = H,(0,0) 
2~ 1/4 — 3/4 
+5 cos o 11508) (p? « q?) -2 (1+ dn + an‘) 1,(8) (Pp? + 4?) 
(C2) 
where 
ry = eg CP) 
q 
and where 
1.(%) = ee 2: 
‘1 r , _ 3/2 cos@ dé¢ 4 > > 
E + 2x sin 5+ 1) ‘ x + 2x sing +1 
0 
2 : : (C3) 
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4 
ms x dx 


By + 2x” sin g+1 





ms ~ Ca + i sin ¢) 
= - sing 109) + oie eG aap (C4) 
' x + 2x sing+l1 


3 
: (x” + x sing) dx 
dv y" 





x! + 2x? sin O41 















































Panta aM came eee gee _ : 1 
cranes ae 2 
a | - VY © S=S Se 
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. z ~ dx 
= 18) = - sind (8) +5 J ———S 5) 
== \/ 4 De is, HD 2 
= , 0 x + 2x sing@+l 
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Subappendix D 


TI - ' + p ate ' 
Hs q) = h (é") 1 +— | dt 


V(p - by + af 


O | 
Trailing Edge, let 
p= os 6° dai 6 
q = 67 cos @ 


or 


6° Va -p)? +q° 


Be. faq 2 ED 
q 
and diange: variables in (D1) 
glo o= sin’ o 
dé' = 2 sine coso do 


2 cos (én +1)¢ 
ce 


7 2 
un (sino) ' sing 


(see leading edge expansion) 
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neat 
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=4 fo 
4 ——— 
ti J=a 
TTT, 


- §& sin6 
Vics? o a sin 9)? +6 cos 


| cos (2n + l)o cose do 
4 2, 
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Q 2 
—_ 4 ] + 
~ OF 
O 


sin o - 5° sin8 
Vsin’ @ +26 sino sin@ + i) 


le can + 1) i - a) sing do 
4 


Now 


[can + 1) e - o)| 


2 


| cos | cz + 1) | + sin cen + 1) | sin c + 1) | 
= sin en + ]) 2 


sin | cn + 1) | = ie sin [(2n + 1) o | 
The differentiation of the integrand with respect to 6 may be performed easily 
if it is compared with the leading edge case. 


cos zn + 1) = 


; en 


86 


“" 
n 85> “os 6 
Le * Coe 


hea [ (2n +1)o] sino do 
Oo 


2 43/2 
[sin‘o - 26 sin@ sing + i | 
Changing variables of integration once again 


xXx = sing 


cosa dag 


| sin [ (an + 1) o | 


(1 + 2n)x + O(x) 
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and 








] 
Q n 86° esa : Ae dx \M 
—H (6,0) = (-1) (1+ 2n) —— ee 
06 n 7 A 2 2 r 5/2 
O c - 26 xX Sing + 6 | 
n 86° e656 as 4 
= (-1)° (1 +. 2n) —— oe ce a. a3? ee CO Oe) 
Oo Ix - 2x siné an 
De cst 
= (27 (2 + any) 82-282 Coe) 
(13) 
where 
foe) 
4 
Ty (@) = x _x___ | (D4) 
o 6U6[x) + 2x" sind + 1] | A 
and integrating from 0 to 6 we get 
2 3/4 
~ _ oo - iam 8 cos 6 Iz (6+ 7) ae Z 
HO, q) = HG: q) + (-1)° (1 + 2n) 30 hee. p) 
(D5) 
where 
@ = tan EP. ) (D6) 
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3 J . Subappendix E 


Suppose that k(x', y') may be written (as before) 
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a k [x'(g', n'), y(n] = 2 (1) Rte) 
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fo) = h l20- cs ¥) 
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Pad i stb ii 


Then with 


and from (3), (2) and 


where 
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N ay 4 
not) = > fn) HCE 07, 8) 
n=O 





288 


and 
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x (m) + &c(n) - x(n") | | 
: | he) {1+ eee 
o _— Ec(n} - x(n") | n= H Ea 
Vo ea ee |, ee 


We can also define 
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: (p- ect) +q +8 
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Leading Edge, P eo 0, ¢ = oir = leading edge sweep 
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Deh ae ete jee teeter = f=) 
pe + q ecm) ( ) ? i, 
2 
2 
l 7 b (n' - n)* + ¢? cos“¢ 
cos’ ¢ c(7n) 
rap n> 7 


Therefore at the leading edge since @ = tan*(=2) = 
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_ Ne" ae 


p + st + ¢ \i(n' - + ¢? Bane % Fi + eee ® 


sin 


where 


letting s = 7' - 1 we must evaluate an expression of the form 


f-0 


} 
ms : in’ @ sin? y } IT (F ) t «? + r2 operé ' Pir 5 20 
lim 2 f/o= »¥ , ds $= 3 
S 


2 2A A 2 
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where 


sin rf ae (no ¢ dependence!) 
eae + cosa 
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= 8 2 
i, x dx 
0 x + 2x sin @sin J + 1 


The second integral goes to zero like /g and the first is independent 
of £. Therefore 
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E~o af Jt 





292 


ban al 


® 


Subappendix F 
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wr (p - §') 

H (P,q,%) = | h (€') 1+ —__ 

.: ns ee ee: 

O (p- &')' +q +8 
Leading edge, let 
Do te 
P = -6 sin @ sin 


q = 6° cos A sin 
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5° _ n° A a x re 
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introduce a change of variables in the integral 


Et os ano 


dt' = 2 sing cosa do 
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dé 


h, (sin? a) = ho E (1 - cos 20) = h, (20) = 


sin 20 


S Cos [ (2n +1) - lle + COS [ (2n + 1) + 1Jo 
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2 sing cos o 


T sin o 


and 
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2 cos 2ng + cos 2(n + l)¢ 
TT 
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2 


as 60 the only portion of the integrand which is important is in the region 


of gd = 0. 


Therefore let x = sin o, dx = coso do and expand 
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Therefore to first order in x 


2 
J F Beane aya? a 


83° [1 - sin’ 3 sin’ | 


2 
a H (8,6, ¥) * 
3 n + Os ; ae = 
. O [x +26 x sin@sin % +6 


1 
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Integrating from 0 to dé we get 
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. 2 2 
Pp +q 
2, 2 2. 2 
a aa +¢ p*+qs, 
A2 
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Subappendix G 


S 


* re ; 2 
H (p,q) =H, (p,.0) + 2h, (p,) (p- PQ) - hy fe) a intl 


Hyperbolic leading and trailing edges, constant chord. 


= 
° 


[x,Cn"9] ae Paes lat? + 0 an 6 | 








xX (n") = c tan é@ nt? i p° tan’ @ 
nar | 
1 dX, 
' Saree 
n'>> ] ¢ dy tan @ 
A eee, a : 
“ar Curvature at n' = 0 = 5 p 
dy' 
n'=0 
From (21) at 7' = 0 
e ' 
ae X (0) X Cn") 
o C 
or 
pe Re m= = tang NEE +p. tan’@ - p tan >| 
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Therefore from (15) and (20), for small 7n' 
Vee 2 2 
HO (& ,0,7y') = Hh ( £,0,0) - 2h, (&) tang n' +p tan” @ - alee 
t 12 t 
+h’, (6) nt In| n" 


Therefore referring to (17) and (18) we must evaluate 


6 
5 | ee | 
lim “a> ds 
re) 
f—-0 : s* 4 ré 
fe) 
6 
h| ee + = J 
ds 
a s* 4 r? 
6 
\f 2 2 2 2 
=f 108 | ee at], Var e? rs ee Ga 
4 4 'd 
, r Z 2 
5S +a 
O 
- {sale jaa] Ve wt Pee -2t 18 tog a 
c £ 
tVe2s at 67+ a 
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‘ 


and 


¢ 





lim lim 
f—-(0 Soe ag 


a--0 
0 


ptan 6, there is a logarithmic infinity as the radius of curva- 


and since a = 
ture goes to zero and the leading edge becomes kinked. 


J) 
299 


APPENDIX F 


POTENTIAL FORM DRAG 


The section load can be obtained by means of the Blasius theorem as 


follows; 
Z 
a dw 
Fy 1 Fy pif ( dz 
C 
where 
dw é 2 Z 2 2 
(Sx) su - Zuvi + (iv) =u -v ~- 2uvi 
Therefore; 
ul es eee 2 
i if," ze | lauy + 3 (u . yy] [dx + idy] 
C.. 
then | 


pe-te [ | 2uvdy + (u2 - v2) ax | 
C 
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(1) 


(2) 


3) 


(4) 


(5) 


Since all of the singularities are on the chordal plane, for the panels, 
equations (4) and (5) reduce to; 


J 


N,. 
J 
- ae (Yu, mie. Fes oe ae (6) 
a) he ; 
1=1 
N. 
1 | 
= Bo iol en Pes yee 
'y - nd. (%, ; %, | (4, : “il ) 2G (7) 
1=1 7 


where the subscripts u and L indicate upper and lower surfaces, 
respectively. | 


For the two-dimensional lifting case; 


wi 
u, = 3 Y, + V.. COS @ (8) 
el 1 : 
eine eae 
u, hau ae COS a | (9) 
i 
a 
T 
a oe ee K : 
‘a, = L, ee = > me % +V_sin (10) 
K=] 
Therefore, 
N; N, Ny i 
F =-p » Yr. V since oe Eo (11) 
xX io a 2a Xs - Xy 
i=l i=l Kel 
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Since 


. 4 
r..T 
p ivk _ 
On > as (12) 
i=l (1 
ro 
i 
Po=- pV, sina 2 TY, =- Lsina (13) 
1=] 


where L is the lift and I = YAx is the local vortex strength. 
Also; from equations (7), (8), (9), and (10) 


N. 
1 
ie p > r, ve cosa = L cosa | (14) 


i=l 


which demonstrates that the discrete vortex lattice always gives the correct 
chord force F, and normal force Fy, provided the correct circulations is 
obtained, | 


Similarly, for the two-dimensional thickness case; 


N. 
: >> 
1 K 
UL. AT. 27 yy X. -X i ve (15) 
1 1 Kel 5 ae 
vv (=) | (16) 
u. eo\dx/. 
i i 
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i oe 


(Se 
My 7 7 Ye Vax], a7) 


and 
(= 
2 =2V, ey ae (18) 
There fore 5 
N Ne. N. 
2 1 . 1 
p Bi My 7 - 
F. = = -— ———— - pV zZ. = 0 (19) 
X 27 Xs * Xy oe) 1 
i=l 86kFl i=l 
since 
a N. 
1 
=. = O 
i 
i=l 


if the airfoil is closed. 


Also; F, = 0 from equation (7), (15), (16), and (17). This demonstrates that 
the discrete source lattice also gives the correct chord force F, and normal 


force Fy. 


The above equations can be generalized to compute the section potential © 
form drag on a finite wing due to lift and thickness by evaluating the compo- 
nent of force in the free stream direction and by using the three dimensional 


influence equations. 


0 


33 


The section potential form drag due to lift is computed as follows; 





na 
4 =-P > a (29) 
7 i=l 


where i is summed over the section of the panel and wj is the total veloc- 
ity normal to the panel chordal surface at the quarter chord of the ith sub- 
panel. The section induced drag coefficient 





“a, C 
i 
CaVG. 
is then given by; 
C C N, 
c % r W 
i ZAR i i 
C SOS ( | (21) 
AVG. inl S = ed 


If wy is computed at the three-quarter chord of the jth subpanel , 
instead of the quarter chord of the ith subpanel as is done in equation (21), 
the section zero percent suction drag coefficient 


G. . c 
qr20 


CAVG. 


is obtained. The section leading edge thrust coefficient (Cy C)/(Cy, ) is 
equal to; | ; 
ee 
to ees ee es nea 
C C C 


AVG. AVG. AVG. 
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The section potential form drag due to thickness is computed by a similar 
procedure. 





L 1 
: Vx, 
or aes > ay ae ae 2 (= Je (23) 
1 % . 20 
1=1 i=] 


where i is summed over the section for both the quarter and three-quarter 
chord stations of each subpanel. £2 j and uj are the source strength and 
total velocity in the free stream direction, respectively, at either the 
quarter or three-quarter chord point of the subpanel. The section induced 
drag coefficient 





Ce C 
1 
Cave. 
is then given by; 
C C 2N. 2N 


Veo V., 
i=] 121 


p—5~—/z.\ /U: : ————__—_= Vi ; 
es eee wees S 1+ can) es ? = » A: waa, (8) on 


ie 6) 


For the special case where the chordal surfaces of the panels are planar 
and parallel to each other the integral of the section induced drag 





over the span can be shown to be identical with the value of induced drag as 
computed in the far field, provided all of the lifting elements (bound 
vortices) are parallel and the lateral widths of the horseshoe vortices are 
equal for the complete system. | 
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The total drag of a wing as computed in the near field is given by; 





N N | 
= gue oe ibe ot ie 
“D, b Ee ) GE ae 
j=l Kel Jk j 
where 
fr: 2 2 ee 
Ww _ (Wa) K Be oa ol ey 
() Py i ! 
: (Mn) 5) Gc) 
2 2. 2 
ie aes) Ge (T +P )y, 
2 ad 2 
OK Min) Be ry, = (¥5, ’y) 
X., + 7 
2 1 ie JK IY 
MgK * X +ty,)? +87 (v6 
| ( jK 7 jk "y 
X., - Ty 
=e Lt (26) 
jK “vy (x - ty, }? ae(y_- 
jK v jK es 
and 


Rig ™ Ry 
Y, : 
jK j 
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Qe 0 


(X., Y.) 
J J 


g 


Equation (26) can be oe 
die Tithe elere: 


stagger of 


Therefore ; 


is the influencing point 
is the noint being influenced 


is hait of we soanvuse lattice spacing 


‘s the taneraut of the vortex line sweep 


oe i parts; that due to the near field 
2 that due to the linit of integration at 





Ee 
|= 
, a . = 
gf 
ee 
~ 


— 
= 


. = E + £E 
v_) oe oO , 
iK jk jK 


(27) 


The contribution to equation (26) from the near field limits of integration or 
lifting element stagger is given by Ess 














ee 2 2 2. 2 
Bg hg ee By, BY, + MH. - 1 +B) wy 
E a a ci et Cae mss —— an - 
Say eens es = > 
jk . 2 =z ) a a ra Z (y . } 
Xe ° TY5y) J ae * Ty) tB° (Yi Wy (X TY) t( Xx Typ)" +BY (Yay > ¥, 


; ree 2 
nee P(e 


K+ Vy 
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That due to the limits of integration at infinity is given by E,, 3K 





1 1 
Be tee (29) 


The contribution to the total drag Cp;, given by equation (25), from ESiK is 
seen to be exactly zero’ for all planar wings and loadings provided both, 

T and yy are constant everywhere on the wing. This is due to the fact that 
there is no contribution to the drag from Es;,~ when Xjj = Yjj = 0 or when 
Xxx = Yen = 0 due to taking the Cauchy principal value. Also, when j # K 
the drag from ESjx is zero because the mutual interference drag due to the 
stagger is zero.” This is seen by interchanging the influencing point and the 
point being influenced and observing that ES sx = Esk; : 


Therefore; 





AR - 2 ) r 
ne 2 2 ro weak (Va), om) 


N N 
AR ua ae re: 
2° $l LOO ey) © 


© ae 
jel Kel Kk” 4\) Ihe. oY 


This equation is identical to that obtained from the standard Trefftz 
plane analysis. From reference (47); 


PV. ag 
a amas II een 7 
- SWAKE 


where g@ is the velocity potential in the far field. 
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In the case of a pienar wing the vorticity trace in the Trefftz plane can 
be replaced by a slit and equation (32) replaced by; _ 








b/2 
2 
ae Og ) 
DF 5 AoW) ay (Y) dy (33) 
-4/2 
where 
K 
Ag) =a MW (34) 
b/2. 
= Oat = % (35) 
aN Qn (Y - > 
r-b/2 


and K/V_Y) is the total circulation at a given lateral station. Therefore, 


3 (E) 0 


1=1 


and Nj is the mmber of vortices per chord. Also; 


ao 1 1 ; 1 
ay on ; (g } 6 : Yn * Yp. 2 Yn = 


n=1 


where N, is the number of vortices per span. 
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After substituting equations (34), (35), (36), and (37) into equation (33) 
and letting N = Nj x N,; 


N N 
1 > LH 


j=l Kel K 


Since a = b/2 ikl 


06 


j=l Kel °K 7 on Os 


which 1s the same as equation (31). 


The far field calculation of induced drag for a complete configuration 
composed of lifting bodies and thick lifting panels is done by representing 
the wake, from all of the bodies and panels, by an equivalent horseshoe vortex 
system where the bound segment of the horseshoe vortex is tangent to the trace 
of the wake in the Trefftz plane and the trailing legs are in the free stream 

direction. The section drag associated with the equivalent system, which in 
- general is not equal to the actual configuration section induced drag is 
given by; 


q, = p Ws xX r; AS; (40) 
where 
W.= Vv. Dower (41) 
J J J 
and 
cid A AL 
rel @mJ+T XK) (42) 
FIN 


where Ns and T;. are components of the unit vector tangent to the trace of 
J 


the wake at the jth section. 
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T ) il 1 
wo) (poo oa" (38) 
ro) (, Tey te hy 7 } 


NN 
- AY | he a (eee ee 
“y= amb > 2 G} (7 ( +y ¥ ye" ae 


| 





LD 


v,2 > K (gon, +k Ny (43) 
= a + = 
j ~ T . ‘ Yi. uP k k 
- Ry r . . (44) 
Va =f Ie V + ' 
J T jk zy " k Y. | 
where Ty and Ts are the components of the unit vector tangent to the 


k k 
trace of the wake at Re station. Also, Ny and Nz, are components 


of the unit vector normal to the trace of the wake at the Kt" section. N is 
the total number of sections along the trace of the wake for the complete 


configuration. 





“ik “5k 
E. = a _———————— (45) 
V 2 . i 2 De tes 2 
3k Zs + (Fi, + FAS.) Zi, + (Fy 7 FASy) 
Yin As y.. +24 
jk k jk * 
E = - (46) 
W 2? l i eo 1 y, 
_i =A 
MG ha Oe oe ‘i Ge 
where 
Yo (0. 2.) Nee UZ, SF | 47 
Zig. es (Ys aX) Te: 27, -7Y TT. 48 


and the indices j and k refer to the section being influenced and the 
influencing section, respectively. 


Sil 


therefore; 


= p(V..T - W0T. ) AS, (49) 
4; erm Ti,” j 


the total configuration induced drag is then given by; 


NN 
= fe 3 
Dae Lu >, rt, (Ey E + oe aie 


(50) 


therefore; 


> 
=) 
fi 
aE 
M-= 
wie 
a 
< 
s | 
“eee” 
, 
—, 
ca 
eH 
ee 
8 
“——_ 
tri 
<! 
* 
-—< 
~*~ 
+ 
tri 
i 
rs! 
me, 
Saae| 
m™ 


(51) 


=r =) Ey AS; (52) 


ARS » (5) (Z)(—. an, 
1 4b oe) ke] \Y/, Vv ; Ys + Vy Yeu 7 yi J (53) 


which is the same as equations (31) and (30). This demonstrates that the same 
induced drag is obtained whether the equality of work and kinetic-energy 
increment, the equivalent far field horseshoe system, or the near field horse- 
shoe vortices with the Kutta-Joukowsky theorem is used. 
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COMPUTER PROGRAM LISTING 
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